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Konvergenzuntersuchungen und Fehlerabschdtzungen 
fiir ein verallgemeinertes Iterationsverfabren 


JOCHEN W. SCHMIDT 


Vorgelegt von L. COLLATZ 


§ 1. Einleitung 
Zur Lésung der Fixpunktgleichung 
(4.4) DG Le 
in einem abstrakten Raum £ kann man ein verallgemeinertes Iterationsverfahren 
benutzen, das sich mit Hilfe einer Folge von Naherungsoperatoren T,, beschreiben 
1aBt. Die Elemente x, einer Naherungsfolge fiir einen Fixpunkt von T bestimmen 
sich ihrerseits, ausgehend von einem %)€£E, als Fixpunkte der Operatoren T,: 


(4.2) in Ba N13 (ANGE 
(s. J. WarcA [12]). Der Operator T, kann von den vorhergehenden Elementen 
%g, +++, X,—1 abhangig sein, womit er sich erst nach Durchfiihrung des (m — 1)-ten 


Iterationsschrittes benennen laBt. Wahrend T nichtlinear sein kann, wird man 
die T,, oft als linear oder zumindestens einfacher als T wahlen. 

Die Willkiir bei der Wahl der ZJ;, wird eingeschrankt durch eine Approxi- 
mationsbedingung, 

(1.3) lima | Pa ald; 

die sich fiir die Konvergenz von (1.2) als hinreichend und notwendig erweist. 
_ Sie wird von einer Reihe der bekannten iterativen Naherungsverfahren, die sich 
durch passende Festlegung der T,, dem Verfahren (1.2) unterordnen lassen, erfiillt. 

Es soll nun im folgenden gezeigt werden, daB man fiir (1.2) — ahnlich wie 
beim gewohnlichen Iterationsverfahren (s. J. WEISSINGER [13], L. Cotiatz [2], 
J. ScHRODER [9], [10]) — zu Satzen kommen kann, welche die Konvergenz 
sichern, eine Fehlerabschatzung erméglichen und meistens auch Existenz- und 
Eindeutigkeitsaussagen, enthalten. 

Die funktionalanalytischen Hilfsmittel, dargestellt im §2, sind in dieser Zu- 
sammenstellung auch fiir sich von Interesse. Bei Verwendung eines archimedi- 
schen Vektorverbandes N, dessen Elemente zur Beschreibung einer Norm in der 
Grundmenge £ dienen, und der gleichmaBigen Konvergenz in N ergibt sich in 
sehr einfacher Weise die hier benotigte Stetigkeit von inversen Operatoren der 
Gestalt (I — P)+, wobei P ein linearer, positiver Operator in N sein soll. Auch 


Auszug aus der Dissertation des Verfassers, TH Dresden, Marz 1959. Es sei an 
dieser Stelle den Referenten Prof. H. Hernricu, Dresden, und Prof. K. MAruun, 
GieBen, sowie Prof. P. H. Mtrier, Dresden, fiir das anhaltende Interesse an dieser 
Arbeit herzlich gedankt. 
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aus weiteren Griinden ist die gleichmaBige Konvergenz im Falle einer allgemeineren 
Norm als Ersatz fiir die Konvergenz von Zahlenfolgen, die doch bei der Zahlen- 
normierung herangezogen wird, sehr vorteilhaft, wie z. B. der Beweis zu Satz 1 zeigt. 

Der §3 enthalt Satze iiber das verallgemeinerte Iterationsverfahren, wahrend 
im §4 einige Verfahren angefiihrt werden, die in diesem enthalten sind. Der 
§5 schlieBlich bringt numerische Beispiele. 


§ 2. Funktionalanalytische Hilfsmittel 


Es werde zundchst die Menge N beschrieben, deren Elemente zur Festlegung 
einer Norm in der Grundmenge E verwendet werden. Ein halbgeordnetes Linear- 
system N ist ein Linearsystem tiber dem Kérper der reellen Zahlen, in dem eine 
Halbordnung erklart ist, welche also die Ordnungsaxiome erfiillt (s. G. BIRK- 
HOFF [1]) und monoton beziiglich der Addition und der Multiplikation mit reellen 
Zahlen ist, d.h. ausa<b, cE N,A>0, A reelle Zahl, folgt at+cSb+c, unddaSdb. 

Existiert fiir jedes a€ N das Element sup(a,0), so heiBt N Vektorverband 
oder Rieszscher Raum. Dort sind zu je zwei Elementen stets das Supremum 
und das Infimum vorhanden. Man kann somit in Vektorverbanden den Betrag 
eines Elementes einfithren, |a| =sup(a, — a); es ist also |a|€N, und es gilt wie 
bei den reellen Zahlen 

a) |a|=0 genau dann, wenn a=0 
(rep Aahee Altay A reelle Zahl 

c) [a+] Jal + [5]. 
Setzt man weiter a*=sup(a, 0), a =sup(—a, 0), so erhalt man a*=0, a =O 
und a=a*—a, |a|=at+a-. 

Ein Vektorverband heiBe archimedisch, wenn na<b fiir a,b€N und jede 
natiirliche Zahl die Ungleichung a0 impliziert. Als Konvergenzbegriff in N 
wollen wir die gleichmaBige Konvergenz verwenden. Eine Folge {a,,} von Ele- 


menten aus N konvergiere gleichmaBig gegen aC N, wenn es ein }€ N und eine 
monoton fallende Nullfolge {/,,} von reellen Zahlen gibt mit 


|a,—a| <A,b, Wize Aly Braise 


Hierfiir schreiben wir kurz a=lim a, (s. G. BrrKuorrF [7]). Wie man leicht tiber- 
priift, gelten folgende Rechenregeln 


a) aus a,=—a folgt lima, =a, 


b) falls lim a, =a ist, so auch lim a,,=a fiir jede Teilfolge {z,,} der natiir- 
lichen Zahlen, 


(2.2 c) gilt lima,=a, limb, =), so existiert lim(a,,+,)=a-+, 


d) aus lima,—=a, lim«,=a« (a,, « reelle Zahlen) folgt die Existenz von 


lim ea, oa; 


e) aus lima,=0, 0Xb,<a,, ergibt sich lim b,=0, 
f) existiert lim a,=a und ist a,=0, so ist auch a=0*. 


n= 


* Einen abstrakten Konvergenzbegriff mit diesen Eigenschaften benutzt J. SCHRO- 
DER [9]. Aus (2.2) allein folgt jedoch nicht, wie es hier erforderlich ist, die Stetigkeit 
von linearen, positiven Operatoren in N. Auch die von L. W. Kantorowi1rscH [6] 
verwendete Ordnungskonvergenz hat diese Eigenschaft nicht. 
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Zum Beispiel ergibt sich c) aus der Ungleichung 
\(@, + 6,) — (a+ 6)|S|a,—a| +|b,—b| SA, c+yu,4d< (A,+4u,) sup(c, d), 


wahrend beim Nachweis von f) benutzt wird, daB N archimedisch ist. Es gilt 
wegen |a, —a| </,c, 4,20 und weiter, da {/,,} eine reelle Nullfolge ist, 


eg SA 
m 


also m(—a)<c fiir jede natiirliche Zahl m. Das bedeutet aber a=0. Ebenso 
bestatigt man die Eindeutigkeit des Grenzwertes. 

Als Beispiele fiir archimedische Vektorverbande sollen, hier der k-dimensionale 
Vektorraum Rk, und der Raum der stetigen Funktionen C erwahnt werden. Die 
za a€C gehérige Funktion |@| wird durch |a| (¢)=|a(é)| erklart, und die ab- 
strakte gleichmaBige Konvergenz bedeutet in C die gleichmaBige Konvergenz im 
Sinne der Analysis. 

Um Lipschitzbedingungen fiir Operatoren in F angeben zu kénnen, bedienen 
wir uns linearer, positiver Operatoren in N, sog. Lipschitz-Operatoren. P ist 
ein solcher Operator, wenn gilt 


a) PAa+ub)=AP aur od, 
b) aus a2 0 folgt PaZo. 
Es sei hervorgehoben, daB P dann stetig ist, d.h. aus lim a, =a folgt lim Pa, = Pa. 


Zum Beweis wird die Ungleichung |Pa|<P|a| bendtigt, die sich folgender- 
maBen ergibt: 


(oN | Paw a |S | Pad --|Pa |e Pat Pa = Pal . 
Ist nun |a, —a| SA,,c, so wird damit 


|Pa,— Pa| =|P(a,—a)|<Pla,—a|SA,Pc; 
also ist P stetig. 


Konvergiert die zu P gehérige geometrische Reihe 
(2.3) DuPis (an, die Folge MATES ¢}) 
»=0 v=0 


fiir jedes aC N gleichmaBig, so existiert der inverse Operator (J — P)+ und ist 
an der Stelle a gleich der Reihe (2.3) (s. J. SCHRODER [9]). Insbesondere ist 
(I — P) + linear und positiv, da mit P auch alle P” positiv sind und wegen (2.2). 
Also ist (I — P)+ nach dem Vorhergehenden stetig. Es sei noch erwahnt, daB 
jedoch IJ — P nicht positiv zu sein braucht. 

Die Grundmenge £, deren Elemente mit 1, y, z,... bezeichnet werden 
sollen, sei ein N-normiertes Linearsystem beziiglich der reellen oder komplexen 
Zahlen. Jedem x€E wird also ein Element ||%||€ N zugeordnet, so daB die 
iiblichen Normenaxiome erfiillt sind. Ist E selbst ein Vektorverband, so kann 
man E durch die Festsetzung ||«||=|x| normieren, wie (2.1) zeigt. 
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Die Konvergenz in E fassen wir im Sinne der Normenkonvergenz aut, d.h. 
limx,=x bedeutet lim || x, —x||=0. & sei beztiglich dieser Konvergenz voll- 
standig, also konvergiert jede Cauchy-Folge aus E gegen ein Element aus E£. 
Die Konvergenzsatze des §3 gelten auch, wenn man in E mit Hilfe von N statt 
der Norm eine Metrik einfiihrt, die den folgenden Axiomen geniigt (s. H. Enr- 


MANN [4]) 
a) r(x, y) =0 genau dann, wenn *= y 
b) r(x, y) S7r(x, 2) +7(y, 2) 
c) r(x, v) =r(x +2, y+2) fiir jedes z€ E. 


Hieraus folgert man, daB 


v(—*,0)=7(%,0), 7(%¥+ 4,0) S7(x, 0) +7(y, 0) 
ist. 
§ 3. Konvergenzsatze 

Wir beginnen mit einem Satz, der bis auf einige fiir die Fehlerabschatzung 
vorteilhafte Umstellungen im Falle der Normierung durch die reellen Zahlen 
von J. WARGA [12] bewiesen wurde. Wir iibertragen ihn hier auf die allgemeinere 
Normierung. 

Satz 1. Es set N ein archimedischer Vektorverband, E ein N-normierter voll- 
stindiger Raum, T, T,, (n=1, 2, ...) Operatoren in E mit folgenden Eigenschaften: 


a) die Gleichung (1.1) tst durch x€ E losbar, 

b) |Tv —T,2z||SM||y—z]| fir vy, 2€E, n=1, 2,...; M Lipschitz-Operator 
mit fiir jedes aC N gleichmaBig konvergenter geometrischer Rethe (2.3), 

c) \\(T —T,) y-—(T—T,)z|| SK ||vy —2|| fir y, z€E, n=1, 2,...; K Lipschite- 
Operator, 

d) die geometrische Rethe zu P=(I—M)1K konvergiert fiir jedes aE N gleich- 
maprg. 

Behauptung. Die Gleichung (1.1) ist in E eindeutig lésbar durch x, und fiir 
die Konvergenz von (1.2) gegen x bet beliebigem Anfangselement x)€E ist (1.3) 
hinreichend und notwendig. Es gilt die Fehlerabschitzung 


(3-1) ||~—ml| SU —P)> Pil — xoll+ (2 —P)*(L —M)2||F % — T x9||*. 
Beweis. Um zu zeigen, daB (1.3) hinreichend ist, fithren wir folgende Ab- 
schatzung durch: 
Il%n4a— #1] S[Zn42 2nd — T4||S [Tha Mga — Tha || + 
+ |[(P= Th) pS a) ll a T x,|| 
SM |l¥u41—*l) + K|lx— 211+ [Toe — 7 Hall 


* Zu besseren Fehlerschranken kann man kommen, wenn man die Operatoren 
nicht durch lineare Majoranten (Lipschitz-Bedingungen), sondern durch nichtlineare 
abschatzt (s. J. ScHRODER [J1]). Hier ware es sinnvoll, fiir T —TZ, nichtlineare 
Majoranten heranzuziehen, was aber nicht weiter ausgefiihrt werden soll. 
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oder da (J — M)* positiv ist, 
l|%n44— x|| = Pao || Ty, 


wobei {a,} mit a,=(I—M)+]|T,,.%,—T%,|| wegen (1.3) und der Stetigkeit 
von (I—M)*+ eine Nullfolge ist. Mit der Abkiirzung n’ gleich gréBte ganze 
Zahl <n/2 folgt hieraus weiter 


(3.2) |] Hae a x || = ees || %o— x|| =f > pr a > Ra Sha 
v=n' +1 v=0 
Wir zeigen jetzt, daB jeder der drei Summanden gegen Null aus WN strebt. 
Wegen d) ist 
n+1 n 
lim P"* || % — x|| =lim >)-P?|| x) — x|| — lim >) P” || x — x]| =0. 
v=0 v=0 
Es war {a,} eine Nullfolge, also gibt es ein b€ N und eine monoton gegen Null 
fallende Folge von reellen Zahlen {/,,} mit 
Gy | a aap = AO, W=20, Asc 


A sei eine Schranke fiir alle A,. Somit erhalt man fiir den zweiten Summanden 


von (3.2) 
6) 


IA 


YEP ae Pr b= 1a Pe, 
=0 v=0 


v=n'+1 v=n/+1 » 
d.h., wegen d) strebt auch dieser Term gegen Null. Weiterhin ist 4,5 1,05 1,6 
fiir v=n’, also folgt 
Ole ae fe A FOS pa OSA ae 
v=0 v=0 v=0 


Damit ist insgesamt gezeigt, daB (1.3) fiir die Konvergenz von {x,} gegen x 
hinreicht. Gleichzeitig sieht man, daB (4.1) nur eine Lésung in FE haben kann, 
denn sonst miiBte {x,,} auch gegen diese weiteren konvergieren. 


Die Notwendigkeit von (1.3) ergibt sich aus 
aia ee, || | Lote a eel alld ate lige 
SM ||%n 41 — %nl] + [na — 41] + Bll n— ¥|] + Mn — % 


d 


wenn man nur *=lim%, beachtet. Die Fehlerabschatzung (3.1) erhalt man nach 
geringer Umformung aus (3.2) ftir ~=0. 

Fiir praktische Belange ist es meistens wichtig, die Lipschitz-Bedingungen 
an J, und J — T, nicht auf dem ganzen Raum £, sondern nur auf einer méglichst 
kleinen Teilmenge von EF zu fordern, wie es etwa fiir das gewohnliche Iterations- 
verfahren in [2], [9] durchgefiihrt ist. Fiir zwei Falle soll dies beriicksichtigt 
werden, zunichst fiir solche Verfahren, die einer scharferen Approximations- 
bedingung geniigen, namlich 
(3.3) Eup he, n=0,1,.... 


Die Existenz eines Fixpunktes von T braucht jetzt nicht mehr vorausgesetzt 


zu werden. 
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Satz 2. Es sei E ein N-normierter Raum, F eine geeignete vollstandige Terl- 
menge von E, auf der die Operatoren T, T,, (n=1, 2,...) erklart sind, wober die 
Bildbereiche der T, (n=(, 2, ...) wieder in F legen sollen. Weiter gelte neben (3.3): 
—2z|| fir y,z€F, n=1,2,...; M Lipschitz-Operator, 

b) \l(T —T,) y —(T—T,) 2|| SK ||y —2|| fur y, 2€F, n=1, 2,...; K Lipschate- 
Operator, 

c) die geometrischen Reihen (2.3) zu M und P=(I—M)*K konvergieren fir 
jedes AEN gleichmafpig, 


— x, || </f} mit 


d) ausgehend von einem geeigneten x»CF set D={y€E; 
f=(I—P)+P||x,— %|| in F enthalten. 

Behauptung. Die Gleichung (1.1) ist in F eindeutig losbar durch x, und das 
Verfahren (1.2) konvergiert gegen diese Losung. Auferdem legt x im D, d.h. 


(3.4) le — «ll Sf=(2—P)7 Pll — |. 


Beweis. Auf Grund der Voraussetzungen sind die Gleichungen (1.2) ein- 
deutig lésbar in F (s. auch L. Coriarz [2]). Man erhalt 


ca — Vall Seca na — a elle ee ee 
+\\(2—T) %— (T= B) all +L. a — Ta 
SM || %,41— %,|| + Kl 4, — %n-1| 
ee Pl ee ee al 


(3.5) 


woraus sich sofort 
Ne cy | Cd tg ie 20 3 |e aes | 


ergibt. Wegen der Stetigkeit von (J — P) + und lim P*|| x, — x9|| =0 folgt hieraus, 
daB {x,,\ eine Cauchy-Folge ist, deren Elemente auBerdem alle in D liegen. (Man 
setze m=1 und dann k,=n.) Das Grenzelement «=limx,, ist ebenfalls aus D, 
wie aus der Ungleichung 


I|* — sal S lle — *al| + @ — P)> Play — | 


durch Grenziibergang hervorgeht. Also ist (3.4) bestatigt. Die Fixpunkt- 
gleichung (1.1) wird von x erfiillt, denn es ist 


O<||Tx—x||S||Tx—T, x,||+]|x,—-|| 


Syn T,,) ea T,) % Xn— ur ee qT, %q|| + || %, ro Ea Tn il te 
=K||x—~,-4||+M||x %,|| + ll % 4. 


Ware schlieBlich y=+ x eine weitere Lésung von (1.1) in F, so ergibt sich fol- 
gendermaBen ein Widerspruch: 


|x = yl] =P x — Tyl| S2K lx — al] + Ke — y|| 4 M29, 

also wegen x=lim x,, 
| 
les ey || See 


woraus wegen c) ||x — y|| =0 folgt. 
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Anmerkung. Sind die Operatoren TJ, (n=1, 2,...) fast linear, d.h. T,y= 
W,, v+2, fir yCE mit linearen Operatoren W, in E, z,€E, so kann man die 
Existenz der Fixpunkte x, von J, in F sichern, wenn man, was bei solchen 
Operatoren vielfach keine Einschrankung bedeutet, a) aus Satz 2 fiir alle y,zC€ E 
verlangt und dafiir die Forderung, daB die J, F in sich abbilden, fortlaBt. 

Zunachst wei8 man namlich nach dem Fixpunktsatz tiber das gewohnliche 
Iterationsverfahren (s. J. SCHRGDER [9]) damit, daB T,, in E genau einen Fix- 
punkt x, hat. Um den Anschlu8 an den vorherigen Beweis zu gewinnen, muB 
man aber bestatigen, daB die x, aus F sind, was hier induktiv geschehen soll. 
Es mégen also %,, yn, aus F sein. Dann gilt in gleicher Weise (3.5), da jetzt 
T,,4, tiberall lipschitzbeschrankt ist und T — T,, nur Argumente aus F hat. Daraus 
folgt wie vorhin x, .,€D, also auch %,,€F. 

Nachdem wir speziell Verfahren mit der Approximationsbedingung (3.3) be- 
handelt haben, wollen wir uns jetzt einem weiteren Sonderfall von (4.2) zuwenden, 
namlich dem Verfahren, das H. EHRMANN [6] ausfiihrlich untersucht hat. Bei 
Benutzung eines anderen Beweisgedankens und etwas abgewandelter, zum Teil 
geringerer Voraussetzungen kann man fiir das Verfahren 


(3.6) a ee Lys n = 4p Dy, eee 
folgenden Satz angeben. . 
Satz 3. Es mége E ein N-normierter, vollstandiger Raum sein, F eine geeignete 


Teilmenge von E, auf der die Operatoren T, S,, (n=0,1,...) erklart sein sollen, 
und x,CF. Weiter gelte 

a) ||Ty—Tz||SP|ly—2zl| fir y, zEF; P Lipschitz-Operator mit gleichmapig 
konvergenter geometrischer Rethe (2.3), 

b) es bestehe die Approximationsbedingung 


(3.7) lim:||S,ja, Pow 09 
und auperdem ser iho ilpoereanaa ery nee 
c) schlieBlich moge D={yEE; ||y — || S(Z— P) +P ||% — xo|| +2 — P) 43} 
in F enthalten sein. 
Behauptung. Das Verfahren (3.6) konvergiert gegen den einzigen Fixpunkt x 
von T in F, und eine Fehlerabschatzung erméglicht die Ungleichung 


(3.8) |x — %|| S (J — P)2 Pla — |] + FZ — P) AIT % — So oll. 
Mit der Festsetzung 
(3.9) LV = Spo pi Vers m=A,2,... 


1aBt sich (3.6) als Verfahren (1.2) auffassen. Es wird dann M=0, K=P (Lip- 
schitz-Bedingung fiir T); die Approximationsbedingung (1.3) geht in (3.7) und 
die Abschatzung (3.1) in (3.8) tiber. 

Kommen wir jetzt zum Beweis, der im wesentlichen durch Vergleich mit 
dem gewohnlichen Iterationsverfahren 


(3.10) “a= LU 3, =A, 2,.2-, Up=%o 


gefiihrt wird. 
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Es sei D={yEE; lly —|| S(— P)> P |\u, — %9||} und z€D, dann ist 


\|2—44 || S ||2—ea || +142 SP) 4 P (4-44 || + || 2 —20l)) + || —a|| 
(3.44) = (I — P)4P|| x, — %|| +1 — P)4||T % — So *ol| 
<(I— P)4P||x,— || +2(0—P)0, 


d.h. z€ D oder DCDCF. Nach dem Fixpunktsatz tiber (3.10) (s. J. SCHRODER [9]) 
hat T damit genau einen Fixpunkt x=lim wu, in F, und auBerdem sind wu, €D, 
VE Ng Oe BO0 

Man kann weiter induktiv nachweisen, daB x,¢D, n=1, 2,...; denn wegen 
x, (Induktionsannahme), u,€F fiir ySn gilt 


|| #041 — Mga] S [| 22 — T Mal] + [l2nl] S Pll %n — Mall + [2] 
SP*|[2o|| +P" all +--+ Pll4n—all + [lznll Ss —P) 74. 
Da pee D) folgt weiter nach (3.11) 


(3.42) 


|| %n4a— %1|| S|] %e4a— Mn-eal] + |] en ta— %||S (I—P)7 P||% — %|| + 2U—P)*6, 
also ist auch x, ,,€D. 

Wie im Satz 1 fiir die rechte Seite von (3.2) beweist man mit Hilfe von (3.12), 
daB wegen lim||z,||=0 auch lim||x, —w,||=0 ist, insgesamt lim x,=limu, =. 


Die Fehlerabschatzung (3.8) schlieBlich ergibt sich wegen u,€D und (3.11) 
durch Grenziibergang aus 


le—aall Sll2—xll + [eee + E-P)A Pll 4 —al] + IP YAIT x0 Sol. 
Damit ist der Beweis vollstandig erbracht. 


Wir wollen jetzt Bemerkungen iiber die Durchfiihrung der Fehlerabschatzung 
sowie tiber die wiederholt verwendeten geometrischen Reihen von Lipschitz- 
Operatoren anschlieBen. Zur Berechnung oder weiteren Abschatzung der GréBe 
f=(I—P)+P||x,—%o|| hat J. ScHRODER [9], [70] bei bekanntem P eine Reihe von 
Methoden angefiihrt, von denen wir die in unseren Beispielen benutzten erwahnen 
wollen. Zunadchst kann man f im Falle P=(J—M)+K aus der linearen Glei- 
chungen 2. Art 
(3.43) f—(M+K)f=K||n— xl 


berechnen. Gibt es weiter ein b€N, b=0, und reelle Zahlen « und y, so daB 
die Ungleichungen 


(3.14) || x1 — xo|| Sad, 
(3.15) Pb<yb 


mit y<1 bestehen, so erhalt man die fertigen, gegentiber (3.13) gréberen Ab- 
schatzungen 


(3.16) ee Nees a 
19 eh Sy 


Dabei ermittelt man das Element e= Pb zweckmaBig aus der linearen Gleichung 


(3.17) e—-Me=—Kob. 
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Kommen wir jetzt zu den geometrischen Reihen der Form (2.3) zuriick. Thre 
gleichmaBige Konvergenz fiir jedes a€ N wird in den Satzen dieses Paragraphen 
fiir die Operatoren M und P gefordert. Hierfiir soll eine einfache hinreichende 
Bedingung angegeben werden. 


Ist P ein Lipschitz-Operator und }=0 ein Element aus N, so daB neben 
(3.15) mit y<1 weiterhin fiir jedes aC N 


(3.18) |a| <Bb 


mit einer Zahl B=£(a, b) gilt, so konvergiert die geometrische Reihe zu P fiir 
jedes aE€N gleichmaBig, wenn wir noch zusadtzlich die wenig einschrankende 
Forderung stellen, daB auch N vollsténdig ist. Der Beweis ergibt sich wegen 
der Ungleichung | P’a| <f y"d aus 


kn 
De 


denn hiermit ist gezeigt, da8 die Partialsummen eine Cauchy-Folge bilden. 


es 


kn kn 
Se er yb, 


Wir werden weiterhin folgendes nachweisen. Ist fiir obiges } 
(3.19) MO=00,> KOS 6b, O+ex<i4 


und gilt (3.18), so konvergieren die geometrischen Reihen zu M und P= (J —M)1K 


fiir beliebiges a€ N gleichmaBig, und in (3.15), (3.16) kann man y= ; e 5 setzen. 


Wegen 6<1 gilt nach dem vorhergehenden die Behauptung fiir @. Um sie 
auch fiir P zu bestatigen, miissen wir (3.15) tiberpriifen. Es ist 


Pose(I—M)1b=cb+elim> Mbseb+elim> 0b 
y=1 y=1 
SOM yyema -f oi 
b a yb 


ate} 
anges 


§ 4. Iterationsverfahren der Gestalt (1.2) 


Es werden hier einige bekannte iterative Naherungsverfahren zur Lésung 
der Fixpunktaufgabe genannt, die sich zwanglos in die Vorschrift (1.2) einfiigen 
lassen. Man kann noch eine Reihe weiterer Verfahren anfiihren oder auch kon- 
struieren, wobei man im wesentlichen Bedingung (1.3) zu beachten hat. 


4. Das gewohnliche Iterationsverfahren x, — Tx, _, erfordert die Festsetzung 
Tigi 2%, 27 Oflitway GE; 


die Operatoren T,, sind also konstant zu wahlen. Die Bedingung (3.3) ist erfiillt. 


2. Beim Verfahren in Einzelschritten, falls es ttberhaupt auf ZT anwendbar 
ist und die J, dann sinngemaB erklart werden, besteht die Beziehung (3.3). 


3. Das gewohnliche Newtonsche Verfahren fiir die Fixpunktgleichung lautet 


44 = if Xy—-1 =f Toe (Xn Fay %n—1) , 
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wobei die Ableitung von 7 etwa im Fréchetschen Sinne zu verstehen ist. Hier 
bestimmt sich der Naherungsoperator T;, also zu 


LV HL Kyat tl Gig eh Ohne eek OY ae: 


Die Gleichungen (3.3) sind erfiillt, ebenso auch fiir folgende Vereinfachungen 
des Newtonschen Verfahrens: 


Tye Ky oy he Ly oe Gye alg Vee 
bzw. wenn man die Ableitungen durch lineare Operatoren S, in EF ersetzt, 


Riyss TAS AV 44s yEE, 


oder speziell durch Zahlen «,,, 
Di a He eae yeE. 


4. Eine Reihe von Verfahren héherer Ordnung (s. z.B. R. Lupwie [8]) ge- 
ntigen ebenfalls (3.3). 

5. Hat T die Gestalt T=R-+S, wo R linear oder jedoch einfacher als S 
sein soll, so ist oftmals die Iteration x,=—Rvx,+S-x,_, zweckmaBig. Bei der 
Festsetzung 

Ly=Ry+S%1, ~~ yeE, 


gilt die Approximationsforderung (3.3). 

6. Das von H. EHRMANN [5] untersuchte Verfahren mit veranderlichen Ope- 
ratoren (3.6) erfiillt nur die Approximationsbedingung (1.3), wenn man zusatz- 
lich (3.7) fordert. Dies entnimmt man der Darstellung (3.9). Ebenso kann man 
in der Iterationsvorschrift 5. die Operatoren R, S durch veranderliche R,, S,, 
ersetzen und hat dann Bedingungen der Art (3.7) zu verlangen, damit (1.3) 
erfiillt ist. 

§ 5. Beispiele 
1. Das Newtonsche Verfahren fiir Integralgleichungen 2. Art 
Die nichtlineare Integralgleichung 


(5.4) x(S) =e (s) +J6(s t, x(t)) dt 


soll durch das vereinfachte Newtonsche Verfahren gelést werden: 
1 1 
(5.2) x, (s) =g(s) + [G(s t, %,r()) dt f 5 (s, t, %9(E)) (u(t) — %y-1(t)) dt. 
0 0 


Dabei seien g(s) fiir OSs<1 sowie G(s, t, x) und 26 b, %) TOPO Ss, ete 
- =5,15 


u(s)SxSv(s) erklart und stetig. Fa8t man alle stetigen Funktionen, die fiir 
OSsS1 zwischen u(s) und v(s) liegen, zur Menge F zusammen, so vermitteln 
die Operatoren 
1 
Py(s)=s(s) + SG(sty()) dt, yer 


0 
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| bzw. A 
ir cio + [4 Cea Natt f Fe (s,t, xo(t)) (y(t) —%-1()) dt, yeC 
| 0 


| Abbildungen von F bzw. C in C. Mit Hilfe der stetigen Funktionen 


& (5,4, xo(2)) | 


| kann man Lipschitz-Bedingungen an die Operatoren J,, und T — J, formulieren. 


(5.3) M(s,4) = [2 (s.4,x0(0)], K(s,t) =max 


0G 
ay Re hire (s, t,x) — 


A. Im Falle N=R, und ||-|| = max 6H mit einer festen, fiir 0<s<1 


stetigen und positiven Gewichtsfunktion w erhalt man 


f && (s, 4, x9(t)) (9 () — 2 (0) dt 


0 


| |Z.7 —T,2|| = max 


0Ss<1 w(s) 


al 
[mee. 0 POZO a < Mly—z fiir y,z€C, 


(FT, Tak || fei) tiny, 26F 
mit 


(5.4) at, 


0<s<1 w(s 0<s<1 w(s 


Erweisen sich nun M<1, ee und sind alle stetigen Funktionen y, 
die der Ungleichung 


P | %1 (t) = %o (2)| 
a & Uke BEd Sis 
(5.5) [Ue 19 ee ay ee el 
gentigen, in F enthalten, so konvergiert die gemaB (5.2) ermittelte Funktionen- 
folge {x,,(s)} gleichmaBig gegen die in F eindeutige Lésung x(s) von (5.1), welche 
gleichzeitig in der Menge (5.5) liegt (Fehlerabschatzung). Dies folgt aus Satz 2 
und der Anmerkung dazu. 


B. Es sei N=C und ||x||(s)=|x(s)| fir OSs<1. Bei der Ordnung, ab 
falls a(s)<b(s) fir OSs<1, ist C ein archimedischer Vektorverband. Wegen 
(5.3) erhalt man 


1 
|Z, v—T,2|| (s) SS M(s, 4) |y(Q—z (yl dt= fs, 0) |]y— z|\|(@) dé fir OS s=1, 
0 


also ||TZ,, y — T,, 2|| SM||y —2||, wobei M der durch 
M a(s =f Ms, a(idt, acéeC 


erklarte lineare, positive Operator in C ist. Entsprechend definiert man K. 
Um Satz 2 anwenden zu kénnen, bendtigen wir z.B. die Funktion / aus (3.13) 
oder auch e aus (3.17), d.h., wir haben die lineare Integralgleichung 
i 


(5.6) f(s) — f (M(s, 4) + K(s,4)) F@) og) t) | (t) — x9(t)| at 


0 
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bzw. } , 
(5.7) e(s) — f M(s, 2) aly ai Joe t) b(t) dt 


zu lésen. Die weiteren sich aus Satz 2 und der Anmerkung ergebenden Bedin- 
gungen wollen wir hier tibergehen. 


Numerisches Beispiel. In der Integralgleichung 
1 
x(s) =1—Is+s?+ fstarctg x(t) dt 
0 


ist die Konstante 1=0,485 38 so bestimmt, daB x=1-+s? eine Lésung ist (ent- 
nommen aus L. W. KAnTorowIrTscu [7]). Das Newtonsche Verfahren (5.2) fiihrt 
mit der Naherung x,=1,5 auf eine Integralgleichung mit ausgeartetem Kern, 
deren Lésung x,=1+0,00671s-+s? ist. Der Fehler von x, betragt also 0,00671s. 

A. Im Falle N=R, erhalten wir mit f= {ye C.0,92=¥ (s) 52,5 ir VSs= 1} 
aes M(s, t) =0,30770st,  K(s,t)=0,22879 st 
und damit fiir w(s)=1 in (5.4) die Fehlerabschatzung 

|x(s) — %(s)| $0,079, OSsS1, 


wie man aus (5.5) entnimmt. Alle weiteren Voraussetzungen iiberpriift man 
leicht. Insbesondere ist die durch die letzte Ungleichung beschriebene Menge 
ganz in F gelegen. 

Es sei noch vermerkt, daB es hier nicht méglich ist, durch Wahl einer ge- 
eigneten Gewichtsfunktion w den Fehler als zu s proportional anzugeben. Dies 
gelingt erst bei der Verwendung von N=C. 

B. Wenn wir N=C setzen, kénnen wir F etwas verkleinern, F= 
{yEC; 0,999< y(s)<2,5 fiir 0<s<1}. Nach (5.3) ergeben sich dann 


M(s,t)=0,30770st,  K(s,t)=0,19281 st. 


Mit 5(s) =4 sieht man sofort, daB in (3.19) 6+e<1 ist, also konvergieren die 
zu M und P gehorigen geometrischen Reihen fiir jede stetige Funktion gleich- 
maBig. GemaS (3.16) erhalt man nach Auflésung der Integralgleichung (5.7) 
— hier mit ausgeartetem Kern — die Abschatzung 


|x(s) — %(s)|S0,061s, O<s<1, 


und mit etwas mehr Miihe aus (5.6) — jetzt ebenfalls eine Integralgleichung mit 
ausgeartetem Kern — die wesentlich bessere Fehlerschranke 


|x(s) — %1(s)|S0,029s, OSs<1. 


Da die durch diese Ungleichungen beschriebenen Mengen ganz in F liegen, 
konvergiert also das Verfahren (5.2) gleichmaBig gegen eine Lésung x(s) der 
behandelten Integralgleichung, die in F eindeutig ist, und es bestehen die eben 
notierten Fehlerabschatzungen. 

Ohne Miihe findet man auch, daB die Integralgleichung in ganz C eindeutig 
lésbar ist; denn fiir F=C werden die Zahlen M und P=K/1—M aus (5.4) fiir 
w(s) =1 ebenfalls kleiner als 1. 
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2. Eine Randwertaufgabe 
L. CoLtatz [3] lést die Differentialgleichung der Knickbiegung 
es) 4-1—0, x(+1)=0 
| naherungsweise durch das Iterationsverfahren 
(5.8) — x — 4% =14s?%, 1, Kp (a4) = 0; 


Mit dem Ansatz %)=/ (1—s*) kommt man iiber (5.8) zu 


%(s) = — 1+1 (26 — 1352+ sf) +4—" coss, 
wo / sich bei der Forderung %)(0)=%,(0) zu 1=0,933505 ergibt. Der Fehler 
von %, ist an der Stelle s=0 dann genau gleich 0,001 45, wie man bei L. CoLtatz 
 findet. 
Zu einer bequem durchfiithrbaren Fehlerabschatzung von %,(s) gelangt man 
mun, wenn man die Vorschrift (5.8) als verallgemeinertes Iterationsverfahren 
(1.2) schreibt. Mit der Greenschen Funktion 


en 


— 


eee EDN 4 fire Nees <Y 
G(s, t) = as 


(t+1) A—s) fir t<s<i1 
2 = =——~») 


die zu — x’’=r(s), x(+1)=0 gehort, laBt sich (5.8) darstellen als 
if 1 
i4(8) = RB xq (s) + S%qa(s) =f C(s,2) x (0) dt+ f G(s,t) (1+ 2 xy a0) dt, 
—3! —1 


wobei die Operatoren R bzw. S durch das erste bzw. zweite Integral erklart 
sein sollen. Dann vermitteln die Operatoren T, JT, (n=1, 2,...), die durch die 
Gleichungen 

(5.9) TVR ie Son 2G, Y BRYA SO Kp. EC 


festzulegen sind, Abbildungen von C in sich. 
A. Im Falle N= R, und || «|| = “max Acau) sind 


M = max a dt = max a4 oe 


—1<=s<1 —I<ssl 2 


bzw. 


= max "el (s,t) ?@di= max L (i= s4):= 


- -Isssi —tss=1 12 42 


Lipschitz-Operatoren fiir rake bzw. T —T,, (n=1, 2, ...) fiir alle y, z€C. Da weiter, 
wie aus (5.9) ersichtlich ist, die Approximationsbedingung (3.3) besteht und d) 
aus Satz 2 wegen F =C trivialerweise erfiillt ist, sichert dieser Satz die Konvergenz 
des Verfahrens (5.8) gegen die in C eindeutige Lésung x(s) der zur Differential- 
gleichung der Knickbiegung gehérigen obigen Integralgleichung. Gleichzeitig 
erhalt man, da || x, — %9|| <0,0124 ist, nach (3.4) die Fehlerabschatzung 


|x(s) — x,(s)| S$ 0,00248, =4= 551. 
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B. Zu besseren, keinen gréBeren Aufwand erfordernden Fehlerabschatzungen 
gelangt man, wenn man N=C setzt. Mit den Lipschitz-Operatoren fiir TJ, bzw. 
7—J,, gultig fir alley,.7EC, 


1 1 
Ma(s)=fG(s,t)a(t)dt, Ka(s)=JG(s,t)Pajdt, aEC 
—1 —1 
und b(s) =1 geht (3.17) in die Randwertaufgabe 


an (fh [ee e(+1)=0 


iiber, deren Lésung e=2—s?—1,8508 cos s<0,149 ist. Mit «=0,0124 ermég- 
licht somit (3.16) die Abschatzung 


|~(s) — x(s)| S 0,0147 e(s) [ 0,00219, —15s<1. 


Von der erforderlichen Konvergenz der zu M und P gehorigen geometrischen 
Reihen kann man sich gemaB (3.19) wegen 6=4 und e= ,, iiberzeugen. 

Die Gleichung (3.13) zur Bestimmung der Fehlerschranke / fiihrt auf eine 
Randwertaufgabe, die schwieriger ist als die urspriinglich zu lésende. Wir wollen 


diese Méglichkeit daher nicht weiter verfolgen. 


3. Ein nichtlineares Gleichungssystem 


Zur naherungsweisen Lésung des Gleichungssystems 


Hey — 2%ay + Vey + 24 — 4,25 =0 
16 Xa X{2) — 4% ay Xe) Vea) + 6 — 8 xq) Xe +1=0, 


das von %q)=0,5 und x(2,=1 erfiillt wird, ist das Verfahren 


Xn(1) = 0,5 ae qa) + 0,5 eas + 24 — 2,125 


1 ; —3——— 
OS 6+ 0, 
RE Ee 5 x2 Lijit OS.035 


(5.10) 


Xn (2) = 
gut geeignet. Mit Hilfe von 


0,5 YA) = 
Cree Oh So (eee 24 ce y= (Jen, 


16 V1) V(2) 


bestimmen sich T und T,, (w=1, 2, ...) wie in (5.9), sind jedoch jetzt alle nicht- 
linear. Aus den Anfangswerten %9q)=0,55 und X%9(2)=0,95 ergibt sich iiber 
(5.10) die Naherung 


Xa) = 0,495 O2F X71 (2) = 1,001 56, 


deren Abweichung von der genannten Lésung wir jetzt abschatzen wollen. Da 


fiir (5.10) die Approximationsbedingung (3.3) gilt, kénnen wir hierzu Satz 2 
heranziehen. Definiert man die Teilmenge F durch 


0,48 = X(1) =< 0,55 ) 0,95 = X(2) < 1,03 5) 
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so bilden die J, F in sich ab, und weiter ist x»¢F. Im Normenraum N=R, 


mit der Norm || *||(,=|* |, 7=4, 2, ergeben sich die Lipschitz-Operatoren M 
und K fir J, und T—T,, zu 


Map| 0 
/ 0,55 0) 
ESE AON a | mm | 2 A f 286 0 ia 
16 V&) Vey 16 Vay Ve) 
mals — marl ae ga Omen) (2, 01504 
yCF FEF (10,25 yay(y+O|[-F 0 0,05 5lamOn 3a 


Aus (3.19) folgt fiir 6=1, da die Summe der Maxima der Zeilensummen von 
M und K kleiner als 1 ist, daB sie zu M und P gehérigen geometrischen Reihen 
konvergieren. Wir k6nnen also (3.13) zur Fehlerabschatzung heranziehen und 


erhalten 
(eee oa = pee 
| Xe) ie 1 (2)| = 0,005 86 }’ 


wahrend der genaue Fehler 
0,004.98 
0,001 56 


betragt. Durch die Ungleichung wird ein Quader beschrieben, der ganz in F 
liegt. Damit sind alle Voraussetzungen des Satzes 2 iiberpriift. Das Verfahren 
(5.40) konvergiert also gegen die einzige Lésung des Gleichungssystems in F, 
und es gilt die genannte Abschatzung des Fehlers von *,. 


Stellt man das Verfahren (5.10) unmittelbar als gewohnliches Iterations- 


verfahren %,= WEY , dar, so wird ip wegen der Nichtlinearitat von R recht 
kompliziert, und die Ermittlung des dann benotigten Lipschitz-Operators fiir 


dieses T macht erheblich mehr Miihe als die Bestimmung von M und kK. Die 
Umschreibung von (5.10) als gewodhnliches Iterationsverfahren laBt sich noch 
allgemein durchfiihren, da hierzu nur die Auflésung von quadratischen Glei- 
chungen notig ist. 

Auf das in diesem Beispiel beschriebene Vorgehen zur Fehlerabschatzung ist 
man in solchen Fallen angewiesen, in denen sich *,=Tj,x, nicht allgemein in 
die Form des gewoéhnlichen Iterationsverfahrens bringen laBt. In der Regel kann 
man dann auch %, als Fixpunkt von J), nicht mehr exakt angeben, doch kennt 
man mit x,,_, meistens eine gute Naherung fiir x,,, die noch durch eine geeignete 
Iteration ,,im kleinen“‘ auf die gewiinschte Genauigkeit gebracht werden kann. 
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Schranken fir die Lésungen von Randwertaufgaben 
mit elliptischer Differentialgleichung 


A. G. MEYER 


Vorgelegt von L. COLLATZ 


Viele Probleme in Physik und Technik fiihren auf Randwertaufgaben mit 
Differentialgleichungen vom elliptischen Typus. Oft ist man auf Naherungs- 
funktionen angewiesen, die entweder die Differentialgleichung oder die Rand- 
bedingungen oder beides nicht exakt erfiillen, sondern nur bis auf angebbare 
Fehler, die sog. Defekte. Um trotzdem exakte Aussagen machen zu kénnen, 
méchte man Schranken fiir die Abweichung der Naherung von der Lésung 
kennen. 

Ausgehend vom Randmaximumssatz sind Fehlerabschatzungen fiir die erste 
und die dritte Randwertaufgabe bei gewissen elliptischen Differentialgleichungen 
von L. Corratz [1] und H. J. Grtnscu [6] durchgefiihrt worden. AuBerdem hat 
L. Cotiatz [2], [3] gezeigt, wie man bei Aufgaben monotoner Art Schranken 
erhalten kann. DaB auch in komplizierteren Fallen, wie z.B. bei nichtlinearer 
Differentialgleichung oder bei nichtlinearen Randbedingungen, die Aufstellung 
solcher Schranken oft médglich ist, will diese Arbeit zeigen. 


In Nr. 2 wird fiir eine gréBere Klasse von Randwertaufgaben gezeigt, daB 
sie Aufgaben von monotoner Art bilden. Das bedeutet, daB zwischen zwei 
Funktionen eine Ungleichung besteht, wenn diese Ungleichung fiir ihre Defekte 
beziiglich Differentialgleichung und Randbedingungen erfiillt ist. Schranken fiir 
die Lésung der Randwertaufgabe erhalt man mittels Funktionen, deren Defekte 
einheitliches Vorzeichen haben. 


Die Konstruktion von Funktionen, mit deren Hilfe eine Fehlerabschatzung 
nach Nr. 2 méglich ist, wird in Nr. 3 behandelt. Die Konstanten empfehlen 
sich dabei als besonders einfach zu handhabende Abschatzungsfunktionen. Ande- 
rerseits lassen Funktionen, die dem speziellen Charakter der Aufgabe angepaBt 
sind, scharfere Schranken erwarten. 


Nr. 4 schlieBlich bringt einige Beispiele. 


1. Bezeichnungen 


Alle betrachteten Funktionen sind reellwertige Funktionen von reellen Ver- 
anderlichen. Die Verdnderlichen sind die Koordinaten von Punkten im n- 
dimensionalen euklidischen Raum R,, mit rechtwinkigem Koordinatensystem. 
Als Abstand der Punkte voneinander wird der euklidische genommen. 
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Eine nichtleere Punktmenge im R, wird ein Gebiet genannt, wenn sie 
offen, beschrankt und zusammenhangend ist. Der Zusammenhang soll end- 
lich sein. Die abgeschlossene Hiille eines Gebietes wird abgeschlossenes Gebtet 
genannt. 


Als Rand I" einer Punktmenge M im R, wird die Menge der Punkte und 
Haufungspunkte von M bezeichnet, die nicht innere Punkte von M sind. Eine 
Punktmenge heiBt Hyperfliche F im R,,, wenn es m Funktionen 9;(&, -.., €,—1) 
(i=1,..., ”) gibt, so daB die Gleichungen x,=9,(&, -.., &,-1) eine umkehrbar 
eindeutige und stetige Abbildung der Punkte (&, ..., €,1) eimes einfach zusam- 
menhangenden abgeschlossenen Gebietes D des R,,_, auf die Punkte (4%, ..., %,) 
von F im R,, vermitteln. Die Hyperflache F im R,, wird glatt genannt, wenn 
auBerdem die Funktionen 9;(&,...,&,—1) (¢=1,..-.,) stetige partielle Ablei- 
tungen erster Ordnung nach allen Argumenten in D haben und von den ” Unter- 
determinanten (m —1)-ten Grades der Funktionalmatrix (@¢,/0&;) an jeder Stelle 
von D mindestens eine von Null verschieden ist. Die Hyperflache F im R, 
wird stetig gekriimmt genannt, wenn F glatt ist und die Funktionen @;(&, ..., €,_3) 
(t=1,..., n) stetige partielle Ableitungen zweiter Ordnung nach allen Argument- 
paaren in allen Punkten von D haben. 


Die Hyperkugel (oder kurz Kugel) mit dem Radius R um den Mittelpunkt 
A=a(q,,..-7,@,) im Ke ist die Menge der Ponkte P=(4,, ...,.4%,.) occ ae 


n 


>» (%;,—4a,)?=R?. Gelegentlich wird statt ,,Hyperkugel‘‘ auch ,,Oberflache der 


i=1 
Hyperkugel™ gesagt. Mit dem Inneren der Hyperkugel wird die Punktmenge 
n 
mit >) (%;—a,)?< R® bezeichnet. 
i=l 
Die partiellen Ableitungen einer Funktion werden meist kurz durch ange- 
hangte Indizes angegeben: 


Mit 7, (oder auch @//du) ist die Ableitung der Funktion f(x,,..., x,) in der 
Richtung einer Halbgeraden mit den Richtungskosinus 4;= cos (u, x,) (i=14, n) 


gemeint. Auch der Einheitsvektor mit den Komponenten yu; wird mit ju bezeichnet, 
so da man schreiben kann: 


fy = (1, grad f) =D gee 


Dabei kann yw selbst wieder ein n-tupel von Ortsfunktionen [gj (Maa oe 5) Hp) 
. . 2 
(7=1,...,) sein, wenn nur fiir alle betrachteten Punkte Dip, 1st, 

j=l 


2. Monotoniesatze 


2.1. Differentialgleichung und Randbedingungen linear 


Den Betrachtungen wird ein Gebiet G des Maifeet, 


x,-Raumes zugrundegelegt. 
Der Rand von G heiBe I. 


PL 


eS = 


a 
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In G sei ein linearer elliptischer Differentialausdruck vorgelegt : 


n n 
(2.1) L[u]=— >) a,,;%,,— >i b,u,+cu* mit den Eigenschaften: 
j=1 


(2.2) a,,, 6; und c sind gegebene beschrankte Funktionen in G+J°. Die Matrix 
der a,; ist symmetrisch in G+J’ und es gibt eine positive Konstante a, 
so da8 iberall in G+I" >i a,,&,¢,2a>)& ist fir beliebige reelle &é, 

k 


; 7 
,=1 (0 


U=1,...,%). Esast c20 im Gl. 
Auf dem Rande J’ sei ein linearer Differentialausdruck 
(2.3) R[w]=A,u—A,u, definiert mit den Eigenschaften: 


(2.4) A, und A, sind beschraénkte Funktionen auf J’, es ist A,=0, A,=0 und 
A,+A,.>0 auf I. 


Von den Funktionen, auf die die Differentialausdriicke L und R angewandt 
werden, wird vorausgesetzt, daB sie in G+" stetig sind, in G stetige partielle 
Ableitungen bis zur zweiten Ordnung einschlieBlich haben und an denjenigen 
Stellen von J’, an denen A,>0 ist, noch stetige partielle Ableitungen erster 
Ordnung nach allen Argumenten haben. 

Der Monotoniesatz laBt sich sehr einfach formulieren, wenn man voraussetzt, 
daB der Rand J’ eine stetig gekriimmte Hyperflache ist und «=» der ins Innere 
von G weisende Normalvektor auf J. In den Anwendungen liegen aber meist 
Gebiete mit ,,Kanten“ und ,,Ecken‘‘ vor und als Vektor w wird haufig nicht 
der Normalvektor vorgegeben, sondern z.B. der sog. Konormalenvektor. Um 
auch diese Falle und iiberhaupt méglichst viele Anwendungen zu erfassen, miissen 
kompliziertere Voraussetzungen tiber J’ und tiber ~ gemacht werden **. 

Zunachst wird der Rand J’ in Teilmengen aufgeteilt je nach Art der vor- 
gegebenen Randbedingung: 

I sei die Menge der Punkte von J’, in denen A,=0, A,>0 ist, 

I sei die Menge der Punkte von J’, in denen A,>0, A,=0 ist, 

I? sei die Menge der Punkte von J’, in denen A,>0, A,>0 ist. 

Dabei kénnen auch eine oder zwei der Mengen J” leer sein. Im letzteren Fall 
erhalt man die sog. erste (bzw. zweite, bzw. dritte) Randwertaufgabe, wenn 
Jet (bzwile = 1, baw. Raat?) ist. 

Voraussetzungen fiir den Rand I’: 

(2.5) I‘ bestehe aus endlich vielen Hyperflachen Ij’ (j=14, ..., mj). 

Die Hyperflachen J}? (j=1,..., m,) seien stetig gekriimmt, 

Die Hyperflachen 17? (j=1, ..., mg) seien glatt. 
Auf jeder der Hyperflachen J}? und J; ist also der innere Normalvektor y ein- 
deutig und stetig. 


* Die Indizes 7, k, ] dienen bei den Funktionen a@;;, b; und ¢ zur Unterscheidung, 
wahrend sie bei der Funktion u (und entsprechend bei allen Funktionen, auf die der 
Ausdruck L angewandt wird) die betreffende partielle Ableitung angeben. 

x* Fiir die Ubertragung auf nichtlineare Randwertaufgaben mu8B auBerdem auf 
die Voraussetzung der Stetigkeit von mw auf J’ verzichtet werden. 


20* 
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Voraussetzungen fiir den Vektor pu: 

(2.6) Der Vektor yu sei eindeutig auf jeder der Hyperflachen f? und 173 und es 
gebe eine Konstante d> 0, so daB (v, ~) 2d>0 dort ist mit dem zugehorigen 
inneren Normalvektor ». 

SchlieBlich miissen noch besondere Voraussetzungen fiir »y und mw in den- 
jenigen Punkten gemacht werden, die gleichzeitig zu mehreren der J? und J} 
gehéren (in den Kanten) : 

(2.7) Sei P ein Punkt von I’, der nicht auf J liegt, aber auf s verschiedenen 
Flachen J}? und J} mit s21. Seien rq), ..-, %) und fq), ---, 5) die damit 
in P vorgegebenen Vektoren. Dann soll gelten: 

a) Es gibt in P Linearkombinationen 


s Ss 
N= »> A; Vi) und M= 2, Bi My 
4=1 ps 


mit positiven Koeffizienten «;, 8; (#=1,...,s) derart, daB (N, %))>0 
und (M,%,)>0 ist fir 1=1,...,s. O.B.d.A. sei |N| =|M]=1 ange= 
nommen. 


b) Wenn P nicht zu J gehért (also nur auf s verschiedenen der Flachen J}? 
liegt), dann gibt es auBerdem eine Zah] e> 0 derart, daB (N(P), M(Q)) =e 
ist fiir alle Punkte Q von J” in einer ganzen Umgebung von P. 


N und M sind also nach (2.7)a Vektoren, die ins Innere von G weisen. Sie 
iibernehmen in den Kanten die Bedeutung, die » und w in den glatten Rand- 
punkten haben. Fiir s=1 sind die Voraussetzungen (2.7) von selbst erfiillt mit 
N=v und M=w wegen (2.6). In dem haufig auftretenden Fall w=y kann 
M=N genommen werden und damit (2.7) b entfallen. Die Voraussetzung (2.7) a 
laBt z.B. ohne weiteres zu, da8 sich in einem Randpunkt mehrere Hyperflachen 
auch unter spitzen Winkeln schneiden, schlieBt aber aus, daB zwei der Rand- 
flachen einander in einem Punkt von J’, der nicht auf J™ liegt, tangieren, d.h., 
da8 in einem solchen Punkt zwei innere Normalvektoren existieren, die zueinander 
entgegengesetzte Richtung haben. 


Monotoniesatz 1. Unter den Voraussetzungen (2.2) und (2.4) bis (2.7) ist die 
Randwertaufgabe Llu|=g in G, R[u|=y auf I’, von monotoner Art, d.h. wenn 
fiir zwei Funktionen u und u®, die sich nicht nur um eine additive Konstante 
unterscheiden, gilt 

Liu] < L[u®] in G, R[u®) < R[u®)] auf I, 
dann ist 
u®) < u® tiberall inG +T. 
M.a.W.: Wenn fiir eine Funktion w gilt 
Llw}S0 mG, R[w)S0 auf I, 


w nicht durchweg konstant in G+I, 
dann ist 
WSO therallinG+T’. 
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Im Beweis des Monotoniesatzes wird eine Eigenschaft benutzt, die eine ein- 
_fache Folgerung aus den angegebenen Voraussetzungen darstellt. Damit der 
_ Beweis nicht zu uniibersichtlich wird, sei zundchst diese Folgerung in einem 
) besonderen Hilfssatz bewiesen. 


| Hilfssatz. Se: P* ein Punkt des Durchschnitts D, der s stetig gekriimmten 
| Hyperflichen I)? (o=1,...,s) von I? mit s=1. Sei weiter K die Kugel mit dem 
| Radius R um den Punkt R= yee N* des von P* in Richtung N* aus- 


gehenden Halbstrahls, wobei N* = by Xe ViG) die nach (2.7)a in P* definierte Linear- 


_ kombination der inneren Ee vie) der I;? ist mit |N*| =1 und (N*, v6) >0 
fiir c=1,..., 8. Dann gibt es eine Zahl R, >0, so dap fir 0<RSR, de Punt. 
menge oie ganz auBerhalb K liegt ausgenommen der Punkt P*. 


Beweis..-1, Sel s—1." Dann 1st 
N*= v4). Der Einfachheit halber sei vor- 
genera der Koordinatenanfangs- 
, punkt in den Punkt P* und die Rich- 
_ tung der positiven x,-Achse in die 
| Richtung N* gelegt. Mit @ sei der Ab- 


CSO EE ee 


4 n—1 
» stand o=|/ >’ x? von der x,-Achse be- 
fat 


zeichnet. Die Ebene %,,=0 ist Tangen- 
_ tialhyperebene an J;? in P*. J;? hat in Fig. 1 
einer Umgebung U von P* die Dar- 
stellung x, =/(%,,...,%,—1), wobei f stetige partielle Ableitungen bis zur zweiten 
Ordnung einschlieBlich in U hat mit f=/,,=0 fiir k=1,...,~—1 in P*, wenn 
mit dem Index x, die Ableitung nach x, angegeben wird. Sei Z der Zylinder 
oSR,, |x,|<R,, wobei R, so klein gewanlt sei, daB Z noch ganz in U liegt. 


Nach dem Taylorschen Satz ist in Z: 


f(%,---, % n— Cigeesy + fal dee) 1X; +45 fonl® Hy, XY 
=0 at Se 


mit P= (9x,,...,0%, 1), @ noch von der betrachteten Stelle abhangig, aber 
0<%<1. Wegen Stetigkeit der zweiten Ableitungen von f/ in Z gibt es eine 
n—1. AuBerdem 


AES: 
ist |%,%;| 07 fir k,/—=1,...,%—1. Man erhalt also 


(Oye 0) ei: 


|f(%4, trey | ye Weel) ‘ o 
ee ies 1 
Dann hat Ry Min ( Mire = 
mit 7(P) der Abstand des Punktes P vom Mittelpunkt R=(0,...,0, R) der 
Kugel K mit dem Radius R, 0< RR, bezeichnet wird, dann erhalt man fiir 
die Punkte von J}?: 


P=o+(R—x,)?= 


also fiir 2-0, PEZ mit 2|x,|<F-(n—1)? ‘ound 0<RS— cee [7r(P)P>R, 
deh. 7(P)SR. 


1) s) die gewiinschten Eigenschaften. Denn wenn . 
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Die Punkte P von J}? haben also von J/ einen Abstand gréBer als R, aus- 
genommen P*. J}? liegt ~alst auBerhalb der Kugel K um P, mit Radius R, aus- 
genommen nur pe 

II. Sei s>1. Nach I gibt es dann zu jeder der Hyperflachen J;; 2 eine Zahl Ro, 
so daB die Hyperfliche J}? ganz auBerhalb der Kugel K, mit Radius Ri.) um 
den Punkt Py = P¥+4 Rey hy liegt fiir 0< Ros Ko, ausgenommen der Punkt P*. 
Seien die positiven Zahlen n,=(N*, »f)) eingeftihrt und 


(x) k= Min [sgl | 


Fig. 2 


R, hat bereits die gewiinschten Eigenschaften, d.h. der Durchschnitt D, der 
Hyperflachen Jj”, ..., J}? liegt auBerhalb der Hyperkugel K mit Radius R um 
den Punkt R= P*+RN* fir 0<RSR,, ausgenommen der Punkt P*. Das 
sieht man folgendermaBen ein: 

R sei fest gewahlt mit O0< RS Rj. 

K sei die Kugel mit Radius R um R= P*+R- N*. 

K, sei die Kugel mit Radius R,.)= si um Py = P*+ R,) v6). 

Wegen («) ist dann 0<Ry<R3. 
Mit 7(P) werde der Ortsvektor des Punktes P beziiglich P, bezeichnet, mit 


0,(P) sein Abstand von Ry, also 0,(P)= Z Vig) —RN*—r(P)]. 


Fallunterscheidung. 1.7(P)>R fir alle PED,, P+ P*. Dann ist der Satz 
richtig, D, liegt auBerhalb K, ausgenommen P*. 


2. Es gibt einen Punkt R mit RED,, R+P*, 7(B)<R. Dann kann & 
in keiner der Hyperkugeln K, liegen, denn Tj? liegt ganz auBerhalb K, wegen 
Rw So, also auch D, ganz auBerhalb ene ausgenommen P*. Es ist also 
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Bs) aR ydiir oa, dexs) duchy 


= (0) —RN*—7, 3 a) —RN*—1)>Ry= 5, 
nN, 
also , 
R2 
oe eta a eee 


o 


2 
(vg) — 2, N*, r) — ae = 
o 


oder 


R 
22 (ney Me N*, 7) <7? — R® 


Wegen |7| <R ist also ° 


a) (% —n,N*, 7r)<0 fiir o=1,...,s, und damit auch 
n—1 

D)- > &, (v%) —n,N*,r)<0, wo die «, die Koeffizienten der Linearkombination 
o=1 


n 
N*= >) a, 6) mit a,>0 sind. 

C= 

ae 


Mit S04 v5) = N* —a, 7%) und n,= (rt), N*) ist 


o=1 


il 
Dias (N* — a, v6), N*)=1— 4,0, 
o= 
und damit 
n—1 


Di % (¥%5) — Me N*, 7) = (N* — a, 1) — N* +a, n, N*, 7) = — a, (vt) — 0, N*, 7). 
Gail 


In dieser Gleichung ist die linke Seite negativ nach b), wahrend die rechte 
positiv ist nach a). Fall 2 ist also nicht méglich, D, liegt ganz auBerhalb K, 
ausgenommen P*, 


Damit ist der Hilfssatz bewiesen. 


Beweis von Monotoniesatz 1. Der Beweis wird indirekt gefiihrt und kniipft an 
die letzte Formulierung an. Die erate Formulierung erhalt man mit w= wu — uy, 

Es sei w>0 an einer Stelle in G+J”. Dann ist auch W>0 fiir den Gr6Bt- 
wert W von win G+I. Sei etwa w(P*) = W>0 in einem Punkt P* aus G+I, 
d.h. P* aus G+ 214-7? 4+ £3. 

I. P*€G. Dann ist w>0 in einer ganzen Umgebung U* von P* und damit 


2 ayn mint D2, w,=cw=0 in U* und w(P)<w(P*) in U*. Nach einem Satz 


on E. Hor 17] (man beachte insbesondere auch die FuBnote auf S. 149) ist 
dann w=W in U*. 

Es ist dann sogar w= W in ganz G+T’, denn sei U* die gré8tmégliche Um- 
gebung von P* in G mit w>0 in U*. Es ist also w=W in U* und entweder 
U*=G oder es gibt einen Punkt P**€G auf dem Rande von U* mit w(P**)=0 
Letzteres kann aber nicht sein wegen w= W in U*. Also ist U*=G und wegen 
Stetigkeit von w sogar w= Win G+I. 

II. P*€I. Dieser Fall ist nicht méglich, denn aus w(P*)>0 und 4,>0 
folgt R[w]—A,w>0 in P* im Widerspruch zur Voraussetzung R[w]<0. 
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Ill. P*cI3. Auch dieser Fall ist nicht moglich, weil dann w in einer von 
P* ins Innere von G weisenden Richtung anwachsen wiirde. Das sieht man so: 
P* kann nicht auch zu J! gehéren wegen II, wohl aber kann P* mehreren ver- 
schiedenen der Flachen J;* und J;? angehoren, die wir kurz Tyo I. nennen. 
Es seien vq), ---»% ) Mit s21 die entsprechenden in P* definierten Normal- 
vektoren und fq), ---, (ys) die dort vorgegebenen Vektoren. Dann gilt in P*: 
R[w] =A, ;-@ —Ao,5° yg SO aut I; fir +=1,..., 5 mit A; ;>0, A, ,;20, aber 
eins der A, ; positiv, da P* aus 1” ist, etwa A, ;,>0. Es ist also 


a ; B20 tir @es ters 
Oni = SAge. — B; -w mit B, +0 
und wegen w(P*)= W>0 er- 
halt man w,,,=2 0 (1=1,..., 5), 
Wyq,) 0. Nun gibt es nach (2.7) 


positive Zahlen B; (¢=1,...,5), 
so daB M= > B; uw ins Innere 
coal 


von G weist. Es ist 
s 
Yu = 26; a) = Wui,) > O- 
+= 


Das bedeutet, daB w anwachst 
in der von P* ins Innere von G 
weisenden Richtung von Mim 
Widerspruch zu der Vorausset- 
zung, daB win P* seinen GroBt- 
wert beziiglich G+ J’annimmt. 

IV. P*c€I. Auch dieser Fall kann nicht eintreten, weil auch er einen Wider- 
spruch zum Randmaximumssatz enthalt. Um das zu zeigen, wird zunachst ein 
Gebiet T konstruiert, auf dessen Rand P* liegt. Sodann wird eine Hilfsfunktion z 
ahnlich wie beim Beweis des Randmaximumssatz gebildet, die ihren GréBt- 
wert Z mit Z>0 im Inneren von T annimmt, obwohl dort L[z]<0 ist. 

P* kann nicht auch zu I! oder J’ gehéren wegen II und III. Sei P* ein 
Punkt des Durchschnitts D, der s Randflachen i (o=4,...;,5) DMt s 14 Danes 
sei s die kleinste solche Zahl, d.h. es gebe keinen Punkt P**€ J, der schon in 
einem Durchschnitt von weniger als s der Flachen J;? (j=1, ..., mg) liegt und die 
Bedingung w(P**) = W erfiillt. 

Nach Voraussetzung (2.7)a gibt es in P* einen Einheitsvektor N*, der ins 
Innere von G weist. Um den Punkt R= P*+RN* als Mittelpunkt sei K die 
Hyperkugel mit dem Radius R. P* liegt also auf der Oberflache von K. K* sei 
die Hyperkugel um P* mit dem Radius R*. R und R* werden noch genauer fest- 
gelegt werden. Das Teilgebiet von G, das im Inneren von K und K* hegt, 
wird T genannt (s. Skizze). Sein Rand 7; setzt sich zusammen aus einer Menge F 
von Punkten, die zu J’ gehéren und aus Teilen K; und K* der Hyperkugeln K 
und Kk*. 


Voritbergehend wird der Ortsvektor beziiglich des Punktes Presale ibe) 


mit 7(P) bezeichnet und seine Koordinaten mit § (jays, 90) ds bee inp 
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(4%, --+,%,) sel 7(P)=(&,...,&)=(%1— 41, ..., %,—28). AuBerdem sei p= 
|r| =/> & eingefithrt. 

Nun werden R und R* festgelegt. 

R>0 wird so klein gewahlt, daB gilt: 


1. F und die Punkte auf der Verbindungsgeraden zwischen P* und P, liegen 
im Inneren von G, 


2. K hat mit D, keine Punkte gemeinsam ausgenommen P*, 


3. F enthalt nur Punkte der Hyperflachen Jj? (o=1,...,s) von I. 

R*>0 wird so klein gewahlt, daB gilt: 

4.w>Oin 7, 

Be (7(Q), M(Q))<o auf F, wo M(Q) die nach (2.7)a im Punkte QCF vor- 
handene Linearkombination der dort definierten i.) (o=1, ..., ) ist. 


Da eine solche Wahl von R und R* méglich ist, muB kurz gezeigt werden. 
Fir genitigend kleines R ist erfiillt 1., da N* ins Innere von G weist, 2. nach 
oben angegebenem Hilfssatz, 3., da Hyperflachen als abgeschlossene Mengen 
definiert wurden. Sei also R klein genug fest gewahlt. Wegen Stetigkeit von 
win G+JI' und w(P*)=W>0 gibt es eine Zahl R,>0, so daB w>0 ist fiir 
0<R*<R, in allen Punkten von G, die im Innern von K* liegen. SchlieBlich 
gibt es nach Voraussetzung (2.7)b eine Zahl e>0, so daB (N*, M(Q)) =e ist 
fiir alle Punkte Q aus J”? in einer ganzen Umgebung von P*. Es gibt also eine 
Zahl R¥>0, so daB (N*, M(Q))=e>0 aut F ist fiir 0<R*< RF. Wenn man 
vortibergehend mit £*(Q) den Vektor von P* nach Q bezeichnet (s. obige Skizze), 
dann ist 7(Q) = —RN*-+ p*(Q) und fiir QEF gilt |p*(Q)| < R*. Sei R; eine Zahl 
mit 0<R;<Min(R#, eR), dann ist (7(Q), M(Q)) <0 auf F fir 0<R*<Rk,, denn 


(7(Q), M(Q)) =— R(N*, M(Q)) + (6*(Q), M(Q)) <— eR+eR=0. 
2e>0 <|p*|<Rs<eR 
Nun sei auch R* fest gewahlt mit 0< R*<Min(k,, R;). Dann erfiillt also R* 
die Forderungen 4. und 5. 

Der so konstruierte Bereich T hat die gewiinschten Eigenschaften. T ist 
nicht leer, denn die Punkte P= P*+aN* mit 0<a<R* liegen in T. Nun 
wird die Hilfsfunktion z=w— W-+ 6h gebildet mit h—=e~*® —e-**®, Die Kon- 
stanten 6 und « werden noch festgelegt. 

Zunachst ist 
L{w|S0 in T, da T Teilbereich von G ist, 

U— Wain 

w=W auf K;, denn K,c(G+T), 

w<W—x auf K* mit geeignet gewahltem ~>0, denn die Punkte von K7 

liegen alle in G oder in einem Durchschnitt von weniger als s der 
Flachen Jj? (o=1,...,8). Dort ist tiberall w<W, es gibt also 


Jo 
wegen Abgeschlossenheit von A} ein x, so daB wlW—x auf 


Toast. 
—wy <0 auf F, denn auf I;* ist R[w]= —A, ,w,,,.0 und A, a> 0 fut o = 
1,..., 8, also —w,,,<0 und damit —wy<0 wegen M= DBote mit 


B,>0 nach (2.7)a. 
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Fiir die Hilfsfunktion h liefert der Differentialausdruck (2.1): 
ee Lh] = — 02 4 Yay E.G + 20 Dd (a, +h,&) Hel —e FP). 
k,l=1 j=l 
Es ist > = 9?> (R—R*)?>0 in T+I7, also nach Voraussetzung (2.2) 
j=l 
4> Oy €,62 44> GS 4a(R—R*?=A>S0 mit A=4a(R— k*)?. 
k, J=1 j=l 
Ebenfalls nach (2.2) sind die Koeffizientenfunktionen von L beschraénkt in 
G+TI und damit auch in T+J;,, es gibt also Konstanten B>0 und C>0, so 
daB 2¥( (a;;,+0;&)<B und c[1i—e*@-<C ist in T+17. Damit kann 


man eae eT ih]<—a?A+aB-4+C. Seiaso groB gewahlt, daB e**L[h] <0 

ist in T+. Wegen e*’>0 ist dann auch L[h]}<0O in T. SchlieBlich sei 

2 ra ~ gewahlt. Dann ist 6>0 (wegen 0<R—R*<R) und 

eu (R— —e-% ‘ 

dh<xaut K*. Auf F ist — hy <0 wegen —hy=— (gradh, M)=2ae-%® (r,M) <0 
nach 5.). Zusammenfassend erhalt man fiir die Funktion 6/: 


L[éh]<0 in T, 64 =0 in P* und auf Ky, 6hSx auf KZ, —dhy<0 auf F. 
Fiir die Funktion z=w—W-+ 6A gilt also: 
Eii<0 “indy 


rey shi) 
2=0 autA, ond i. 
—%Z%y<0 auf FP. 


Sei Z der GroéBtwert von zin T+J/7. Er werde angenommen in einem Punkt 
P**, Es kann nicht P** = P* sein, denn in P* ist z,,>0 und M weist ins Innere 
von T. Wegen z(P*)=0 ist also insbesondere Z>0. P** kann also auch nicht 
auf Ky oder Af liegen, weil dort z<O ist. SchlieBlich kann P** auch nicht 
auf F liegen, denn in den Punkten Q€F ist zy>0 und M(Q) weist ins Innere 
von G, also auch ins Innere von T, wenn nicht gleichzeitig QC Kp oder QC K* 
ist. Auf K; und K} war z<0<Z schon festgestellt worden. P** muB also im 
Inneren von T liegen. Es ist z3=const wegen zy>0. Dann ist Z>0 aber ein 
Widerspruch zum Randmaximumssatz, wie schon im Teil I dieses Beweises fiir 
die Funktion w gezeigt wurde, denn z erfiillt in T+J} die Voraussetzungen, 
die win G+ TI erfiillt. 

Einen Punkt P* mit den verlangten Eigenschaften gibt es also nicht in 
G+I’, die eingangs des Beweises gemachte Annahme w(P*)>0 ist falsch, es 
ist w<0 tiberall in G+J" in Ubereinstimmung mit der Behauptung des Satzes. 


2.2. Differentialgleichung oder Randbedingung nichtlinear 
Es soll die Monotonie der Randwertaufgabe nachgewiesen werden fiir den 
Fall, daB Differentialgleichung oder Randbedingung nichtlineare Differential- 
ausdriicke enthalten. Der Beweis wird durch Reduktion auf den linearen Fall 
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erfolgen, d.h. es wird gezeigt werden, daB8 fiir Funktionen, die den Voraus- 
setzungen des hier formulierten Monotoniesatzes 2 geniigen, auch geeignete 
lineare Differentialausdriicke existieren, so daB diese Funktionen schon dem 
Monotoniesatz 1 geniigen. 

Den Betrachtungen wird wieder ein Gebiet G des x1, ..., x,,-Raumes zugrunde- 
gelegt. Der Rand von G heiBe wieder I’. 

An die Stelle des linearen Ausdruckes (2.1) tritt hier der allgemeinere Diffe- 
rentialausdruck zweiter Ordnung 


(2.8) Ele) = = Ol 445.20, Lg y by Uy vos Uns Mars Mra «+99 Une) 
kurz geschrieben Bia = — shaw) U, U;, Up) - 


Wahrend in (2.2) die Voraussetzungen fiir L in G+J" formuliert werden konnten, 
muB fiir die Formulierung der Voraussetzungen fiir E ein Bereich H des (n? +2" +1)- 
dimensionalen Raumes mit den Koordinaten %;, 4, u;,u,; benutzt werden. H 
habe die Eigenschaften: 


(2.9) a) H ist beztiglich der Koordinaten u, u;,u,; konvex, d.h. mit zwei 
Punkten (€;,«, %;, %,;) und (&;, B, B;, B,1), deren erste m Koordinaten 
iibereinstimmen, liegt auch ihre ganze Verbindungsstrecke in H, also 
alle Punkte 


(&,a+A4(B — a), a; + A(B; — 0t;), Og + A(Baz — &p21)) mit OSA<1. 
b) Der abgeschlossene Bereich ae ist in der Projektion von H in den 
x,-Raum enthalten, d.h., wenn der Punkt (%;, 4, 4;, u,;) alle Lagen 
in H durchlauft, dann soll (x,) alle Lagen in G+J" annehmen. 
Der Differentialausdruck (2.8) erfiille die Voraussetzungen 


(2.10) D(x;, u, u;, 1) hat stetige partielle Ableitungen erster Ordnung nach 
den m?-+n-+41 Argumenten 4, u;, 4; 1n a Es gibt eine Zahl a>0, so daB 


fiir beliebige reelle £ (j=1,... n) SP ansS &2a- Se iiberall in H 
i] 
: a@ 
gilt, wo ®,, = Ze ist. Es ist ®, = 8 <5 in H, 
Fiir (x,,..., %,) auf I” sei ein Differentialausdruck erster Ordnung vorgegeben: 


(2:44)-S| | SV See, Fn My «2+ ) Un)» 
kurz geschrieben S[u] = WY (x;, u, u;). 


In H erfiille dieser Differentialausdruck die Voraussetzungen 


(2.12) Fiir (x;) auf hat Y(x;, u, u,;) stetige partielle Boe. erster Ordnung 
nach den »+1 Argumenten u und u; (j=1,..., 2), es ist Yi (%,, u, 4) 20 
dort und 


WY, (x;,%, %) + I/ Dd Wi? (x,;, u, uj) >O (Wurzel nichtnegativ genommen). 
R=1 


An jeder Stelle (x,;) auf J” 14Bt sich in H definieren 


(2.13) 4;=¥,(x;,4,,), Ao=|/ Dd Pa (*;, 4, u;) und fiir A,+-0 der Einheits- 
ae Wi, (4%, U, uj) 


we (p=4\ geo 
A, > J d 


vektor w=(ty, +--+) My) Mit “y= — 
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A,, A, und ys hangen also an jeder Stelle (x;) von J’ noch von den u, u; ab. 

Bei beliebiger Wahl dieses (n+ 1)-tupels in H gelte 

(2.14) [laBt sich mit Hilfe von (2.13) in Teilmengen J", J und J aufteilen 
(je nachdem A,=0 oder A,=0 oder A,, A,>0), die den Voraussetzungen 
(2.5) gentigen. Der in (2.13) definierte Vektor w und der innere Normal- 
vektor v auf J’? und J" erfiillen die Voraussetzungen (2.6) und (2.7). 


Von den Funktionen, auf die die Differentialausdriicke E und S angewandt 
werden, wird vorausgesetzt, daB sie in G+TJ’ stetig sind, in G stetige partielle 
Ableitungen bis zur zweiten Ordnung einschlieBlich haben und auf J? und J? 
noch stetige partielle Ableitungen erster Ordnung haben. AuB8erdem sollen sie 
in H liegen, d.h. fiir (x;) aus G+" soll (x,;, u, u;, %,_;) em Punkt in # sein. 

Monotoniesatz 2. Unter den Voraussetzungen (2.10), (2.12) und (2.14) 1st, 
wenn ein Bereich H des x;, u, u;, u,,-Raumes nut den Ergenschajften (2.9) zugrunde 
gelegt wird, die Randwertaufgabe E[u|=g in G, S[u|=y auf I’, von monotoner 
Art in H, d.h., wenn fiir zwei Funktionen u® und u®, die beide in H liegen und 
sich nicht nur durch eine additive Konstante unterscheiden, gilt 


E[u] < E[u®)] in G, 
S[u®] <= S[uw®] auf I, 
dann ist UD < yf2 tiberall in G+TI. 

Beweis. Es geniigt hier, lineare Differentialausdriicke L und R anzugeben 
derart, daB L die Voraussetzungen (2.2) und R die Voraussetzungen (2.4) bis 
(2.7) erfillt und L[w®]<L[u] in G und R[u®]< R[u®)] auf I gilt. Dann sind 
alle Voraussetzungen des Monotoniesatz 1 erfiillt, der unmittelbar die obige 
Behauptung ausspricht. 

uw) und u® sind zwei den Voraussetzungen des Satzes geniigende, sonst be- 


hebige aber fiir den Beweis feste Funktionen. Sei wieder w= u® — u® eingefiihrt. 


Dann ist 
Elu) — El) = — @(x;, u® + w, u + w,, u@ + w,)) + 


5 D(x;, ur), ui), ul?) =< 0. 


Der Mittelwertsatz liefert (wegen H konvex und (2.10)): Mit0<8 <1, 8=9(x,) 
ist in G 


E(u] — E[u®) =— w- @,(x;,u+ bw, u” + dw,,u)+9w,)) — 
— Sw, ©, (x, u + 910, ui + 6 10;, uf) + 8 ay)) — 


f= 
— Lin (%;,u? +9 w, uw +9w,,u2) +9 w,)). 


Nun werden fiir (x;) aus G die Funktionen u®)(x,), w(x;) und O(x,;) auf der 
vechten Seite eingesetzt und die Bezeichnungen eingefiihrt: 


eee gs u® (x,) + B(x) w(%),...,-.-), 
ANS suiaeeeee i) 

B,(%;) = fie , UP) (26) + B(%,) W(%), 0.5004), 
C(x) =— B, (;, wo) (x,) (x) w(%,)2..)2.). 
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Damit 148t sich in G ein linearer Differentialausdruck definieren: 


Liu] = = Di Ags Us: — 7s Bu, + Cu. 
$,¢t=1 vy=1 
Dieser Differentialausdruck erfiillt die Voraussetzungen (2.2), denn es gilt: 
A,,, B,, C sind beschrankte Funktionen in G+J’, da DOr, 2, nach (2.10) 


stetig, also beschraénkt sind in H und (x;, u+8uw, u)+8w,, u?)+ Ow, )) in H 
hegt fir (%,) aus G+I’. Die Matrix der A,, ist symmetrisch in G+TJ und es 


gibt eine Zahl a>0, so daB in G+I" s Ay Ge =a ye ist fiir beliebige reelle 


s,t=1 =F 
&,(r=1,...,), da ®, und @, , diese Eigenschaft haben nach (2.10). SchlieB- 
lich ist C=0 in G+J" wegen C=— ®, >0 in 7 nach (2.10). 
Wenn man den linearen Differentialausdruck L auf die Funktion w anwendet, 
erhalt man 


Liw] = — 03 (@,,, + ®,,,) @.:— 2, ®,, vw, —D,w und mit w,,=w,,: 
Gyigeil vy=1 

L{w] =— > ®,,, st ap SOE, w, — B®, w = E[u] — E[u®] <o. 
Syl y=1 


Eine entsprechende Betrachtung la8t sich fiir den Randausdruck S durchfiihren. 
Es ist auf J” 


S[u] = S[a] = P(x, 0 + w, ul + w,) — (x, u, uP) 


}= 


= oP, (x;,,u + dw, ul + 3w,) + 21 5 Fie (x;, vw) + Ow, ul + Pw,) 


mit O< <A, O=V (%,). 
Fiir (x;) auf I” werden wieder die Funktionen wu (x;), w(x;) und #(x,) ein- 
gesetzt und die Bezeichnungen nach (2.13) eingefiihrt: 


A, (x) = ¥%(«,;, uw + bw, u + dw,) , 


are or, wl) On, a bu) 


und insPunkten von /;;in denen A,=-0 istiw=(M, .-«;4,) mit 
aoe (4,, u9+0w, us + 3w,) 


= — TEA 5 FA acy We 


Damit laBt sich auf J’ ein linearer Randausdruck Fk definieren: 
R[u] = A,u— A,u,. 


Dieser Randausdruck hat die gewiinschten Eigenschaften, denn fiir A, und A, 


gelten die Voraussetzungen (2.4), da Y, und yz >» WF die entsprechenden Voraus- 


setzungen (2.12) erfiillen. AuBerdem geniigen r , e und jw den Voraussetzungen 
(2.5) bis (2.7), da (2.14) gilt und (x,;, u+0w, u!”)+0w,) ein Punkt in H ist 
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) auf J. Wenn man R auf die Funktion w anwendet, erhalt man mit 
n a (%;, uw + Ow, us?) ee dw) 

pte = — Ds A, (#1) ; 

k=l 


fiir (4; 


Me 


w, = (grad w, u) = 
k 


R[w] =v A,— w, A,= w+, (x,, u + Ov, ...) se, Wi (x;,u-+ bw, ...) 
= 
= Stu) — S[u®] <0. 
Es ist also R[w®]<R[u®] auf L. 


Damit sind alle Voraussetzungen des Monotoniesatz 1 erfiillt, dieser liefert 
die Behauptung wu) <u® in G+I’. 


1 


2.38. Evfassung nichtbeschrankter Gebvete 


Nun sei B eine offene und zusammenhangende Punktmenge (endlicher 
Zusammenhang, aber nicht notwendig beschrankt). Mit J’ werde wieder der 
Rand von B bezeichnet, d.h. die Menge der Punkte und Haufungspunkte von B, 
die endlichen Abstand vom Ursprung haben und nicht innere Punkte von B sind. 

In Erweiterung der Bezeichnung von Nr. 1 wird hier eine Punktmenge auch 
dann Hyperflache genannt, wenn jede einfach zusammenhangende abgeschlossene 
Teilmenge eine Hyperflache im Sinne von Nr. 1 bildet. 

Fiir die Differentialgleichung wird die allgemeinere (auch nichtlineare) Gestalt 
(2.8) zugelassen, die Randbedingung wird als linear vorausgesetzt: In einem ge- 
eigneten Bereich H des x;,u, u;,%,;-Raumes mit den Eigenschaften (2.9) beziig- 
lich B+ I sei wieder ein elliptischer Differentialausdruck E[u]=— O(x;, u, u;, Up)) 
definiert wie in (2.8), der wieder die Voraussetzungen (2.10) erfiillt. Auf dem 
Rande J’ sei ein linearer Differentialausdruck R[w]=A,u—Ag,u, definiert wie 
in (2.3), der die Voraussetzungen (2.4) erfiillt. 1” geniige wieder den Bedingungen 
(2.5), die Vektoren y und w den Bedingungen (2.6) und (2.7). 

Von den Funktionen, auf die die Ausdriicke F und R angewandt werden, 
wird vorausgesetzt, daB sie in B+ J’ stetig sind, in B stetige partielle Ableitungen 
erster und zweiter Ordnung haben und in den Punkten von J’ mit A,+0 noch 
stetige partielle Ableitungen erster Ordnung. Wenn B nicht beschrankt ist, muB 
bei der Formulierung der Randwertaufgabe zusatzlich angegeben werden, wie 
sich die Lésung uw in den uneigentlichen Randpunkten von B verhalten soll. 

Sei o(P) der Abstand des Punktes P vom Ursprung, K, die Menge der 
Punkte von B+J’ mit @(P)=r. Das Verhalten von wu in den uneigentlichen 
Randpunkten von B kann dann festgelegt werden durch die Bedingung 
(2.15) lim {sup | w(x;) — y,(%)|}=0, 


1—> CO 


wo y,(%;) eine geeignet vorgegebene Funktion auf K, ist (nicht notwendig be- 
schrankt fiir 7—> oo), 


Monotoniesatz 3. Wenn unter den genannten Voraussetzungen fiir zwei Funk- 
tionen u und u, die beide in H liegen und sich nicht nur wm eine additive Kon- 
stante unterscheiden, gilt E[w®] < E[u®] in B, 

R[w®)] < R[u®) auf I; 
lim {sup (w) — u®))\<0, 


1 co” Ky 
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dann ist yO) = 2) in B eras 
Beweis. Sei wieder w=u — y®) gesetzt. Die Behauptung lautet w<0 in 
B+I. Sei w(B)>0 in einem Punkt P, aus B+I. Wegen lim {sup w} <0 gibt 


es ein %, so daB w(P)<w(A) ist fiir alle Punkte P aus B+JI' mit 0(P)=7. 
Wenn es also einen Punkt R, mit w(R)>0 gibt, dann wird das Maximum von 
w(P) in B+J" in einem Punkt P* von B+I' angenommen mit o(P*)<7,. P* 
kann aber nicht in B liegen, denn nach Beweis von Monotoniesatz 2 gibt es einen 
linearen elliptischen Differentialausdruck L mit L[w]<0 in B und in Teil I von 
Monotoniesatz 1 wurde nachgewiesen, daB dann P* nicht in einem endlichen 
Punkt von B liegen kann. Nach den Teilen II, III und IV von Monotoniesatz 4 
kann P* auch nicht auf J’ liegen. Es gibt also keinen Punkt P, mit w(R)>0, 
es ist w(P)SO0in B+. 


3. Zur praktischen Durchfiihrung von Fehlerabschatzungen 
Die Randwertaufgabe laute: In einem geeigneten Gebiet G mit dem Rand I’ 
wird eine Funktion u(x%;) gesucht mit 
(3.1) E[ujJ=g inG, 
R[uj=y aut I, 


wo g bzw. y beschrankte Funktionen sind in G bzw. auf I’. 


Die Differentialausdriicke E[u] = — O(x,, u, u;, u,;) und Riu] =A,u—A,u, 
sollen in einem Bereich H des %;, u, u;, u,;-Raumes die Eigenschaften haben, 
die fiir die Anwendung des Monotoniesatz 2 gefordert wurden*. Die Existenz 
einer Lésung der Randwertaufgabe (3.1) wird vorausgesetzt. 


Nach Monotoniesatz 2 erhalt man: 


(3.2) Wenn fiir zwei Funktionen uw und u®) aus H, die sich nicht nur um eine 
additive Konstante von u unterscheiden, gilt 


Elu)]<g<E[u®] inG, 


Riu] Sy SRW] auf L, 
dann ist 
uD <u < y in-G sd: 


In der Praxis wird man haufig zunachst eine Naherung v fiir die Lésung u 
der Randwertaufgabe bestimmen und dann nach Schranken fiir den Fehler u —v 
fragen. Um (3.2) anwenden zu kénnen, miissen Funktionen gY und g® gesucht 
werden, so daB u=v-+ q die Bedingungen von (3.2) erfiillen. Diese sind dann 
Differentialbedingungen fiir die g@. Um dies auch in der Darstellung zum Aus- 
druck zu bringen, sei der Differentialausdruck 


(3.3) E*(q] = Elv+q] — El] 


* Die Voraussetzungen reduzieren sich natiirlich auf die von Monotoniesatz 1, 
wenn £ linear ist. 
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eingefithrt, der die gleichen Eigenschaften wie E hat, solange v und v+q in H 
liegen, was immer vorausgesetzt wird. Im linearen Fall stimmen £ und E£* 
iiberein, Beim Aufsuchen der g“ empfiehlt sich die Benutzung der Form 


(3.4) E*[q] = — > ®y, (%, 9 +99, 4 +BY, Ye PGs i Gees 


s,t=1 
— 1 @, (4,0 + 9g, 0 + OG; Mp1 + POE) I — 
v=1 


— &,(%;,,0 +99, 0; +99;, U1 + OG 1) -9 


mit 0<#<1, D=0(X%;). 
Zur Naherung v seien positive Konstante ¢;, 6; ((=1, 2) bekannt mit 


(3.5) —éeSEluj]—Elv]) Se, inG, 
— 6, = Riu) —R[v] S6, auf L. 


Dann erhalt man aus (3.2) die Fehlerabschatzung: 


(3.6) Wenn fiir zwei Funktionen g® und g®), fiir die v+g® in AH liegen und 
u—v—q® nicht konstant sind in G+’, gilt 


Be Fin tes Pe | Aa 


Rig’) at = oy" Moyes ig Cant ae, 
dann ist 
qQ<su—vsq -in G+TI. 


Hier sei eingeschoben, daB mit Hilfe der Fehlerabschatzung wieder die Ein- 
deutigkeit der Lésung der Randwertaufgabe (3.1) in H nachgewiesen werden 
kann. Denn wenn es zwei Lésungen u und u* in H gibt, die sich nicht nur um 
eine additive Konstante unterscheiden, dann ist mit w* als Naherung fiir u: 
&= €2=6,=0,=0 und damit q@s<u—u*<g® nach (3.6), wenn E*[g?]|<0< 
E*{g®] in G und Rig] <0<R[q] auf I ist. Dies aber ist schon erfiillt fiir 
gQU=q®=0 in G+I. 


Nach (3.6) erscheint es wichtig, eine Funktion qg(x;) zu kennen mit 


(3.7) EX) ce: en. G. 
R(qj2oé auf L 


bei vorgegebenen positiven Konstanten ¢ und 6. 


Man kann (bei der zweiten Randwertaufgabe) leicht Falle angeben, in denen 
es eine solche Funktion nicht gibt. Allerdings lassen sich zu vorgegebenen E 
und R, falls in der Randbedingung A,>d>0 ist, die Konstanten a, b, w, a; 
(@=1,..., 2) in der Funktion 
(3.8) q=a(b—e*") mit g => (4,—a,)? 

j=l 

so wahlen, daB (3.7) erfiillt ist, vorausgesetzt wie immer, da8 v-+g noch in H 
liegt. Die Fehlerabschatzung mit dieser Funktion ist meist weder einfach noch 
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besonders scharf. Man wird deshalb entweder einfachere Funktionen (z.B. Kon- 
stanten) fiir die Fehlerabschatzung suchen oder Funktionen, die der speziell 
vorliegenden Aufgabe besser angepaBt sind. 

Durch konstante g kann z.B. abgeschatzt werden, wenn der Differential- 
ausdruck FE iiber die Voraussetzung (2.10) hinaus die Bedingung ®,< —c in H 
mit konstantem c>0 und der Raudausdruck R iiber die Voraussetzung (2.4) 
hinaus die Bedingung A,=d mit konstantem d>0 erfiillt. Dann ist, wie man 
aus (3.6) mit (3.4) ersieht: 


(3.9) Min (— ag 3) <u vSMax(, “) inca 


Wenn die Naherung v schon die Randbedingung R[v|=y auf J’ erfiillt, kann 
auf die Bedingung 4,=d>0 auf I’ verzichtet werden und man erhalt 
(3.10) a ee inG+Tr. 


Entsprechend kann auf die Bedingung ®,<—c<0 in H verzichtet werden, 
wenn v schon die Differentialgleichung E[v|=g in G erfiillt. Man erhalt dann 

) 
(3.14) 


OF 
1S4-—vs3 
fe = 


d 


in G+TI". 


In allen drei Fallen (3.9) bis (3.11) wird natiirlich vorausgesetzt, daB wu —v nicht 
konstant ist in G+J" und auBer v auch noch v+q™ in H liegen. 

Bei linearer Differentialgleichung E[u]=L[u|—g in G kann man geeignete 
Abschatzungsfunktionen g® angeben, wenn Funktionen z und z* bekannt sind mit 


(3.12) Eee an, Riz 0. aut £, 
EN AE ine a a al eres EN eg Oe 

Dann ist 

(3.13) — 62—6,2%* Su—vS &2-+ 62% in G+, 


wenn nicht #—v—q™ konstant sind in G+I’. 


Fiir die erste und dritte Randwertaufgabe, d.h. A,;=2d mit konstantem d>0, 
kommt man schon mit der Funktion 2 allein aus: 


(3.14) —az—" Su—vSagz+9 in G+, 


wo a =Max(0, a— ct 2), a,=Max(0, e,—o* ist mit c¥= inf[e(x,)], 


c der Koeffizient von w im Differentialausdruck L[w] nach (2.1). Es ist 
wieder u —v—g nicht konstant in G+TI’ vorausgesetzt. 


Fiir L[w] = — Au eignet sich im Falle der ersten oder dritten Randwertaufgabe 
z.B. die Funktion 
1 n 
(3.15) L=>, ei 4) 
J= 


mit L[z]=1. Es muB nur 0 so groB gewahlt werden, daB R[z|20 ausfallt. 
Die a; wird man so wahlen, da8 z méglichst klein ist. 
Arch, Rational Mech. Anal., Vol. 6 21 
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Fiir L[uj]=—Au-+cwu und beliebige Randbedingungen eignen sich in der 
Ebene z. B. die Funktionen 


(3.16) z=1(1—af,iVeo)) mit e=Vie— 8)? + (y — 7) 
(Getz) a ; [1 — a, cosh Ye (x — &) — az cosh Ye (y — 9)]. 


Es ist wieder L[z]=1. R[z]=0 laBt sich erreichen durch geeignete Wahl der 
a, dy, 43, &,. Im Fall der zweiten Randwertaufgabe muB hier noch eine Funk- 
tion z* mit den Eigenschaften (3.12) bekannt sein. Fiir den Kreis eignet sich z. B.: 


(3.18) z*=afy(tyeo), 
fiir das Rechteck z.B.: 
(3.19) z* = a, cosh Ve (x — €) + a, cosh Yc (y — 9) 


mit passend gewahlten a, a, a2, &, 7. 


Zum AbschluB sei noch ein Hinweis fiir die Behandlung nichtlinearer Diffe- 
rentialgleichungen gestattet. Gelegentlich fiihrt die Aufspaltung des Differential- 
ausdruckes in einen linearen und einen nichtlinearen Anteil leichter zur Auf- 
findung geeigneter Abschatzungsfunktionen : 


(3.20) Elu| = Liu] + N[u]. 


Wenn z.B. speziell N[w| nicht von den Ableitungen von uw abhangt, d.h. N[w| = 
— Y(x;, uv) ist, dann kann wegen der Voraussetzung (2.10) die Aufspaltung immer 
so vorgenommen werden, daB L die Bedingungen (2.2) erfiillt und ¥%<0 in H 
ist. Dann gilt 


E*[q] = L(g) + N[v +9] — N[o] = Lg] — 9 (x; v +89), 


0 


lA 


also 
= Lig) ming SO: 


Ee* 
ia] = Lig|— tur g= 0, 


so dal} sich die fiir lineare Differentialgleichungen gegebenen Formeln (3.13), 
(3.14) und gegebenenfalls auch die Funktionen (3.15) bis (3.19) verwenden lassen. 


4. Beispiele 


Die Existenz einer Lésung wird bei allen behandelten Randwertaufgaben 
vorausgesetzt. 


4.1, Minimalfliche (Plateausches Problem) 
Flu] =— (1+ 4) 4, +24, Uy Uy — (1+ 02) Uyy =O > ING: |e ye 


sgn x « 2 fir |a[(—=41, lyls4 
R{u| =u4= 3 


s 


seny : fur apace slg) ed 


» 
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Der Differentialausdruck E erfiillt die Voraussetzungen (2.10) fiir beliebige 
Argumente x, y, u, U,, Uy, Uz~, Uzy, Uyy- Die Voraussetzungen (2.4) bis (2.7) fiir 
den Randausdruck R sind trivial erfiillt. Fiir den Bereich H brauchen deshalb 
nur die ersten beiden Koordinaten beschrankt zu werden: |x| <1, |y| <4. 


Aus der Eindeutigkeit der Lésung ergeben sich die Symmetrien 
u(x, y) = u(y, x) =—u(x,—y) unddamit w(x,0)=«(0,y) =0. 
Es gentigt also, Schranken aufzustellen fiir die Aufgabe: 


Hiav=0inG*? 0<x2<15 0<y=1 


= free 41,9 4S 
R[u] =u= = fie yA) he 
OF NiO ,e (Ve 12 under y— ty [a = 1. 


Fiir die Naherung v seien Funktionen verwandt, die schon die Randbedingungen 
erfiillen: 


p= SE (At (1 — 9) (1 = y) [iy + (22+ 9?) + ay 2? 9? +» Sh. 
Bei Verwendung nur eines Parameters a, erhalt man fiir « die Schranken 


SU — 0,25 (1 — #4) 1-4] Sus. 
Bei Verwendung von mehreren Parametern wird man zunachst den Defekt 
|E{v]| méglichst klein machen und sodann ay variieren, bis E[v] einheitliches 
Vorzeichen in G* hat. So erhalt man hier die Abschatzung 


E = {4+ (4 — x?) (1 — y?) [— 0,085 + 0,004.55 (x? + y2) — 0,109 75 x? y}} 
< 0,065 ~” (1 — 2) (1 — y?). 


Die gréBte Abweichung liegt an den Stellen x?= y?=4. Dort ist 


0,148 <u(/%, /4)<0,167 bei Verwendung von ay, 
0,1548 < w(//4, 4) S0,1646 bei Verwendung von ay, a, 42. 


4.2. Geschwindigheitspotential der kompressiblen Strémung um eine Kugel 


Nach Transformation auf dimensionslose GrdBen erhalt man fiir das Potential 
gy die Randwertautfgabe: 


E(g] = —Agv+4or{grad p- grad K + (%—1)KAg}=0 in B: r>1 
R[g]}=—9,=0 auf lj: r=1 
lim ( sup |y—rcos#|) =0. 


Te OOS 
Dabei ist K=| grad y|”. « und x sind Konstante. 


21* 
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Die Voraussetzungen (2.10) sind erfiillt, wenn man sich auf einen Bereich H 
<_!_.1 _¢ mit einem positiven e*. H geniigt den 
4 +H a . . . . 
Voraussetzungen (2.9). Die Voraussetzungen (2.4) bis (2.7) sind wieder trivial 
erfiillt. Die Bedingung tiber das Verhalten in den uneigentlichen Randpunkten 


ist angegeben. Monotoniesatz 3 kann angewandt werden. 


beschrankt mit |grad p?S 


Aus der Eindeutigkeit folgen wieder die Symmetrien y(7,0)= (7, —0)= 


2) 9 


: =0=5 zu betrachten, wenn man 


—(r,x—*#) und damit insbesondere or, —F)=0. Es geniigt also, 


das Problem in dem Bereich B*: 7>1, 


noch auf dem hinzukommenden Rand J} (0= =e = : 
Rig] =gv=0 stellt. 


Fiir die Konstanten wurde x =1,405 gesetzt und «=0,4**. 


r= 1) die Randbedingung 


E. Lamia [8] hat eine Naherung @* aufgestellt, die bereits die Randbedingung 
erfiillt: 
pt =RF+RA+# BF 
mit 
ig 


Pizacos 0; = — (5 cos3#+ 3 cos#), 


Pi oe (63 cos 5% + 35 cos 30 + 30cos%), 


R= Ya, (2 =1,3,5). 


Tabelle der a;,: 


=I + 1,00000000 — 2 
2, + 0,563 78018 —0,058600 53 + 0,003 047 62 
4 _ ++ 0,100 29679 — 0,013 963 64 
5 — 0,041 843 81 — 0,061 002 67 +0,018 54781 
6 os a= + 0,003 53812 
7 — 0,003 231 58 — 0,014 52218 —0,01994805 
8 -+0,016033 63 + 0,016 73200 + 0,009 494 86 
10 + 0,001 689 24 -++0,004070 94 + 0,001 44623 
dA — 0,004 464 41 — 0,005 30080 — 0,001 464 44 
14 ++ 0,000 587 64 -+0,00045900 + 0,00007968 
Fiir die Fehlerabschatzung sind hier die Funktionen q= 5 fiir 0<b<1 geeignet. 


Sei b= 4. 


1 
ice - hat die Stromungsgeschwindigkeit die lokale Schall- 


geschwindigkeit erreicht. 


* Fiir [grad o|?=. 


** Fiir «a=0,56 ist bei v=1, oe — die lokale Schallgeschwindigkeit bereits iiber- 
schritten. 
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se gD = — 02008 ata) 0,04 


Vv Vr 


und » = p+ q" erhalt man E[g] <0 < E[p®] in B*, 


— ot) = — 0,004 <0 <0,02= — g”) auf J, 
(4) — 0,908" 03s (2) 
Y V7 << 7 Pp auf 5, 
lim | sup_(p—reos#)|=o=lim{ sup (rcos# — o®)}. 
¥—> co —nr<0<n f 1 \_ nega 


Der Monotoniesatz 3 liefert damit die Aussage 


: QQ) <p <q in B*¥++4+]3, 
ltt 
8 . 
om = ya hive ag BY Bese ee 


Fiir den Punkt v=1, #=0 erhalt man also z.B. die Schranken 
1507 == OU, O) 455507. 


Eine entsprechende Abschatzung laBt sich fiir die Stromfunktion Y dieses 
Problems durchfiihren. Fiir die von E. Lamia [8] angegebene dreigliedrige Nahe- 
rung ¥* erhalt man: 


— 0,013 77 Prsin? dS YW — W* < 0,00427 |/r sin? 3. 


4.3. Randwertaufgabe mit nichtlinearer Randbedingung 
SAw=4 , in.G:, |x). Ay, 


lz] =4, |y| S14, 


gradu| =u auf I: 
| nl =1, |e] St. 


Die Voraussetzungen (2.10) sind ohne Einschrankungen erfiillt. Der Rand- 
ausdruck lautet hier: 


S[u] =u —|grad u| =u — Vu 


Die Voraussetzungen (2.12) sind erfiillt wegen 


Woa40 ound 1 Vee +2 = 20. 


Die nach (2.13) definierten GroBen sind 4;=4,=1, w= seed up Det Rand 
co) 


I’ gehért also ganz zur Menge J® (A,, A,>0) und erfiillt die Voraussetzungen 
(2.5) (stiickweise glatt). Der Vektor ~ und der innere Normalvektor » geniigen 
den Voraussetzungen (2.6) und (2.7), wenn grad w iiberall auf J’ ins Innere von 
G weist. Dies ist eine Einschrankung fiir den Bereich H. 
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Sei H die Punktmenge des 4%, Y, 4, Uy, Uy, Ux, Uny» u, ,-Raumes mit 
Jen pS oa fie | el ee 
—seny-u,2a fir |y|=41, |x| S41, 


wo a>0 eine vorgegebene Zahl ist, die beliebig klein gewahlt werden kann. 
Fiir H gilt (2.9). Damit sind alle Voraussetzungen von Monotoniesatz 2 erfillt. 
Wenn man, was hier leicht méglich ist, eine Naherung v benutzt, die schon 
die Differentialgleichung erfiillt, dann liefert Monotoniesatz 2 die Abschatzung: 
Aus |S[u]—S[v]| <A auf IP mit konstantem A folgt |w—v| <A in G+J’, wie 
man mit =v —A, qg®=v- A leicht sieht. 
Die Eindeutigkeit der Lésung in H liefert die Symmetrien: u(x, y)=u(y, x)= 


wet v2 
e + Ay + 4,01 + Ago 


verwendet mit v,—Re[(x+7y)*”]. Mit a)=0,848044, a,=—0,019562, ag= 
0,000529 erhalt man fiir den Randdefekt |S[w] — S[v]| $0,005 auf /’ und damit 


u(—x,¥). Fir eine Naherung sei deshalb der Ansatz v = 


|«—v| $0,008 inG+T. 

Insbesondere wird also z. B. 
0,8430 < u(0, 0) <0,8531, 
0,4297 S u(1, 1) < 0,4398. 


Diese Arbeit ist ein Teilauszug aus der von der Mathematisch-Naturwissenschaft- 
lichen Fakultat der Universitat Hamburg im Januar 1960 genehmigten Dissertation 
des Verfassers. Herrn Prof. Dr. L. Corratz danke ich fiir die Anregung und sein 
forderndes Interesse. 
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1. Introduction 


This paper considers supersonic spatial flow fields which are conical in the 
sense originally introduced into aerodynamics by BUSEMANN [J]. In such a 
flow the velocity and the conditions defining the state of the gas, e.g., the pressure 
and temperature, are constant on rays through one point of the physical space, 
called the center of the conical field. Generalized conical flows, being flows 
in which these quantities are homogeneous of degree higher than zero (e.g. [2]), 
and conical fields used to construct flows of incompressible fluids (e.g. [3]) are 
thus excluded. 

The treatment of conical flows within the frame of the linearized theory 
was initiated by BusEMANN [4] and has been given much attention since by 
many authors (e.g. [5] and [6]). Also, higher order approximations were con- 
sidered, where either linear theory ([7], [8]) or the non-linear solution for the 
axially symmetric flow around a circular cone ([9], [10]) were chosen as a starting 


point. 
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For a long time the development of non-linear theory has been restricted 
to the study of particular examples, such as the solution for the axi-symmetric 
flow around a circular cone, given by BusEMANN [J], [//] and by TayLor & 
Macco tt [1/2], and other types of conical flow [73]. In the non-linear theory the 
flow around a specific body is obtained as a numerical solution of the differential 
equations. Methods of construction of such a numerical solution have been 
discussed by MASLEN [14], FoweE zt [15] and FERRI, VAGLIO-LAURIN & NEss 
[16], [17]. However, in solutions found by these methods, for example for the 
flow around a flat delta wing with supersonic leading edges as given by MASLEN 
[14] and Fowe t [15], certain discrepancies arise. It is of interest therefore 
to consider in more detail the non-linear equations governing conical flow. This 
has been done by BULAKH in a number of papers [/8]—[22], partly commenting 
on papers cited above [2/], [22]. In the present paper the properties of non- 
linear isentropic conical flow are studied through a different approach, con- 
sidering surface elements of integral surfaces of the non-linear equation from 
the point of view of differential geometry. For this purpose the hodograph 
transformation of isentropic conical flow, as studied first by BUSEMANN [11], [23] 
and later by Giese [24], NrKOL’sKiI [25] and RyzHov [26], appears to be par- 
ticularly useful. 


2. Irrotational conical flow analyzed on the unit sphere 
around the center of the flow field 


In the physical space let a right-handed co-ordinate system x, y, z be fixed 
with the origin at the center of the conical field, and let u,v and w be the 
components of the velocity along the axes, respectively. The coefficients of 
viscosity and heat conduction of the gas are assumed to be zero. If it is assumed 
moreover that the flow is isentropic, the three conservation laws (mass, momentum 
and energy) yield the following equation: 


uP v2 w 
u, (1 =| o,(4 75] | w,(4 =] 
UU 


ad a (ut, Uy) 


(1) 


where a is the local velocity of sound, related to the velocity components by 


ia he tere Bet ‘ 1 =f 
Pm FE a — 2 ot wt) = FEN gt eh oa, (2) 


ax is the critical velocity of sound, and y is the ratio of specific heats (y= 2 ; 
Cy 


If the flow is free of rotation, a velocity potential may be defined in the 
usual way such that 


Pe Uc, Py = 4, a (3) 
Equation (1) then becomes 


u? ( v2 2 
raft‘) +m,{1-2) +0) 
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This equation may now be specialized for conical flow, using the property 
that the velocity does not change along rays through the center of the field. 
The velocity and the state of the gas depend therefore on two length co-ordinates 
instead of three, as in the general case of spatial flows. These co-ordinates may 
be taken arbitrarily to be x and y, and the flow may be considered in a plane 
z= const. It is convenient also to analyze the flow on a unit sphere with center 
at the center of the conical field. The plane z=1 is then a plane tangent to this 
sphere at the point (0,0,1). By rotating the x, y, z axes, any point on the 
sphere can be taken as this point. The conical properties of the flow may be 
expressed by putting 

p=zF(E,n), (5) 
where = x/z and y= x/z, and from (3) then follows 
“w= Sb eS = § 7 F: (6) 


With the aid of (5) the following equation may be written for equation (4): 


F.\1+e@— oe ! 2F.,|En eee eney) | 


(7) 


+ F,, [t+ 9-2-2") =o. 


ar 


In order to determine under what conditions this quasi-linear homogeneous 
partial differential equation is elliptic, parabolic or hyperbolic the character- 
istic directions may be determined from the equation 


pee Set" —aGp_ fev enstpcens 


dE | char. @? dé a 


E [tpg — Coa") —o. (8) 


a2 


For the local investigation of the flow it is convenient to use a co-ordinate 
system in which velocity components are measured along and perpendicular 
to the radius under consideration. The x, y and z axes are therefore rotated in 
such a way that the z axis has the direction of the radius under consideration 
and the x axis is in the direction of the velocity component perpendicular to 
the radius. The rotated system may be indicated by X, Y and Z, and the 
velocity components by U, V and W, respectively, U being the velocity com- 
ponent perpendicular to the radius and W the velocity component along the 
radius. If the flow is analyzed on the unit sphere, U is the velocity component 
tangent to the sphere. A streamline on the sphere may be defined as the inter- 
section with a streamsurface, which may be constructed as a cone with the 
vertex at the center of the conical field and going through a spatial streamline. 
A streamline on the unit sphere is thus directed along the X axis or the velocity 
component U and is named a conical streamline. The characteristics given by 
(8) may be drawn on the unit sphere and are called the conical characteristics. 
For the point (0, 0, 1) in the X, Y, Z system the conical characteristic directions 


read dz i 
dH char. 


Lg a (9) 


a2 
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The conical characteristics thus subtend the Mach angle, defined in terms of the 
velocity on the surface of the unit sphere, with a conical streamline. Let us 
call this Mach angle the conical Mach angle 7, and the Mach number defined 
in terms of the velocity U the conical Mach number M,(=U/a); then equa- 
tion (9) can be written as 

a= ete 

(FFF non <= stan ge. (10) 
In analogy with two-dimensional plane flow the velocity normal to the conical 
characteristics is equal to the speed of sound. The conical characteristic direc- 
tions are real and have two different values for M,>1; the equation is then of 
hyperbolic type, and the flow will be called conical-supersonic flow. For M,=4 
the two conical characteristic directions are coincident, real and perpendicular 
to the conical streamline; the equation is parabolic, and the flow will be termed 
conical-sonic. If M,<1, the conical characteristic directions are imaginary; the 
equation is of elliptic type, and the flow will be called conical-subsonic flow. 
Points on the unit sphere where U=0 will be called conical stagnation points*. 
It is of interest to consider the relation between the characteristic surfaces in 
the spatial flow and the conical characteristics thus defined. A disturbance 
generated at a point of the flow field travels along the characteristic surface 
through that point; the characteristic surface starts as a characteristic or Mach 
cone. All characteristic surfaces emanating from points on the ray of the 
point considered may be constructed by a similarity transformation of the given 
characteristic surface with respect to the center of the conical field. The envelope 
of all characteristic surfaces so obtained is thus a conical surface, which intersects 
the unit sphere along the conical characteristics going through the point of inter- 
section of the ray under consideration. This can be seen in the following way. 
Since the velocity component normal to the envelope, being a surface to which all 
characteristic surfaces are tangent, is equal to the velocity component normal 
to a characteristic surface, this velocity is sonic. Since the normal to the conical 
envelope is perpendicular to every curve on the envelope through the point 
considered, the normal at a point of intersection with the unit sphere also is 
perpendicular to the ray through that point and the intersection of the 
envelope with the unit sphere. The normal to the envelope thus is tangent 
to the unit sphere and perpendicular to the intersection of the envelope with 
the unit sphere; the velocity component along it is sonic, which shows that 
the intersections are conical characteristics. A conical disturbance may be defined 
as a disturbance generated with equal strength all along one ray. The conical 
characteristics are then the lines along which conical disturbances travel. An 
analogy with two-dimensional plane flow may be interpreted in the following 
sense. When a two-dimensional sound source is moving in a plane, the sound 
signals emitted by the source are propagated by sound waves which, if the velocity 
of the source is supersonic, form an envelope, being characteristics or Mach lines. 
In conical flow the disturbances travelling along the characteristic surfaces may 
be thought of as propagating on the unit sphere along curves which are the 


* A disadvantage of this term is that the source-like character of a conical stagna- 
tion point (U =0, thus W +0, because g +0) is not expressed by it. 
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intersections of the characteristic surfaces with the unit sphere. When the flow 
is conical-supersonic, the curves of intersection belonging to the characteristic 
surfaces of one ray have an envelope which consists of conical characteristics on 
the sphere. 


It may be noted that the characteristic surface emanating from the center 
of the conical field is of special interest. This surface coincides with the Mach 
cone with its apex at the center of the field and intersects the unit sphere along 
a conical-sonic line. This may be seen from the fact that a characteristic cone 
at some point of this conical characteristic surface is tangent to it along the 
radius through this point. The velocity component normal to the radius is 
therefore equal to the velocity component normal to the characteristic surface 
or Mach cone. Since the latter velocity component is sonic, the velocity com- 
ponent normal to the radius is sonic, and M,=1. 


The velocity at the center of the conical field is in general multivalued; 
hence the Mach cone is not necessarily circular. Also it may be noticed that 
the influence of the center of the conical field is not restricted to the down- 
stream interior of the characteristic cone from the center of the field. Actually 
it is confined to a cone, which may be constructed by connecting the center of 
the cone field with all points on the unit sphere by a curve formed by the 
conical characteristics which envelop the conical-sonic line. This may be seen 
in the spatial flow field by assuming that the disturbances originating in the 
center of the field travel initially over the conical characteristic surface from 
the center of the field. According to HuyGENs’ principle each point reached 
on this surface in turn acts as a source of disturbances which propagate along 
the characteristic surface of that point. The latter surface is tangent to the 
conical characteristic surface from the origin but does not necessarily lie inside 
the downstream interior of it. The envelope of all characteristic surfaces starting 
at the characteristic surface from the center of the field thus bounds the region 
of influence of that center. This envelope intersects the unit sphere along the 
aforementioned curve. 

In order to illustrate the quantities defined on the unit sphere, parallel flow 
throughout the physical space may be considered as an example of conical flow. 
The center of the field may be chosen arbitrarily at any point of the flow. The unit 
sphere is sketched in Fig. 1. The conical streamlines, being the intersections with 
the sphere of meridian planes through the diameter connecting the conical stagnation 
points, go from one conical stagnation point to the other. Conical-subsonic and 
conical-supersonic regions may be distinguished, separated by circular conical-sonic 
lines which are the intersections of the Mach cone from the center of the field with 
the sphere. The envelope of characteristic surfaces emanating from points of one 
radius consists of two planes which pass through this radius and are tangent to the 
Mach cone through the center of the field. The conical characteristics are the inter- 
sections of these planes with the sphere. 


The difficulty in trying to determine the physical properties of conical flow 
lies partly in the fact that regions of conical-subsonic and conical-supersonic 
flow may occur simultaneously in a flow field, in which case (7) is of the mixed 
type. Since partial differential equations of the mixed type are in general 
difficult to handle, we use an approach which lends itself well to this specific 
problem. 
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The properties of conical flow are given by the states of motion that a gas 
particle in a conical flow might have. , 

The motion of a gas particle may be determined by its velocity, acceleration 
and higher derivatives of the velocity up to an arbitrary order. It is seen from 
equation (6) that when for a point on the unit sphere the velocity is given by 
its components along the axes, the co-ordinates of the integral surface [(é, 7) 
representing a flow where such a situation occurs and the tangent plane at 
some point on the surface are determined. Similarly, the magnitude and direction 
of the derivative of the velocity up to some order determines the surface element 
of F at that point, up to that order. 


chorocteristic surfoce 


through the centre of the field chorocteristic surfoces emanoting 


7, from points on radius OA 


s conicol chorocferistics 
conical stagnation | ff 


point 


conical—supersonic flow Mp >| 
conical-sonic line Mp=1 
conical—subsonic flow Mp< / 
— conical stognation point U=t 


Fig. 1. Description on the unit sphere of parallel flow throughout physical space 


To all geometrically possible surface elements that the differential equation 
permits to be surface elements of an integral surface thus correspond possible motions 
of a gas particle in a conical flow. The physics of conical flow is thus reduced 
to the differential geometry of integral surfaces of the differential equation for 
conical flow. In combination with (7), it is useful from the point of view of 
differential geometry to consider an equivalent equation obtained by the 
Legendre or hodograph transformation, which will now be discussed. 


3. The hodograph transformation of irrotational conical flow 


The hodograph transformation of irrotational conical flow is obtained in the 
usual way by introducing the Legendre potential 


z(u,v) =uE+un—F(é,n), (11) 


hea tote (12) 
Comparing (6) and (11), we see that 


so that 


1 (u,v) = — w(u, v). (13) 
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The opposite of the velocity component along one axis may therefore serve as 
the Legendre potential, being a function of the other two velocity components, 
and the conical flow may be represented by a surface in the hodograph space. 


By use of (13) in equation (12) follows 


oD, = — &; vW,=— YY. (14) 
Differentiating (14) yields 


d&=— dw, = — (w,, du+w,,a0v), 


(15) 
dn = —dw,=— (w,,du+w,, dav). 


If the Jacobian determinant 4 =w,,,w,,—w%, is finite and different from zero, 
(145) may be solved for du and dv, and we have 


du= pe. (— Wyy GE ae Wy» 4N) , 


(16) 
Me Ti (@,, aE — W,, an). 

The transformation is then locally one-to-one; that is, to one point on the unit 
sphere in the physical space corresponds one point on the surface in the hodograph 
space, and vice versa. Singularities in the transformation occur for 4=0 or 
A->oo. Furthermore, since 


du=u,dé+u,dn =F, dé +E, dn, 


(17) 
dv =v, 45 + v,dy = f,,dé + Edn, 


comparing this equation with (16) yields the second derivatives of F, and the 
following differential equation may be written instead of equation (7): 
(u+w wy) pteee ai 


2 (u-+w wy)? 
v a2 


a2 


| 
Cis 


2 Wy v Be W, 
18 
(u+w wy)? bs) 


az 


== (i) 


2. 
Ee 2541 + w2 


It may be noted that if w(w, v) is a solution of this equation satisfying given 
boundary conditions, —w(—wu, —v) is also a solution. Therefore, as in all 
isentropic flows subject to boundary conditions in the form of prescribed stream- 
lines, the flow may be reversed. The significance of this fact will be discussed 
later in relation to the limit cones. 

Equation (18) is of the same type as (7), and the characteristics of the equation 
may be written as 

2 ) 

Pe a i 


char. 


(19) 


1 L we ee =(. 
T a? 
For the conical characteristic directions on the hodograph surface in the 
U, V and W co-ordinates we have 
dV errr 
cake Ss = Vices 20 
oe + I i : ( ) 
Elliptic and hyperbolic regions again may be distinguished on the hodograph 
surface, having two different imaginary and real conical characteristic directions, 
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respectively, and corresponding to conical-subsonic and conical-supersonic flow. 
If the flow is conical-sonic, the two conical characteristic directions coincide 
along the U axis. 

As in two-dimensional plane flow the velocity along the conical characteristics 
on the hodograph surface is sonic. The angle between the velocity vector q 
and a conical characteristic on the hodograph surface is therefore the angle «, 
defined from the local Mach number M(M-=gq/a) and given by 


tana = +)/M?2—1. (21) 


It is equal to the angle between the velocity vector ¢g and the characteristics 
in the hodograph plane for two-dimensional plane flow, which are the well- 
known Prandtl-Meyer epicycloids. Thus, when a cone passing through a conical 
characteristic on the hodograph surface and having its vertex at the origin of 
the hodograph space is developed onto a plane, a Prandtl-Meyer epicycloid is 
obtained [24]. 


4. Differential-geometric description of the hodograph transformation 
for irrotational conical flow 


From (14) it follows that the radius in the physical space is perpendicular 
to the surface element at the corresponding point on the hodograph surface. 
The sphere obtained by collecting at one point the unit vectors along the normals 
to the hodograph surface is the unit sphere in the physical space, as discussed 
before. From differential geometry the transformation from the hodograph space 
to the physical space may be recognized as being the spherical or Gaussian trans- 
formation of the hodograph surface. 

An analysis of the geometry of the hodograph surface may start by investi- 
gating the properties of its curvature. If further exploration is of interest, 
third and higher derivatives may also be taken into consideration. 

A study of the curvature of the hodograph surface can be based upon DUPIN’s 
indicatrix. Duprn’s indicatrix is formed by laying out along the normal sections 
of the surface distances equal to the square root of the absolute value of the 
radius of curvature of those sections. A curve related to Dupin’s indicatrix is 
obtained as the intersection with the hodograph surface of a plane parallel to 
the tangent plane at the point under consideration at a distance C such that 
higher order derivatives may be neglected with respect to the second order 
derivatives when the shape of the curve of intersection is determined. If the 
second derivatives are continuous, a Taylor expansion yields the following 


equation for this intersection in the U, V, W system, attached to a radius for 
which U=U, and W=W;: 
W=W,-C, 


2, 
Wou(U — U,)? + 2Woy(U — U,) V + Wy, V2 20 =0. a 


This is a conic*, being an ellipse when (Woy Wy — Wov),=Ko,>0, a hyper- 
bola when K,,<0 and a parabola degenerated into two parallel lines for Kg =0. 


es Dupin’s indicatrix is similar to this conic for elliptic and parabolic points. 
For hyperbolic points Dupin’s indicatrix consists in two conjugate hyperbolas, 
which are similar to this conic if C is given two equal and opposite values. 
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Points on the surface are called elliptic, hyperbolic and parabolic points, respec- 
tively. For an elliptic point the surface is curved in the same sense in all directions; 
for a parabolic point the same is true, while in one direction the curvature of the 
surface is zero; and for a hyperbolic point curvatures of the surface of opposite 
signs occur. For a hyperbolic point there are two directions for which the curva- 
ture of the surface becomes zero. 

The axes of the conic are in the principal directions, and the corresponding 
radii of curvature are called the principal radii of curvature 0, and g,. Through- 
out this paper @, will be chosen as the major principal radius of curvature and 
0, as the minor principal radius of curvature. The principal directions are given 
by the angle «,. with respect to the U axis (%,, measured positive in counter- 
clockwise direction), where «. may be deduced from 


2Wouy 
tan BOS ayaa (23) 
Lines on the surface which at each point are tangent to one of the principal 
directions are called lines of curvature. The radius of curvature of the curve 
of intersection of a plane through the normal to the surface and making an 


angle « with the major principal direction is given by EULER’s theorem, 


+ — | cos?% + 

Og 

or in the U, V, W system, ™ a 
1 


mais = Wyy cos?« + 2Wyy sina cosa+ Wry sin?a, (25) 


where « is the angle with respect to the U axis. 


1 


sin? %, (24) 


The normal curvature x, of a curve on the hodograph surface is the opposite 
of this value, whereas the other intrinsic second-order parameter of a curve 
on a surface, the geodesic torsion, may be expressed as 


Tz = Wyy cos? « — (Wyy — Wry) sina cosa — Wyy sin? a, (26) 


where T, is positive if the normal to the surface turns to the right when moving 
along the curve. 

From (23) and (26) it then follows that the geodesic torsions of the lines of 
curvature are equal to zero. If transformed to the unit sphere, the lines of 
curvature have the same direction there as on the hodograph surface. 

The ratio of corresponding line elements along the lines of curvature on the 
unit sphere and the hodograph surface may be obtained from (15); the result 


is the Rodrigues equations 
ire Goma. ds. (27) 


where 9, and og, are the radii of curvature in the principal directions 1 and 2, 
respectively, and the indices h and ph refer to the hodograph surface and to 
the unit sphere in the physical space, respectively. 

A direction making an angle «, with principal direction 1 on the hodograph 
surface (%, being measured positive in the counter-clockwise direction) makes 
the angle «,, with the direction corresponding to principal direction 1 on the 
unit sphere. From (27) these angles are related by 


tan &p, = ra tan a. (28) 
2 
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The area of a surface element dA,, on the unit sphere and the area of the 
corresponding surface element on the hodograph surface dA, are connected with 
the Gaussian curvature K, by the relation 


dA pp eee ey 
dAy, Carer dsm me : (29) 
Since 
Ke =Wou Ww — Wiv = A, (30) 


the Gaussian curvature or the Jacobian determinant is thus seen to be equal 
to the ratio of magnitudes of corresponding surface elements, being positive for 
an elliptic point, zero for a parabolic point and negative for a hyperbolic point. 
In relation to equation (28) it may be deduced that the image on the unit sphere 
of a closed curve on the hodograph surface is traversed in the same sense as 
the curve on this surface if K,>0, and in the opposite sense if Kg<0. Singu- 
larities in the transformation are to be expected for Kg =0 and Kg—o (A=0 
and A—oo). 

In addition to the Gaussian curvature, the mean curvature may be defined 
as the sum of the principal curvatures and given by 


S 1 1 
Ky =— +4+-— =— (Wou + Wy), (31) 
Q1 Q2 
where 0, and @, are chosen to be positive if the hodograph surface is convex 
towards the direction of the positive W axis. 
Other directions of interest are conjugate directions. The directions defined 
by the angles «, and a, are said to be conjugate if 
Qo 


’ 


-. 
tan o, tan %, = — 
Q1 


or 
Wou + Woy {tan (0% + a) + tan (a + a) }+ Wy tan (+ o,) tan (% + o},) =0. (32) 
By use of (28), equation (32) may be written as 

tan a, tan a, = — 1. (33) 


The image on the unit sphere in physical space of a direction on the hodograph 
surface is therefore perpendicular to its conjugate direction on the hodograph 
surface. 

Directions which are self-conjugate are asymptotic directions, which trans- 
form perpendicularly to their images. They may be obtained from (22) or (32): 


= Woy + V- (WouWyy — Wav) ; G4) 


(a 
dU asympt. : Wry 


Thus for an elliptic point (Kg>0) the asymptotic directions are imaginary ; 
for a hyperbolic point (Kg<0) there are two real asymptotic directions; and 
for a parabolic point (Kg=0) the two real asymptotic directions coincide along 
the axis of the parabola, which is in the principal direction where 0+. 

The differential-geometric properties of a surface, thus summarized, may be 
used to express the physical quantities of interest in terms of geometrical 
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properties of the hodograph surface, which is the particular surface under con- 
sideration. This surface is characterized by differential equation (18), which at 
a point of the surface, since Wy =W,=0, yields the relation 


Wou + Wry (1 — M?) =0, (35) 


first given by BuSEMANN [23]. The radii of curvature Ry and Ry in the directions 
of the U and V axis, respectively, then satisfy the equation 


Ry (1 — M?) + Ry =0. (36) 


From (35), (20) and (32) it may easily be seen that the conical characteristics 
on the hodograph surface form a conjugate net. The conical characteristics on 
the unit sphere of one family are therefore perpendicular to the conical character- 
istics on the hodograph surface of the other family. The + (—) sign in equation 
(10) corresponds to the + (—) sign in (20), and the same situation is encountered as 
in two-dimensional plane flow, as has been shown by GIESE [24] and RyzHov[26]*. 

Another system of conjugate directions is formed by the streamlines on the 
hodograph surface and the lines of constant speed (or a, f and M are constant). 
Partial differentiation of Pay te 37) 


with respect to v, for a point X = Y=0, yields 
Ode ie (38) 


The lines g = constant on the hodograph surface therefore intersect the U axis 
perpendicularly and are thus perpendicular to the conical streamline on the 
unit sphere. Furthermore, the lines g= constant on the unit sphere are normal 
to the acceleration, which has the same direction as the conical hodograph 
streamline. The streamlines and the lines g= constant therefore form a con- 
jugate system on the hodograph surface. Further it is seen that the lines g = 
const. bisect the angle between the conical characteristics on that surface [24]. 
The variations of g,a and M are found by partial differentiation with respect 
to wu of equations (2) and (37). For a point X=Y=0 


dq _ oq M, 


Sone big aaah ig (39) 
da Oe = pA M 
Be ae ) 


dM eM _ M, (1 


+ 1M} 4 
Gees Mig E® aie am) 


The variations of the velocity component normal to the radius U, and of 
M, are also of interest. For a point X=Y=0 it can be shown, after some 


calculation, that eU, 
aU =1+WWoy, (42a) 
Be 
=WWoy, (42b) 


* See also [28], p. 483, note 9. 
Arch. Rational Mech. Anal., Vol. 6 22 


310 J. W. REyN: 


and aM, 1 y—1 are 
ey eS ier th + WWouy 5 (43a) 
eM, 1 
= Wo: (43 b) 


from which the direction of a line M, = constant is found to be given by 


eat 
ian i M?+WWouy 


aU M,=const. W Wov 
The direction of the conical streamline on the hodograph surface, which also 
is the direction of the acceleration, may be seen from (32) and (38) to satisfy 


Dy EN ee Woy 
tan B Cae ee (45) 


(44) 


The magnitude of the acceleration along a streamline in the physical space 
may be expressed in terms of the curvatures of the hodograph surface in the 
following way. Along a streamline in the physical space 


(=) ee 
dz],  w” 
(2 set 
ie).= w’ 


where the index s refers to conditions along the streamline. By differentiation 
of (14) we have 


and (46) 


(dx), ee Oy (dz), er z(dw,),, 
47 
(dy), = oc? w, (dz), — 2(dw,),. 
or from equation (46) 
(d x), ued 2 [Wun Se (48) 
lease 
—2 [Wy (dU), + Wyy (20), 
(ay), = —# an Eero eae 49) 
1 AR 
Uv 
Again using (46) leads to* 
dv a= (y+ W Wy) Wy 4 — (U+W Wy) Wuy 
d (aa |= (w +w Wy) Qyy— (v+w wy) Wyy ’ et 
an 
a4x\.=— = Wuu Wyy— Cy an 
(dx), ae tree eee (du), , 
dy), = —zy——__Wuu Woo 
Os — 29) tpg O40) Dey Gy 
dz). = — > Wuu Wyv— Who 
oe on (U-+ W Wy) Wyy—(U+W BW) Wyy Natl 
For the « component then follows 
(du\ __ 1 (U+wwy) ey y»—(v+w2,) wy, 
% | ds )= q2 Pe Aes SAA ie ee (52) 


* Equation (50) may be used as a check on (45). 
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and for the v component 


(Z)=-2 (ut+ww,) Wy y4—(U+W Wy) Wy» 5: 
ds qé Wy y Yyy— Wav ( 3) 
Furthermore, 
‘dw du dv 
(ae) = (ae htt Ca). Po 
so that from equations (52)—(54) follows 
(= oa. 1 (w+ww,) (Wy Wyo Wy Wu») — (U+ ®,) (Wy Wu v— Wy Wan) 55 
ds I q 2 Wy y Wyy— ny Gs ( 
By use of the relation 
d d d ) 
(ae) (as) + (as). ° (Ge), pe) 


and of equations (52), (53) and (55) the acceleration along the streamline is then 
found to be 


dq @(54) = 1 (UW Wy)? Wy y—2 (UL Wy) (UW Wy) Wyvt(ULW Wy)? yy (57) 
dt ds G2 Wy y Wyy— Wie : 
or, by (18) 
dq = a” (1 +435) Wy y— 2Wy Wy Wy y+ (1 +9) Wyv 
ee aL Be (58) 
qGé Wuu Cpe Cuy 


In the U, V, W co-ordinate system (58) reads at the point under consideration 


dq us, a? WoutWyv (59) 
dt qv Wyv Wyv —Wi,’ 


where 7 is the distance measured along the radius. When (30) and (31) are used, 
the acceleration along the physical streamline becomes 


pom one! Ky BS, 80 
§s My Ke Mr (01 =fr 02) : (60) 


This acceleration is therefore seen to be simply related to the principal radii 
of curvature of the hodograph surface. 


From this result the pressure gradient along the streamline may be found: 


hy ies oa Ky Lod 
ds Mr Ke My ert @2), 


(61) 


where g is the density. The acceleration normal to the streamline may be derived 
from equation (45) and (60): 


so =— = [o.+ cal 48 (1 1+ (14 Fe) 1 (62) 


From this result follows the radius of curvature of the streamline in the physical 
space: 


R= tt . (63) 


| 1+ @2| Ea 
| 1 


Die 
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5. Analysis of the hodograph surface when the transformation is regular 


Once the co-ordinates at a point of the hodograph surface are given by the 
velocity g and the tangent plane at that point by the direction of the normal, the 
geometry of the surface may be further specified by the curvature, given by the 


Fig. 2a. Direction of the conical hodograph streamline 


value of the three second derivatives. The hodograph surface is described by the 
differential equation (18), which yields one relation for the second derivatives in 
two perpendicular directions, 7.e. (35). Two additional data are then required to 
determine a surface element to the second order. For these it will be convenient 
to choose @,/0, and @,+ 2. In addition the sign of a must be given, since 


without loss of generality «, may be chosen in the interval between — = and 5 ; 


When the velocity g and the tangent plane (determined by the direction of the 
normal or the radius in the physical space) are given, the value of M, is fixed. 
It may then be asked how @,/o, and 0,-+0, determine the surface. For a 
given value of M,, the directions of the characteristics are determined, and since 
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they are conjugate, to every value of 03/0, corresponds one value of the angle 
of the major axis a of the given sign. Thence follows also one value of B giving 
the direction of the conical hodograph streamline (or the acceleration), since it 
is the direction conjugate to the direction of the V axis. 


| ce; 
IU (Mea) t 
2 


Mp=V2 


a arctanyMZ-/" Mp=ye’ 


S 


1<Mp<Ve" 


a roton-Me- 1) 


MpoVF 


Mo co 


x 
a ( j2— jf)! 


Fig. 2b. Direction of the major principal axis at a point of the hodograph surface 


The parameters M, and 9,/0, therefore determine «, and 6. The parameter 
0,;+@2 may be seen from (60) and (62) to determine the magnitude of the ac- 
celeration and can be varied independently of 0,/0,*. 


With the aid of (20), (23), (32), (35) and (45) the following relations may 
then be derived for «, and f: 


mp1 
tang,=+ - — (64) 


* Except when 9,+0,=0; then @,/e,=—1. 
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a | G2 (me—1)b2—1— £2 
E | (m2—1)\Jmz—1— 8 
tan B =— I mastic as ee (65) 


Ate 123 
a 


where the -- signs are associated. These functions are illustrated in Fig. 2. 


In order to classify the possible motions of a gas particle in a conical flow, 
the shape of the hodograph surface at elliptic and hyperbolic points may now be 
investigated. For these points the Jacobians are finite and different from zero, 
so the transformation will be regular. In addition, parabolic points may be 
considered as limiting cases of these points when A,->0, and conical points 
and points on an edge surface as limiting cases when Kg—oo. The transfor- 
mation then becomes singular. Such points will be examined in more detail in 
later sections. 

It can be seen in Fig. 2 that M, and K, cannot be chosen completely inde- 
pendently of each other. This may be seen from (35) in the following manner. 
For elliptic points all radii of curvature at these points, and in particular Ry 
and R;, have the same sign. From (35) it then follows that M, >1, and the flow 
is conical-supersonic. Conversely, if the flow is conical-supersonic, it follows 
from equation (35) that Ry and Ry have the same sign; the point may then 
be either elliptic or hyperbolic, with the asymptotic directions lying in the same 
quadrant. If the flow is conical-subsonic (M,<1), Ry and Ry have opposite 
signs, the point is thus hyperbolic, with the asymptotic directions lying in 
different quadrants. Conversely, at a hyperbolic point, the flow is conical- 
subsonic (M,<1), conical-sonic (M,=1) or conical-supersonic (M,> 1), depending 
on the relative positions of the asymptotic directions with respect to the U axis. 


These results may also be obtained by forming the product of the two 
asymptotic directions. From (34) and (35) then follows* 


dV dV ss 
aac | arled 2 = as s ee 


The same conclusions may then be derived by noting that at an elliptic point 
the asymptotic directions are conjugate complex, while they are real at a hyper- 
bolic point. 


«) Conical-subsonic flow 

Consider first conical stagnation points. At such a point U=0 and M,=0, 
if a0 (M is finite). From (35) it follows that W,y-+W,y=0; thus the curve 
given by (22) for +C consists of two conjugate orthogonal hyperbolas, which 
are similar to Dupin’s indicatrix. From Fig. 2 it can be seen that 6 may have 
any value, so at the stagnation point conical streamlines from all directions 
can come together. This may also be concluded by remarking that at the point 
itself the direction of the U axis cannot be defined, since U=0. For the same 


* With the aid of equations (20) and (66) the following property may be derived: 


(0 ink (aes aaa av 
WON /asymp. 1 dU asymp, 2 Gane A (Fer en wy 0 x (66 a) 
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reasons the directions of the principal axes can also have any value. Since 
02/0,= —1, the acceleration along the streamline, as given by (60), is zero, since 
M=+-0. This result can be seen immediately by noting that the streamline falls 


along the radius through that point. 


Conversely, an orthogonal hyperbolic 
point does not necessarily represent a 
conical stagnation point. If the point is 
orthogonal hyperbolic, 0,/0,=—1, and 
the acceleration along the streamline is 
zero. The direction of the conical hodo- 
graph streamline is therefore perpendicu- 


lar to the U axis | = + Sie The conical 


streamline thus becomes perpendicular 
to its image, and the V axis therefore 
coincides with an asymptote. The U axis 
is then also an asymptote; thus Wyy= 
Wyy=0. From (35) it then follows that 
the conical Mach number M, may take 
any value, while for M,=£0 the principal 
directions bisect the angles between the 


U and V axes (m= +). From equa- 


tion (28) it is seen that at an orthogonal 
hyperbolic point the transformation be- 
comes conformal. 


For conical-subsonic flow (M, < 1) 
points on the hodograph surface are 
hyperbolic. From BuSEMANN’s relation, 
equation (35), it follows that for «,=0, 
02/0, has the value M?—1. Since the 
physical conical streamline falls along a 
principal direction, the hodograph coni- 
cal streamline also does so and is there- 
fore along the U axis; thus B=0. The 
acute angle between the asymptotes is 
bisected by the U axis. The situation 
is sketched in Fig.3a for accelerating 
flow and in Fig.3e for decelerating flow. 
According to (60), for accelerating flow 
0, + 02> 0, and, since | @,| >| Qo], it follows 
that 0, >0 and e,<0. From the equations 


(4), (+) V 
L PES 


if Ss 


4 


Bh 


e@ B=0;9<0 
Fig. 3a—e. Typical conditions at a hyperbolic point 
of a conical-subsonic flow 


of RopRIGUES, equation (27), it may be deduced that points on Dupin’s 
indicatrix ABCD are mapped onto a figure A’ B’C’D’ on the unit sphere 
which is compressed in the direction of AC and stretched in the direction of 
BD. The images of points in the directions of the asymptotes coincide with 
the image of the point considered. The figure is then turned over along the 
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major axis AC. For decelerating flow @,+@2<0 and since | @;| >| @s|, it follows 
that o,<0 and 9,>0. Dupin’s indicatrix ABCD is than compressed and 
stretched in the same manner, but turned over along the minor axis AC. In 
both cases the geodesic curvature of the physical conical streamline is zero. 
Further typical cases may be obtained when the physical conical streamline 
has a geodesic curvature different from zero and the acceleration along the 
streamline is positive, equal to zero or negative, respectively. Typical conditions 
at a hyperbolic point of a conical-subsonic flow are given in Fig. 3. The sketches 
are arranged so that the direction of the conical hodograph streamline turns 
in counter-clockwise sense from Fig.3a to Fig.3e. It may be noticed that 
Fig. 3e may be obtained from 3a by reversing the direction of flow, whereas 
Fig. 3b and Fig. 3d are also interchangeable in that manner. Reversing the 
flow direction in Fig. 3c amounts to changing the direction of the positive Y 
and V axes. From equations (64) and (65) and Fig. 2 it may be seen that for 
a given value of M,, as || is increased from zero on, |%| changes in such a 
way that |B|>]o,|. The acute angle between the asymptotes oscillates from 
2 arc tan |/1— M? for B=0 to 2/2 for an orthogonal hyperbolic point (|8| =2/2) 
(Fig. 3c). The conical streamline lies between the asymptotes which enclose 
the acute angle, except for orthogonal hyperbolic points, in which cases it 
touches one of the asymptotes, which then subtend an angle of 2/2 radians. 


From (43a), (43b) and (45) the variation of M, along a conical streamline 
may be obtained: 
a 


dU a 


y—1 y72 Kg | 
TE Sora teller 


or if by repeated use of (25) and (26) the second derivatives are expressed in 
terms of @,, @, and #,, and if furthermore (64) and (60) are used, we have 


dM, 1 avert wars 
(= 4 ones —4e Sy] ie 


When the flow is accelerating, so that g,>0, it follows that M. increases along 
the streamline if 
ace 1 2 Mr 
We (Ae arene 
and decreases if 
Wega a ‘| My , 
ai at aaa Mea?’ 
If the flow is decelerating, so that g,<0, it follows that M, increases along 
the streamline for 
/ Wx 2\ My 
Wei D) VE par oe 
and decreases for 


Meth aya Sly 
Wd sae dg 


The conical streamline is tangent to a line of constant M, for 


c 


Wires (1 ai > 1 M2) Mr , 
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When the acceleration is equal to zero, (43 b) shows that along the conical stream- 
line M, again may either increase, decrease or remain stationary. For comparison 
it may be noted that for two-dimensional plane flow the Mach number increases 
in expanding flow, decreases in compressing flow and does not change if the 
pressure does not vary along the streamline. 


B) Conical-sonic flow 

t a hyperbolic point the flow can also be conical-sonic (M,=1). From 
BUSEMANN’S relation, equation (35), it follows that if M,=1 and W,,+0, 
then Wyy=0. One of the asymptotes of Duptn’s indicatrix thus falls along 
the U axis. The other asymptote falls in the first (and third) quadrant for 
%, > 0 and in the second (and fourth) quadrant for #,<0. The angle of the major 
principal axis a and the angle of the conical hodograph streamline f are seen 
from (64) and (65) to satisfy the relations 


tan 04, == =e | —- Ge (68) 


and rere 


tan B =o Lae A (09) 


These relations are illustrated in the curve for M,=1 in Fig. 2. It may be 
deduced from equations (68) and (69) that if || is increased from zero onward, 
|x,| increases in such a way that 2|a,|>|8|>|o,|. The acute angle between 
the asymptotes equals 2|«,| and increases from zero to 2/2 when || increases 
from zero to 2/2, while |«,| increases from zero to 2/4. The conditions which 
may be encountered at a hyperbolic point if the flow is conical-sonic are similar 
to those in conical-subsonic flow, except that the case 6=0 cannot occur. This 
leaves three typical conditions, all with a curved conical physical streamline 
and positive, zero and negative acceleration along the streamline, respectively. 


They are sketched in Fig. 4. The case B= +=, c= =F _ again corresponds to 


an orthogonal hyperbolic point; thus W,,=0, and the acceleration along the 
streamline is equal to zero. It may be noted again that by reversing the direction 
of flow, Figs. 4a and 4c are interchangeable. It can further be shown that at 
a direction on the unit sphere perpendicular to the physical conical streamline 
the flow is conical-subsonic on the convex side of the streamline and conical- 
supersonic on the concave side. This direction, in fact, coincides with the Y 
axis and maps onto the U axis, since the latter is an asymptote and thus trans- 
forms perpendicularly to itself. On the Y axis the pressure on the convex side 
of the streamline is higher than on the concave side; the velocity gradient is 
therefore in the direction from the convex side to the concave side of the stream- 
line. From (39) it follows that U increases when q increases; thus U increases 
from the convex side to the concave side. From (43a) it may be seen that 
then M, also increases. The conical-subsonic flow is thus on the convex side, 
and the conical-supersonic flow region is on the concave side of the conical 
streamline. The positive U direction corresponds to the Y direction on the 
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Fig. 4a—c. Typical conditions at a hyperbolic point of a conical- 
sonic line 


concave side of the streamline 
and the negative U direction 
to the Y direction on the con- 
vex side. 

It follows from (44) that 
the direction of the conical- 
sonic line is given by 


cae 


Vata 1 


2 WW’ 


(70) 


or, from equation (30), since 
‘uu =0, by 

(Far) 

. dU }]M,=1 

eee oni 1 

GN A or aie 


(71) 


where the sign is chosen equal 
to the product of the signs 
of W,tanf and g,, since Wyy 
has the sign opposite to that 
of the product of tanf and 
g,. It can be remarked that 
if g, and tanf have the same 
signs (or the sign of tan is 
positive when g,=0), the coni- 
cal physical streamline has 
its concave side on the posi- 
tive Y axis. If g, and tanB 
have different signs (or the 
sign of tan # is negative when 
g,—0), the conical physical 
streamline has its concave 
side on the negative Y axis. 
If the concave side of the 
physical conical streamline is 
in the positive Y direction, 
the conical-sonic line on the 
hodograph surface therefore 
lies in the first (and third) 
quadrant for W>0 and in 
the second (and fourth) quad- 
rant for W<0, while it is 
tangent to the V axis for 
W=0 (then M=M,=1). If 
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the concave side of the physical streamline is in the negative Y direction, the 
conical-sonic line on the hodograph surface therefore lies in the first (and third) 
quadrant for W<0 and in the second (and fourth) quadrant for W>0, while 
again it is tangent to the V axis for W=0 (M,=M=1). Insummary, the angle 
on the hodograph surface between the velocity vector U and the part of the 
conical-sonic line on it, on the concave side of the conical physical streamline, 
is acute for W>0, equal to 2/2 for W=0 and obtuse for W<0. 

Since at a hyperbolic point Kg and also a//M?—1 will be finite, for all 
values of M, the conical-sonic line cannot be tangent to the U axis but may 
have any other direction. On the unit sphere the conical-sonic line therefore 
cannot be normal to the physical conical streamline. 


Again it may be deduced from (67) that, regardless of the value of the ac- 
celeration along the streamline, M, may increase, decrease or remain stationary 
along the conical streamline for a point on the conical-sonic line. This also 
follows from the result that the conical-sonic line may have any direction (except 
that normal to the conical physical streamline). The various possible directions 
of the conical-sonic line with respect to the conical streamline are illustrated in 
Fig. 4. 

The variation along the conical-sonic line of U,, the velocity component 
normal to the radius, is of interest because it shows a behavior different from 
that for two-dimensional plane flow. From equations (42a), (42b) and (70) it 
may be shown that for M,=1, 

dU, 1 
( dU laa se , (72) 


‘ 


thus dU,,/dU <0 for y=1.4. The velocity component normal to the radius increases 
along the conical-sonic line on the hodograph surface in that direction which 
makes an obtuse angle with the velocity vector U. In order to investigate the 
variation of U, along the conical-sonic line on the unit sphere, consider first 
the case g,=0, as sketched in Fig. 4b. It may then be shown that along the 
conical-sonic line on the unit sphere U, increases in the direction which makes 
an obtuse angle with the velocity vector U when W>0; that U, increases in 
the direction which makes an acute angle with the velocity vector U when 
W <0; while U, remains stationary when W=0. By considering the properties 
of the mapping when g,+0 (Figs. 4a and 4c), the same conclusions may be 
seen to hold. The first result was also given in [17], where it was tacitly assumed 
that W>0 (v,>0 in the notation of [17]). For the direction along the conical- 
sonic line in which U, increases, since U,=M, a=a, it follows that a increases; 
from (2) it follows that g decreases; thus W decreases, while from BERNOULLI'S 
law it follows that p increases. 

The conical hodograph characteristics, as given by (20), are both tangent 
to the U axis. Their geodesic curvature x, may be obtained by differentiating 
(19) along a characteristic. By use of (35), for the image on the hodograph 
surface of a point X= Y=0 we obtain 
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Fig. 5a—c. Typical conditions at a hyperbolic point of a conical-super- 
sonic flow 


where the + signs in equa- 
tions (73) and (20) are associ- 
ated and x, is positive if the 
concave side is on the positive 
V direction. At a conical-sonic 
point the curvatures of these 
characteristics thus approach 
+ oo. Both characteristics 
curve away from the U axis, 
and their images, the physical 
conical characteristics, both 
curve away from the Y axis. 
The physical conical charac- 
teristics lie on the concave 
side of the physical conical 
streamline, since this is the 
conical-supersonic region. 


y) Conical-supersonic flow 


At a hyperbolic point, coni- 
cal-supersonic flow can also 
occur. From (64) and (65) and 
from Fig. 2, conditions which 
may be encountered at such 
a point may be determined. 
The most characteristic of 
them are sketched in Fig. 5, 
and again 1t may be seen that 
reversal of the flow may be 
used to deduce Figs.5a and 
5c from each other. These con- 
ditions are largely similar to 
those at a hyperbolic point 
of a conical-sonic line, or a 
conical-subsonic flow, as dis- 
cussed above. It can be re- 
marked, however, that the 
conical hodograph streamline 
cannot lie in the region around 
the U axis enclosed by coni- 
cal hodograph characteristics, 
since |f|>|arc tan /M?— 4]. 
Also it follows from equations 
(64) and (65) that |B|/=|o,|, 
so that the major principal 
axis lies between the hodo- 
graph streamline and the U 
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axis. The direction of the conical hodograph characteristics with respect to the 
direction of the major axis depends on the conical Mach number M.. For M. < 2, 
it follows that |o,|=|arc tan /M?—41], and for M,= //2, |o|< | arc tan /M? — 1]. 
If M,—1, conditions at a conical-sonic line as discussed before are again found. 
The other limiting case occurs when M,->oo. 

Since U remains finite, a0 when M,—>oo, and since g remains finite, 
M=gq/a—oco. The angle of the major axis « and the angle of the conical hodograph 
streamline f are obtained from (64) and (65): 


a oer, 
tana = +4 |/— %, (74) 
Or 
and 
tan B > + 0. (75) 


These relations are illustrated by the curve for M,+o in Fig. 2. From 
equation (66) it is seen that one of the asymptotes falls along the V axis; thus 


W,y=0. The other asymptote then makes the angle + (2 |o,|— = with the 


U axis. Since g, and @, remain finite and a— +0, it is seen from (60) that the 
acceleration along the streamline is equal to zero. When M,-> oo, the acceleration 
normal to the streamline, according to equations (45), (62) and (65), is 


En = as / % ae (76) 
and the radius of curvature of the physical streamline is 
2 See 
R= Pe Ke (77) 


A point in a conical-supersonic flow may also be represented by an elliptic point 
on the hodograph surface. Typical conditions encountered at elliptic points again 
may be deduced from the relations given in (64) and (65), as illustrated in Fig. 2 
and sketched in Fig. 6. Again, it may be seen that Figs.6a and 6d, and 
Figs. 6b and 6c, may be obtained from each other through flow reversal. If 
B=0, the conical streamline remains parallel to itself when transformed and 
thus falls along a principal axis; thus « or % equals zero. According to equation 
(36) Ry/Ry has the value M? —1, so that for M_< 2 the major axis of Duprn’s 
indicatrix, which is an ellipse, falls along the U axis and the minor axis along 
the V axis, while for M>/2 the minor axis falls along the U axis and the 
major axis along the V axis. If Mq= /2, Dupry’s indicatrix is a circle, and 
the point is an umbilical point. The transformation is then conformal. In 
Fig. 6a the situation is sketched for accelerating flow and in Fig. 6d for deceler- 
ating flow. From equation (62) it may be concluded that for accelerating flow 
0,-+0.>0 and, since 9, and g, have the same sign, 9, >0 and 9,>0. From the 
Rodrigues equations (27) it is then seen that Dupin’s indicatrix 4 BCD when 
mapped onto the unit sphere is compressed in the direction of the major axis 
and stretched in the direction of the minor axis into the ellipse A’ B’C’D’, 
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Fig. 6a—d, Typical conditions at an elliptic point of a conical- 


supersonic flow 


which has the same ratio of princi- 
pal axes as ABCD. If the flow 
is decelerating, 0,-+0:<0, and, 
since 9, and go, have the same 
sign, it follows that 0,<0 and 
0,<0. DuPIN’s indicatrix is com- 
pressed and stretched in the same 
manner and then rotated in its 
plane through 180°. Other situations 
that may be encountered at an 
elliptic point are given in the other 
sketches of Fig. 6. It may be noticed 
that the conical hodograph stream- 
line cannot lie in the region inclosed 
by the conical hodograph charac- 
teristics around a line through the 
point considered and normal to the U 
axis, since |B|<|are tan /M? —4|*. 
For finite values of M. the accelera- 
tion along the streamline is there- 
fore different from zero at an elliptic 


point, since then || +=. In con- 


trast to the situation for conical- 
supersonic flow at a _ hyperbolic 
point, at an elliptic point |B| < |o|, 
while furthermore for M <J2, 
|a,|<|are tan /M?—1| and for 
M,=YV2, |o,|=|are tan /M?—1]. 
The case M,->co does not occur 
at an elliptic point, since, as may 
be seen from equation (74), «% be- 
comes imaginary for 9,/0.>0. 


* The difference between elliptic 
and hyperbolic points in regard to the 
direction of the hodograph streamline 
with respect to the hodograph charac- 
teristics may also be understood in the 
following way. The conical disturb- 
ances in the flow travel along the 
downstream physical characteristics 
which map onto those parts of the 
hodograph characteristics that are bi- 
sected by the U axis for an elliptic 
point and by a line normal to the U 
axis and through the point considered 
for a hyperbolic point. In order for 
these parts to be downstream charac- 
teristics the hodograph streamline 
should lie in the regions described. 
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Singular points of the transformation may now be considered as limiting 
cases of elliptic or hyperbolic points as.0,/0,->0. In these cases the direction 
of the conical hodograph streamline and the major axis approach one of the 
directions of the conical hodograph characteristics. Singularities thus occur for 
M,=1. Depending on how g,+@, is chosen to vary in this limiting process, 
different singularities may be obtained. If 0,0, remains finite, the acceleration 
remains finite, g,—>0 when @,/0,->0, while @, remains finite. Thus the Gaussian 


> : 1 " : 
curvature Kg—>-+oo, since Kg = a A region on the unit sphere corresponds 
1&2 


to a conical point or an edge surface, which represent a region of parallel flow 
or a conical simple wave flow, respectively. If 0,+0,—-ce, the acceleration 
approaches infinity, and @,—>co when 9,/0,->0. Thus the Gaussian curvature 
Kg vanishes. The point on the hodograph surface is a parabolic point, 
representing a point on a limit cone or the edge point in the flow around a 
sharp edge. 


These singularities will now be discussed in more detail in the following 
sections. 


6. Limit cones or conical limit lines 


Limit lines or surfaces probably were first discovered in some solutions of 
the hodograph equation for two-dimensional plane flow. A more systematic 
investigation of the properties of limit surfaces has since been given for two- 
dimensional and three-dimensional flow. For an extensive discussion, giving many 
references to the literature, reference may be made to [28]. 


If limit lines or surfaces appear in a solution, regions in the flow are found 
for which the velocity is many-valued. The transformation from the physical 
space into the hodograph space therefore becomes singular. Regions with a 
many-valued solution for the velocity are bounded by limit surfaces, so called 
because the direction of the flow is reversed at these surfaces and the flow thus 
has a limiting boundary which cannot be crossed. Two types of limit surfaces 
may be distinguished. 

For limit surfaces of the kind most studied, the Jacobian determinant is 
assumed to be continuous and to vanish at this surface. In the reversal of the 
flow the acceleration and the pressure gradient then go to infinity. In addition, 
the limit surface in the physical space is the envelope of the characteristic 
surfaces of one family, while, correspondingly, in the hodograph space the 
streamlines become tangent to the characteristic surface of the other family. 
The other type of limit surfaces occurs when the Jacobian determinant changes 
sign discontinuously across a characteristic surface. The acceleration and the 
pressure gradient are discontinuous but remain finite, and no envelope of 
characteristic surfaces is formed, while accordingly the streamline in the hodo- 
graph space is not tangent to the characteristic surface. 

This physically unacceptable behavior of the flow at limit surfaces can appear 
in the solution because it is assumed that the coefficients of viscosity and heat 
conduction of the gas are zero. It may then be assumed that thermodynamical 
processes in the gas are reversible and, if no heat is added, isentropic. In the 
absence of frictional forces the inertia and pressure forces controlling the motion 
of the gas remain in balance when the direction of the flow is reversed, so that 
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to every isentropic solution with some prescribed stream surfaces there corresponds 
4 solution with reversed flow. At limit surfaces, however, the reversibility of 
the flow occurs in such a way as to make obvious its physical impossibility, 
since it results in a multivalued region for the velocity. The occurence of infinite 
velocity gradients in the first type of limit surfaces serves as an indication that 
viscous stresses may no longer be neglected and that the assumption of isentropic 
flow, implying reversibility along the streamlines, is not justified any longer. 
The phenomenon of limit surfaces is analogous to ocean waves breaking on the 
beach, where the continued steepening of the waves is not counteracted by a 
mechanism analogous to the action of viscous stresses in the gas flow. In order 
to obtain solutions of physical value, layers called shock waves have to be 
introduced, these being layers in which there occur large velocity gradients and 
hence considerable effects of viscosity and heat conduction. If in such solutions 
the viscosity is assumed to approach zero, the thickness of the shock waves 
approaches zero, and the shock waves may be considered to be surfaces of dis- 
continuity connecting inviscid solutions. 


One is thus led to the suggestion that to every solution with a shock wave 
there corresponds a solution with limit surfaces along which compression occurs. 
If the direction of flow in these solutions is reversed, expansion occurs along 
the limit surfaces, and these solutions do not have physical significance, since 
expansion shocks cannot be formed. To every solution containing shock waves 
then corresponds a solution with reversed streamlines without physical meaning. 
It should be noted that it is not possible to state in general that to every solution 
with limit surfaces where compression occurs there corresponds a solution with 
shock waves. 


Apart from the trivial case of parallel flow throughout the physical space, 
a stream surface in a supersonic flow may be taken to represent the surface 
of a body immersed in a supersonic stream, which, in general, experiences wave 
drag. Since wave drag is associated with the entropy rise through shock waves, 
shock waves occur in all flows which are supersonic in the sense taken above, 
and the corresponding isentropic flows contain limit surfaces. The singular 
behavior of the pressure waves along these surfaces appears to be an essential 
feature of non-linear isentropic supersonic flow. 


In conical flow a limit surface is necessarily conical, so that we may speak 
of limit cones. The intersection of a limit cone with the unit sphere will be 
called a conical limit line. As will be shown now, lines of parabolic points on 
the hodograph surface in general represent conical limit lines of the first type, 
while the second type of conical limit lines may occur along a conical character- 
istic. 

a) Conical limit lines of the first type 


Consider parabolic points as limiting cases of elliptic or hyperbolic points 
by letting @,->0o, thus 0,/0,>+0. From equations (64) and (65) it then follows 
that the direction of the conical hodograph streamline and the major axis 
approach one of the directions of the conical characteristics. One property of 
a point of a limit line, namely, that the hodograph streamline is tangent to a 
characteristic there, is thus seen to hold also at such a parabolic point. 
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The two asymptotic directions coincide with the major principal direction, 
and DupIn’s indicatrix consists in two parallel lines. All directions different 
from the major principal direction are conjugate to the latter and thus map 
onto the unit sphere along the direction of the minor principal axis. The hodo- 
graph characteristic which is not tangent to the major principal axis therefore 
also maps into this direction. In order to determine how curves on the hodograph 
surface which are tangent to the major principal axis are transformed, it is 
necessary to consider the third derivatives. Two relations between the four 
third derivatives at a point of the hodograph surface are given by differentiation 
of the differential equation (18). 


By partial differentiation with respect to u and v in a co-ordinate system 
where the w axis is taken in the direction of the normal at the point considered 
(that is, the direction of the radius in the physical space), for the third derivatives 
at that point there follow the expressions 
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It appears convenient for present purposes to rotate the axes around the w axis, 
so that the w axis becomes parallel to the major axis and the v axis parallel 
to the minor axis at the point considered, and u>0 at this point. Then, since 


at a parabolic point w,,—=w,,=0, Wy,=——, u/a=1 and vja= + VMP—1, 
(78) and (79) yield Qe 
— (i — 2) yuu + 2 MP1 yyy + ae fa (78a) 
and 
EMA 2) cach 2 VI 10, 4. ME (79a) 


Let 'the co-ordinates at a point be given by u=u,, v=v, and w=w,. A Taylor 
expansion including terms up to the third order then gives at such a point 
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Consider a curve at (u#,,v,,#@,) making an angle a, with the major principal 
axis and having a geodesic curvature x,. The equation of such a curve may 
be written as 
vy — v,= tan a, (4 — u,) + 3%, cos Fa, (u — m4)? +--+. (81) 
Arch. Rational Mech. Anal., Vol. 6 23 
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If this expression is put into equation (80) after the first derivatives of w are 
taken, the image of this curve on the unit sphere may be seen to subtend with 
the major principal direction an angle «», given by 
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If 9, and g, are both finite and different from zero, 7.e., at elliptic and hyperbolic 
points, (28) is found again when the limit is performed. At a parabolic point 
0,->00; thus when «, == 0, equation (82) shows that all directions different from 
the major principal direction map into the direction of the minor principal axis, 
as we have noticed above. When «,=0, that is, when a curve on the hodograph 
surface is tangent to the major axis, (82) yields 


4 
ee 
= ieee 
tan a), = 2. 
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(33) 


From equations (73), (45) and tan B= + M2? —1, the geodesic curvature of the 
hodograph characteristics in a parabolic point is found to be 


Ce \née = ee (84) 
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which is equal to the expression for two-dimensional plane flow if M, is replaced 
by M. 

By use of this result and (78a) in (83) the direction of the image of the 
hodograph characteristic tangent to the major principal axis is found to be 
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which when compared with equation (10) is seen to be the direction of the corres- 
ponding conical physical characteristic. The hodograph characteristics thus map 
onto the physical characteristics. 

It may also be seen that the conical hodograph streamline maps onto the 
conical physical streamline. By differentiation of (50) along the streamline it 
follows in the above-mentioned co-ordinate system that 
y3 
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from which the geodesic curvature of a hodograph streamline in a parabolic 
point may be obtained: 
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Using this expression in (83) yields the correct value (tan Oy) s= m2 —41. 
If wy,.=-0, the return of a conical physical streamline at an image point of a 
parabolic point may be deduced. In this case, since at a parabolic point the 
conical hodograph streamline is in the direction of the wu axis and w,,,,—0, the 
projection of the hodograph streamline on a normal plane through the major 
principal direction has a point of inflection. Since furthermore (t,),=0, the 
normal to the hodograph surface along a streamline thus has an extreme value 
at the parabolic point, and the conical physical streamline returns and exhibits 
a cusp. The velocity at a parabolic point changes sign, so that downstream 
of the parabolic point the solution on the surface w(u,v) continues on the 
surface —w(—wu,—v). The rays in the physical space upstream and down- 
stream of a line of parabolic points cover the same region, so a many-valued 
region on the unit sphere appears, bounded by the image points of the parabolic 
points. This may also be seen by noting that differentiation of (30) along a 
streamline at a parabolic point yields 
dk Oates , 
(Fa een S 

uuu 0. A parabolic line then separates a region 
of elliptic points from a region of hyperbolic points. A closed curve on the 
hodograph surface which does not intersect itself and which is divided into 
two parts by a parabolic line is traversed along its image on the unit sphere 
in the same direction as on the hodograph surface in the region of elliptic points, 
but in the opposite direction in the region of hyperbolic points. A doubly covered 
region on the unit sphere thus exists. All curves tangent to a major principal 
axis, including one of the two families of characteristics and the streamlines, 
when mapped onto the unit sphere return and show a cusp at the image of the 
parabolic line. All curves which are not tangent to a major principal axis, in- 
cluding the other family of characteristics, are tangent to this image line. The 
image line is the envelope of these curves and thus also, for example is the 
envelope of the lines of constant speed and of one family of physical conical 
characteristics. The latter may also be seen by transforming the physical conical 
characteristics, as given by (8), to the hodograph surface by means of equa- 
tion (15). The equation which follows may be written at the point considered 
in the form 7 \2 A 
(G7) — 2 — 1)] =o. (89) 

The conical hodograph characteristics are therefore the representation of the 
conical physical characteristics, as might be seen by equating to zero the ex- 
pression in brackets in (89) and then comparing the results with (20). Further, 
(89) may be satisfied along a line on the hodograph surface where K;=0. This 
‘parabolic line, when mapped to the unit sphere, will have at any point a slope 
equal to that of one family of physical conical characteristics. Since this line 
does not coincide with a characteristic, it will be the envelope of the physical 
characteristics of one family. 

By differentiating equation (30), we may find the slope of the parabolic line, 
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which is different from zero if w, 
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so that for w,,,,--0 and w,,,, finite the parabolic line does not touch a major 
principal axis at any point. 

As has been remarked previously, by equations (60), (61) and (62) the ac- 
celeration and pressure gradient go to infinity, since 9,00. In addition it is 
concluded from (63) that the radius of curvature of the spatial streamline goes 
to zero. 

It is thus shown that the flow at a parabolic point, if w,,,,,==0, shows the 
typical behavior exhibited at limit surfaces of the first type, and that parabolic 
lines are the images of limit cones or conical limit lines. 

The variation of M, along the streamline at a parabolic point is related to 
the sign of the acceleration. From (67) it follows that for Ag=0 


= fr Psa 0» 


Since this quantity is positive, M, increases in an accelerating flow (g,>0) and 
decreases in a decelerating flow (g,<0). This is at variance with the situation 
at elliptic and hyperbolic points, where the increase or decrease of M, along 
the streamline is not determined by the sign of the acceleration alone. On the 
contrary, as with the geodesic curvature of the characteristics, equation (84), 
the local flatness of the hodograph surface (Kg=0) is again expressed, since 
if M, replaced by M, the same expression for the variation of M is found as in 
two-dimensional plane flow. 

Conditions which may be encountered at a parabolic point when w,,,,,=0 
are sketched in Fig.7*. Again Fig. 7d may be obtained from Fig. 7a and 
Fig. 7c from Fig. 7b by flow reversal. Further it is seen that in an expanding 
flow the conical limit line is inclined along the upstream conical physical char- 
acteristic on the concave side of the streamline or along the downstream conical 
physical characteristic on the convex side of the streamline. If the flow is com- 
pressing, the conical limit line is inclined along the downstream conical physical 
characteristic on the concave side of the streamline or along the upstream conical 
physical characteristic on the convex side of the streamline. 


If B=0, as for example in axially symmetric conical flow, it follows at once 


that M,=1. From (78a) and (79a) it follows that then w,,,,—=—27' and 
Wyyuy=O. From (88) it is seen that a Q2 

dK¢ Ad y+1 

Ee ie Rent, Mick aon? (92) 


so that the hyperbolic points are on the side of the origin (U=V=0) with 
respect to the parabolic point, and the elliptic points are on the other side. 
For accelerating motion the flow changes from conical-subsonic to conical- 
supersonic flow, whereas for decelerating motion the flow is from conical- : 
supersonic flow to conical-subsonic flow. From equation (90) it is seen that 
the parabolic line is perpendicular to the streamline (6=+42), as is the conical- 
sonic line, as may be concluded from (44). On the unit sphere the conical limit 
line is also perpendicular to the streamline. If the flow is axially symmetric, 


* The co-ordinate axes X, Y and Z, and U, V and W are defined by the velocity 
components on the upstream side of the parabolic point. 
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the two lines coincide, so that each conical limit line in an axi-symmetric conical 
flow, which is represented by a parabolic line on the hodograph surface, is a 
conical-sonic limit line. 


aaa ely (Ny aU, D 


@-B=orctantye—1 
“A U 


Y Ke>0 
Chor.! 
= TU(-U,-v) 
: eang 
b B#0;9570 conical limit line 
—TMY(-U,-V) 
+)V, 


U 
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d B=0;9,<0 conical limit line 


Fig. 7a—d. Typical conditions at a parabolic point of a conical limit line 


Singularities on conical limit lines may now be investigated by considering 
third and higher derivatives at parabolic points. It may be expected that 
because of the planar character of the hodograph surface at parabolic points 
these singularities will bear a marked resemblance to those of limit lines in two- 
dimensional plane flow. 
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For example, if the third derivatives are taken into account, just as in 
two-dimensional plane flow cusps of a conical limit line may be encountered. 
The slope of the parabolic line may be taken as the parameter which governs 
conditions at the point in this instance. Expressing the third derivatives of 
interest in terms of this quantity 6 (where 6 is the angle with the major principal 
axis, measured positive in the counterclockwise sense), by means of equations 


(78a), (79a) and (90) we have 


A ia yah tan 6 (93) 
ig 40, —(M2—2)tand+2VM2—1 ” 
Wy “uv pals s z: +2 (94) 
40, — (M2—2)tand+2)/M2—1 
1 _ yl — 
2 iS 2 Ss a= Mes 
0) = Vad ee s) 2) ( tag | act ae “2 : (95) 
mee Os —(M2—2)tand+2VM2—1 : 


where the upper and lower signs are to be taken according as the streamline is 
inclined at the positive or negative characteristic angle. 

Accordingly, the geodesic curvature of the hodograph streamline and the 
variation of the Gaussian curvature along the streamline may also be expressed 
in terms of 6, so that with the aid of (87) and (88) may be obtained the formulae 


: a) + VM2—1 tan 6+1 
(%).=7 a se 
a —(M?2—2)tand+2)M?2—1 
ae _ ytt tan 6 (97) 
du ]Kg=0,s 403 —(M2—2)tand+2/M2—1_ 


It may now be seen that if the parabolic line is tangent to the hodograph stream- 
line, when 0=O the conical limit line exhibits a cusp. From (97) it appears 
conicol limit line that the Gaussian curva- 
S thore ture remains stationary 
ae / along the streamline in this 
@=B-orctonty42-; case, and from equation 
(93), that w,,,—0. If then 
Wuuuu 0, the correspond- 
ing physical streamline does 
not return. If it is assumed 
that there exists on the 
surface where w,,,,,= 0a line which does not touch the streamline (since wy, ,=-0), 
the streamline through the considered parabolic point separates a region with 
streamlines which do not intersect the parabolic line from a region with stream- 
lines which intersect the parabolic line twice. Correspondingly, the physical 
streamlines either do not intersect the conical limit line and do not return, or 
reverse twice, at the first and second branches of the limit line successively. A 
sketch of the conditions at a cusp of the limit line is given in Fig. 8. In order 
for the characteristics to map in the way indicated in this figure, it may be 
noted that the hodograph streamline should osculate the hodograph character- 
istic. From (84) and (96) it may be shown that for 6=0 the geodesic curvature 
of the streamline and its tangent characteristic are equal. 


Fig. 8. Cusp at a conical limit line 
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For given values of y, a, M, and g., the quantity 6 uniquely determines 
the third derivatives w,,,,, Wy,, and w,,,, unless the determinant of the system 
of equations (78a), (79a) and (90) vanishes. This determinant is equal to 


+ 2M? —1[— (M?— 2) tan 6 + 2M? — 1]; 


thus it is equal to zero for M,=1 and for tan d= +2/M?—1[M?— 2]. The 
Wate 
4 Q2 
and w,,,,—0. In the latter case the parabolic line is tangent to the other hodo- 


graph characteristic, and the system becomes dependent if 0, ->0v. In accordance 
with the theory of characteristics, it is possible to assign not only 6 and w 
but also one of the other third derivatives. 


= s Dp IU 
case M,=1 was discussed above, where it was shown that 6= 52. Cave 


VVV 


Among the singularities of higher order, the parabolic line that coincides 
entirely with a streamline is of particular interest. This line is a plane curve 
which coincides with a conical characteristic. At all points of it the normal 
to the hodograph surface is also the binormal to the plane curve, so that one 
point on the unit sphere corresponds to a curve on the hodograph surface. 
Since @,->oo, the accelerations along and normal to the streamline are infinite, 
as may be deduced from (60) and (62), while the radius of curvature of the 
physical streamline is zero. The flow is one around a sharp edge, and conditions 
at the edge point are given by the parabolic line on the hodograph surface. 
Such a singularity may be met for example in the flow around a subsonic leading 
edge of a flat delta wing at incidence or in the flow around the tip of a flat 
rectangular wing at an angle of attack with respect to a parallel supersonic flow. 
This type of singularity is different from the actual conical limit lines, since 
no flow reversal occurs, but is mentioned here because, as in the case of conical 
limit lines, the Jacobian of the transformation vanishes. 


B) Conical limit lines of the second type 

In order to recognize the occurence of conical limit lines of the second type 
we now consider how the conical hodograph characteristics determine the 
geometry of an integral surface in the hodograph space. 

Let a region R of an integral surface be given, including the boundary C; 
we ask how this surface can be extended into a region S, adjacent to F& along 
C, in such a way that w(w, v), w,(u, v) and w,(u, v) are continuous across C. In 
the region S the co-ordinates of the points and the tangent planes to the surface 
along C are thus known. Furthermore, the shape of C, being the boundary of 
the region R, is known, and since C is also a curve in S, the geodesic torsion 
and the normal curvature of C, being the two geometric properties determined 
by the second derivatives in S, are known. In addition, S is known to satisfy 
the differential equation (18). Written in the U, V and W co-ordinate system, 
the following equations are thus obtained: 


Woy + (1 — M?) Gv =0, (98) 
— sin « cos «a Wyy + (cos? « — sin?a) Wyy + sina cosaWyy=T,, (99) 


cos?a Wyy + 2sina cosa Wyy-+ sintaWyy = — x,, (100) 
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where t, and %, are the geodesic torsion and the normal curvature of the curve 
C on the surface, respectively, making an angle « (« measured positive in the 
counterclockwise direction) with the U axis. 

Two cases may now occur; either the rank of the coefficient matrix is equal 
to 2 or equal to 3, while the rank cannot be equal to 1. In the first case equations 
(98) — (100) may be solved for the second derivatives in S to obtain the same values 
as those on C in R. Secondly, if «=arc tan VM 4, the rank is 2, and the 
equations are dependent if there is a linear dependency between T, and %,. 
In this case one of the second derivatives in S may be chosen, so that a jump 
in the curvature occurs along C, which is then a characteristic. 


For the geodesic torsion of a characteristic then follows the formula 
(T,)ehar. = are [Wov + VME — 1 Wey], (101) 
and for the normal curvature of the characteristic, 


2VM2—1 


ue Wov + VME — 1 Wry), (102) 


(7, Jehan: =F 


where upper and lower signs in equations (101), (102) and (20) correspond. 
From this dependence of (T,)char, ANd (%,)char, it follows that 


M?2—2 
a bie =+ 2 ae (ae eter. . (103) 
_ fe) =53 


Geometrically, this result may be understood by noting that the conical 
characteristics form a conjugate net on the hodograph surface. If one travels 
along one characteristic, the tangent plane thus rotates around the other char- 
acteristic direction. The geodesic torsion of the conical characteristic is there- 
fore determined by its shape and by the normal to the surface at the point 
considered, since these fix the conical Mach number M, and (%,) char. - 

Also, the geodesic curvature of the characteristic may be expressed in this 
way by means of (73) and (102): 

tt yt WN % 
(%g)char, = xb 2aVMEo ae as V1 (%y)char. (104) 

Ifa jump in Wyy across C is indicated by {Wy}, where {Wy} = (Wry) 5 — (Wrv)p 
(the indices S and R refer to conditions on C in S and R, respectively), the 
discontinuities in Wyy and W,, are obtained from (98) and (101), as follows: 


{Wov} = (M? — 1) {Wy}, 
(Wov} = = VM? —4 {Wy}, (105) 
{Wry} = {Wy}. 


From these results the discontinuity across C of the Gaussian curvature may 
be calculated, and it is found that 


{Ke} =42/Me—1 [Wov + VM2—1 Wry] {Wry}, 
or 


, (106) 
{Kg} =—M? (%») char. (Wpr}. 
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By differentiation of the differential equation (18) on both sides of C, the 
jump across C may be seen to obey along C a differential equation governing 
the propagation of the jump disturbance along this characteristic. If, then, 
somewhere on C a discontinuity {A,} is introduced which propagates along C 
in such a way that a part of C exists where the Gaussian curvature changes 
sign across C, this part is a curve on the surface separating a region of elliptic 
points from a region of hyperbolic points. There is then a doubly covered region 
on the unit sphere, bounded by the physical conical characteristic C’, which 
is the image of C. The acceleration has a discontinuity across C but remains 
finite as long as the principal radii of curvature remain finite. Accordingly, 
the physical streamline returns and has a discontinuity in curvature at the 
characteristic C’ where it has a cusp. The conical hodograph streamlines and 
the lines of curvature on the hodograph surface have corners at C. 

It might possibly be conjectured that the different nature of the two types of 
limit lines corresponds to whether they lie in or on the boundary of an independent 


hyperbolic region (second type of limit line) or a hyperbolic region depending on an 
elliptic region (first type of limit line). 


7. Conical simple wave flow 


Spatial simple wave flows may be defined as flows where the velocity vector 
and the state of the gas do not vary in planes tangent to the characteristic 
surfaces originating at points of these planes. If all these planes go through 
at least one common point, the flow is conical, and the point of intersection 
is the center of the conical field*. The intersections with the unit sphere are 
great circles which are conical characteristics, since each plane is the envelope 
of characteristic surfaces originating on rays in it. The name “‘simple waves’’ 
refers to the fact that the conical characteristics of the other family do not 
carry flow disturbances, because the state of motion of the gas remains constant 
when crossing these conical characteristics along a straight characteristic of the 
first family. 

In order to investigate the behavior of conical simple wave flow in the 
hodograph space, consider first an integral surface R with boundary B. Let 
this surface be divided into three regions R,, S and R, by the characteristics 
C, and C, in a manner as sketched in Fig. 9. The image of & on the unit sphere 
is called R’, while all curves and points on R’ corresponding to their images 
on R are denoted by primes. In order to correspond to physical reality, as 
given by the underlying assumptions of isentropic flow throughout the region 
bounded by B’, it may be assumed that w(u,v), w,(#,v) and w,(w,v) are 
continuous functions on R. The region S may consist of elliptic, parabolic and 
hyperbolic points, and characteristics along which discontinuities in the second 
or higher derivatives occur. It is thus assumed that the Gaussian curvature 
remains finite in S. Chose a point P on C,, let the conical hodograph charac- 


teristic of the other family through P be C, and the streamline s. Construct 


* Spatial simple wave flows which are not conical are excluded from the present 
discussion. An example of such a flow is the isentropic supersonic flow over a curved 
two-dimensional profile in a homogeneous stream. 
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a plane N at P normal to C, and at P' a plane N’ parallel to N. Since the 
normal at P lies in N, O’ lies in N’, and N’ cuts the unit sphere along a great 
circle, which is denoted by P’Q’, while its image on & is denoted by PQ. The 
image of C, on the unit sphere is Cy, and since*at"™P" the line Po is normal 
to C,, Cy is tangent to P’Q’ at_P’. Call U the intersection of C, and Ce, and 
call its image U’. The region S’ may now be considered to become a simple 
wave region if discontinuities in the Gaussian curvature are introduced along 
C, and C, and if R, and S are deformed in such a way that U’ approaches Q’, 
while C; coincides with P’ Q’ if the radius of curvature of the arc PQ approaches 
zero when Q coincides with P, wherever the point P is chosen on Ci. inthis 
process the normals to the surface along C, are chosen to vary in such a way 


Fig. 9. Conical simple wave flow 


that they do not coincide with the normals along C, when C, and C, have come 
together. According to the definition of a conical simple wave flow given above, 
in this way S’ becomes a region with straight characteristics along which the 
velocity vector and the state of the gas are constant. In the hodograph space 
the simple wave flow is then given by a surface S of zero breadth, bounded 
along its length by two conical characteristics C, and C, and in general imbedded 
between the regions R, and R,. Such a surface will be called an edge surface. 
It follows from EuLER’s theorem, equation (24), that if the radius of curvature 
in any direction approaches zero, one of the principal radii of curvature must 
approach zero. On the edge surface S, therefore, 0,>0, and the Gaussian 
curvature goes to infinity. Duprn’s indicatrix of points on S degenerates into 
two isolated points with co-ordinates -+ V| o,| on the major principal axis and 
the point midway between of them, representing )|R|=0. Points on the edge 
surface may be distinguished by their normals, and from the limiting process 
it follows that at all points having the same velocity vector gq the normals lie 
in a plane normal to C; (or C,)*. On S, the images of different physical stream- 
lines in S’ flow together over the same space curve C, (or C,); they have corners 
where they leave or enter S from R, or R,. As has been noted in the previous 


* It may be clear by now that it is not an arbitrary curve on R that may be 
chosen as a sharp edge in the hodograph surface, where w, and w, are discontinuous. 
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discussion of the transformation for finite values of the Gaussian curvature, 
when 0./0,—>0, the major principal axis falls along the direction of one of the 
conical characteristics, which then also is the direction of the conical streamline. 
From (32) it may be seen that when 9,=0, if two directions are conjugate, at 
least one of them coincides with the major principal direction, so all directions 
are conjugate to the major principal axis. Thus all lines having such a direction, 
including the conical hodograph characteristic which is not tangent to the major 
principal axis and the line of constant speed, map on the unit sphere into the 
direction of the minor principal axis, which then also is the direction of the 
straight conical characteristic. In order to determine the transformation of a 
curve tangent to the major principal direction the derivation of equation (28) 
may be repeated in a more general form. Using (15) then yields the relation 
Dypaic Opn eels, 


tan %,, = — 
ee Wy yt Wy» tan a ” 


(107) 


where the w axis again is chosen parallel to the normal to the surface, but the 
wu and v axes are still arbitrary. If they are chosen parallel to the principal 
axes, (28) is found again. If the w axis is chosen in the direction of particular 
interest, tan «,=0, and equation (107) together with (25) and (26) yield the 


expression 
Ta 


whan copes (108) 


un 


Therefore, since the geodesic torsion of the lines of curvature is zero, the direction 
of a tangent to these lines transforms parallel to itself. 


The direction of the image of the characteristic tangent to the major principal 
axis can be obtained by the aid of (103) and (108): 


Te ge SS ee (109) 


so again we see that the hodograph characteristics map onto the physical 
characteristics. 


In order to investigate the transformation of the streamline its geodesic 
torsion has to be determined. By repeated use of equations (25) and (26) along 
the principal directions, the second derivatives may be expressed as a function 
of 0,, 2 and %. With these results and equations (64) and (65) the geodesic 
torsion of the hodograph streamline may be obtained from (26) as 


diy | eee (Mg—1)\ 4M—1— a 


(110) 


so that, when 9,=0, it follows that 


(,),= ++ VMZ—1. (111) 
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If this value is used in equation (108), and if further it is noticed that for 


02=0, (%)s= ' the direction of the image on the unit sphere of the hodograph 


streamline turns out to satisfy 
tan a», = + VM? —1, (4112) 


which is equal to the expression for the conical physical streamline. The conical 
streamlines on the unit sphere and on the edge surface thus map onto each other. 


In order to carry over the reasoning for points on the surface with finite 
Gaussian curvature to points on the edge surface these points may be distinguished 
in a similar fashion as elliptic, parabolic and hyperbolic points. It should be 
noted that then the edge surface coincides with a space curve. Taking the tangent, 
principal normal and the binormal of this space curve as a reference system, 
the normals to the edge surface at points with the same radius vector g are seen 
to lie in the normal plane of the space curve. They can be given by the angle 
@ with the binormal at the corresponding point of the space curve, where 3 
may be chosen as measured from the direction of the binormal, increasing when 
turning to the principal normal over the shortest angle. According to MEUSNIER’S 
theorem, a curve on the edge surface in the major principal direction has 
geodesic curvature and normal curvature given by 

x, =— cost, (113) 
and f 

4,=— x sin ?, (114) 
where @ is the radius of curvature of the space curve. Since this curve also is 
a characteristic, application of (113) and (114) for 0= $a and use of (104) gives 


x 
ta rae 5 5 
G at (115) 
where q,, is the component of g along the principal normal (g, <0 since 9 =0). 
For the major principal radius of curvature of a point on the edge surface 

it follows from equations (114) and (115) that 


— 2 Gn ) 
yt sine tte) 
If now a closed curve on the edge surface is traversed, the end points of the 
unit normal vectors to the surface describe a curve which is traversed in the 
same direction as its image curve on the unit sphere if g,>0 and in the opposite 
direction if @,<0 in the enclosed region. Accordingly, points on this surface 
will be called elliptic edge points if e,>0 or hyperbolic edge points if 0,<0. 
If 9,> 0%, the point will be called a parabolic edge point. From (116) it follows 
that poimts for which 7<#<2z are elliptic edge points, and those for which 
0<0<a are hyperbolic edge points, while for parabolic edge points #=0 or a. 
It may be superfluous to remark that the edge surface does not necessarily 
extend over the entire range of #. 
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The acceleration along the streamline follows from equations (60) and (116): 


2s TNE 
Pet lie Feng) Silane 


Ree (117) 
so that at elliptic edge points the flow is accelerating, and at hyperbolic edge 
points the flow is decelerating. As may easily be seen from the transformation, 
in an accelerating conical simple wave flow the straight conical physical char- 
acteristics are inclined in the upstream direction on the concave side of the 
streamline and in the downstream direction on the convex side of the stream- 
line, while for a decelerating flow the straight conical characteristics are inclined 
downstream on the concave side of the streamline or upstream on the convex 
side of the streamline. 


The transformation from the unit sphere to the edge surface exhibits a 
singularity at a parabolic edge point. Two lines of parabolic edge points (#=0 
and #=z) may exist on the edge surface, separating a region of elliptic edge 
points from a region of hyperbolic edge points. The geodesic torsion of a line 
of parabolic edge points is by definition equal to the torsion of the space curve 
which coincides with the edge surface. 


Consider first the case when tT, the torsion of the space curve, is not equal to 
zero. Since @,—>0ce at a parabolic edge point, it follows from equations (111) 
and (103) that the geodesic torsions of the hodograph streamline and of the 
characteristic in the streamline direction are zero. When crossing the line of 
parabolic edge points these curves thus enter a region with a different sign of 
o, and therefore go in either direction from a region of elliptic edge points to 
a region of hyperbolic edge points. A doubly covered region on the unit sphere 
appears, and, since (T,),=(Tg)char. =0, the physical streamlines and the family 
of curved physical characteristics return at the image of the line of parabolic 
edge points. According to (108), since t,=1t=-0 and x,=O, this image line is 
everywhere tangent to the direction of the minor principal axis. According to 
(28) a characteristic of the family of straight physical characteristics is also 
tangent to this direction. Since the parabolic edge line does not coincide with 
such a characteristic, the image line is the envelope of the family of straight 
physical characteristics. As can easily be seen, the curved characteristics cannot 
form an envelope. The reversal of the streamline is made possible by a change 
of sign of the velocity. The part of the edge surface downstream of a line of 
parabolic edge points is thus mirrored with respect to the origin. It may be 
noted that as a consequence of the convention adopted for the sign, 0, changes 
sign when mirrored; thus elliptic edge points become hyperbolic edge points 
and vice versa. Along with the reversal of the streamline, the acceleration goes 
to infinity, as may be concluded from (117), and the radius of curvature of the 
spatial streamline goes to zero. It may be noted that on the edge surface the 
geodesic torsion of a curve tangent to the major principal direction is a pro- 
perty which has an equivalent significance as the direction of a curve at points 
on the hodograph surface for finite K,. It is seen from equations (103) and (411) 
that only when 9,—>co are the geodesic torsions of the hodograph streamline 
and of the characteristic tangent to it equal to each other. This is equivalent 
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to the tangency of the hodograph streamline and characteristic at an ordinary 
parabolic point. If t=-0, the line of parabolic edge points is seen to be the 
image of a conical limit line of the first type in a conical simple wave flow. Since 
the image of the line of parabolic edge points is tangent to a straight charac- 
teristic, as is the case for ordinary parabolic points, the following situations 
may prevail. In an expanding simple wave flow the conical limit line is inclined 
along the upstream straight characteristic direction on the concave side of the 
streamline or along the downstream straight characteristic direction on the 
convex side of the streamline. In a compressing simple wave flow the conical 
limit line is inclined along the downstream straight characteristic direction on 
the concave side of the streamline or along the upstream straight characteristic 
on the convex side of the streamline. From equation (67) it is seen that the 
same conclusions with regard to the variation of M, along the streamline hold 
as for an ordinary conical limit line of the first type. 

Singular points on a conical limit line in a conical simple wave flow arise 
at those points where t becomes zero, and singularities similar to those for 
ordinary conical limit lines may be found. 

Of special interest again is the situation where t vanishes all along the space 
curve, which thus is a plane curve. The line of parabolic edge points then co- 
incides with a streamline along which all normals to the edge surface fall in the — 
direction of the binormal to the space curve. A point on the unit sphere is thus 
mapped onto a line on the edge surface. From equation (117) it follows that 
at this point the acceleration along the streamline becomes infinite. Also the 
acceleration normal to the streamline becomes infinite, and the radius of curvature 
of the physical streamline becomes zero. The physical streamline has a corner 
at the image point of the line of parabolic edge points, and the flow may be 
recognized to be a centered simple wave flow. The image point is the center 
of the straight characteristics. Such a singularity may be met, for example, 
in the flow around a supersonic leading edge of a flat delta wing at incidence 
in a homogeneous supersonic flow. 


A conical simple wave flow may be bounded by an analytically different 
region along a characteristic, which, incidentally, may be a conical limit line of 
the second type. 


Along a straight physical characteristic the simple wave flow may be bounded 
by a conical-supersonic parallel flow, since also in such a flow the characteristics 
are straight. The edge surface in the hodograph space then ends at a (conical) 
point of the hodograph surface. We discuss parallel flow in the next section. 


If the conical simple wave flow is bounded by a curved characteristic, the 
edge surface is along a conical characteristic connected with a hodograph surface, 
which may have elliptic, parabolic or hyperbolic points along this line. Stream- 
lines may enter or leave the edge surface from or to the hodograph surface, 
corresponding to flow into or out of a region of simple wave flow. They may 
be obtained from each other by flow reversal. Consider now a streamline at 
an elliptic point on the hodograph surface going onto an edge surface, and for 
example let the flow be as given in Fig. 6a. The streamline on the hodograph 
surface has a direction given by B=0 and deflects to the characteristic angle 
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arc tan + M?—1 when entering the edge surface. Since g,>0, from (60) follows 
0,>0 and g,>0 for an elliptic point; thus (x,)enar > 0, and at the point of the 
edge surface g,>0. The acceleration in the conical simple wave flow is thus 
positive, which is in accordance with the direction of the hodograph streamline 
when the edge surface coincides with a downstream hodograph characteristic. 
If the edge surface is taken along an upstream part of such a characteristic, 
the flow should decelerate; thus g,<0, and 9,<0. This conclusion may be 
obtained by mirroring the edge surface with respect to the origin. The edge 
surface point thus becomes a hyperbolic edge surface point, and a conical limit 
line of the second type occurs. By inspection of Figs. 5 and 6 and by similar 
arguments involving the signs of 0,, @, and g, it may be shown more generally 
that the following rule holds. 


If the streamline at an elliptic or hyperbolic point on the hodograph surface 
enters a region of conical simple wave flow, the edge surface of which coincides 
with a downstream hodograph characteristic, the transition from one type of 
flow to the other will be continuous, while if this surface coincides with an 
upstream characteristic, a conical limit line of the second type occurs. It may 
easily be verified that for elliptic points the interchange of downstream and 
upstream characteristics involves a change in the sign of the acceleration. Con- 
tinuous transition then occurs if the simple wave flow and the neighboring flow 
are both expanding or compressing flows, while limit lines appear if one of the 
types of flow is expanding and the other compressing. 


8. Regions of parallel flow in a conical flow field 


Another singularity in the transformation occurs when the flow is parallel 
in a region of the conical flow field. Since the flow is conical, this region will 
be a cone, and the shape of this cone and the flow adjacent to it are of interest. 
If the parallel flow is conical-supersonic, it may be bounded by straight conical 
characteristics or the envelope of such conical characteristics, being a circular 
conical-sonic line. If the parallel flow is conical-subsonic, the boundary may be 
obtained by considering the hodograph transformation. 

The region of parallel flow is mapped into the hodograph space at a point, 
and A->oo. The hodograph surface, which is the representation of the adjacent 
flow field, exhibits a point where, depending on the way this point is approached 
along this surface, a different normal is found; this point is a conical point, 
and K,—>co*. The cone of the normals at a conical point is congruent with 
and has the same position in space as the cone bounding the region of parallel 
flow. Now if the conical point is approached on the hodograph surface in a 
direction corresponding to the U direction of an arbitrary normal at this point, 
it is seen that Ry vanishes with respect to Ry, and it follows from the curvature 
relation, equation (36), that M,=1. The tangent cone at the conical point thus 


* An exception occurs if the parallel flow is entirely bounded by straight conical 
characteristics, in which case the normal at the point takes on different values 
because edge surfaces end at the point. The point is then not conical, and M, +1 
in general. An example of such a situation is found in the interaction of conical 
simple waves, as studied by GresE & CoHN [13]. 
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coincides with the characteristic cone of that point in the hodograph space, 
the axis of which is aligned along the velocity vector g and the semivertex angle 
cone or a part of it, aligned along the velocity vector g and with a semi-vertex 
angle equal to w=are tan(M@?—1)~*. In the physical space the boundary is 
thus the upstream or downstream pointing Mach cone of the velocity in the 
parallel flow with apex at the center of the conical field, or a part of such a 
circular cone continued by a plane tangent to it along the generator where the 
two figures meet. Obviously, the planes corresponding to a figure composed 
of straight conical characteristics may also be a boundary of the parallel flow. 
The parallel flow may either be downstream or upstream of this boundary. 

From EvuLrr’s theorem (24) it follows that if the radius of curvature in 
any direction approaches zero, one of the principal radii of curvature must be 
zero; thus 9,->0, and DupIn’s indicatrix degenerates to the point Vial on the 
major principal axis and the point VR] =(. The direction of the major principal 
axis jis along the U axis, since @,/0,=0 and M,=1; thus «,=f=0. As for conical 
simple wave flows it may be shown by (108) and (111) that the conical stream- 
lines map into each other. The same may be shown by equations (108) and (103) 
for the conical characteristics. 

In order to investigate the properties of the flow adjacent to the parallel 
region, it will be assumed that the surface in the neighborhood of the conical 
point may be expressed by a series expansion, the first three terms of which 
can be written as 


w=q+A(P)p+ BY) p*. (118) 


A co-ordinate system is thus used which may be obtained by rotating the u, v 
and w axes until the w axis is along the velocity vector g in the parallel flow 
and wu and v are replaced by w= cosW and v= sin¥. The velocity component 
normal to q is thus denoted by #. In these co-ordinates the differential equation 
(18) reads: 
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Substituting equation (118) into (119) and collecting terms of the same 
order yields the following differential equations for the functions A (WY) and B(W): 


Das By OB 1 w2 


[A+ A”] [4-21 ng) eo, (120) 
and 


[1+ A®(1— M2)] B” — 24 A'(1— M2) B+ 
+ 2[2+ (3.4? + A’?+4+ 2A A (4 — M2] B— (121) 
a : M? A(A + A”) [A?{2 + (y — 1) M4 2] =0, 
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where M is the Mach number in the parallel flow, represented by the point 
w=q, p=0. 


Solution of equation (120) yields two results, namely, 


A(W) =Asin(W—%), (122) 


corresponding to a point on the hodograph surface which has a tangent plane, and 


A(W) =A=+—1 12 
which corresponds to a point on the hodograph surface with a tangent cone. 
As before, we again see that the tangent cone is a circular cone with semi-vertex 
angle arc tan /M? — 1 and with its axis along the velocity vector g of the parallel 
flow. 


For a conical point equation (121) simplifies by use of (123) to 


= ely Aree Ee 
BP) =B ia Weave (124) 


Thus to a second approximation the hodograph surface is axially symmetric 
in the vicinity of the conical point. Correspondingly, the flow adjacent to the 
parallel flow is axi-symmetric in the immediate neighborhood of the Mach cone 
of the parallel flow. An analogous result was obtained by BuLAku [18], [20]. 


Furthermore, it is seen from (124) that a plane through g cuts the hodograph 
surface along a line of intersection which is concave towards the velocity vector 
q for |w|<g and convex towards this axis of the tangent cone for |w|>g. From 
this fact the possible types of conical flow with a region of parallel flow may 
be obtained. They are sketched in Fig. 10. In this figure the possible meridian 
cuts of the hodograph surface are given. The corresponding physical flow patterns 
may be constructed by noting that the normal to the hodograph surface should 
be taken positive in the direction such that passing along the hodograph stream- 
line and the corresponding physical streamline does not lead to contradictory 
results. It is seen that continuous transition from or into a parallel flow is 
possible across a downstream Mach cone if the adjacent flow is expanding (type I) 
and across an upstream Mach cone if the adjacent flow is compressing (type II). 
Transitions from or into a region of parallel flow across a downstream Mach 
cone if the adjacent flow is compressing (type III) and across an upstream Mach 
cone if the adjacent flow is expanding (type IV) are possible only if the hodograph 
surface is mirrored with respect to the origin. Thus in these cases conical 
limit lines of the second type occur. In Fig. 10 flows into a parallel flow region 
are denoted by a and flows out of such a region by 6. Furthermore, it may be 
noticed that the types of flow IIb, Ila, IVb and IVa may be obtained from 
the types Ia, Ib, IIIa and IIIb, respectively, by flow reversal. 


Further properties of the adjacent flow may be derived from the actual 
value of B, which yields for the major principal radius of curvature at a conical 
Arch, Rational Mech. Anal., Vol. 6 24 
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point the formula 


eee er (125) 


Mat! 


where the positive sign should be taken in expanding flow and the negative 
sign in compressing flow. From the approximate axial symmetry of the hodograph 
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Fig. 10a and b. Types of conical flow with a region of parallel flow. a) Continuous transition from or into a region of 
parallel flow, b) transition from or into a region of parallel flow across a limit cone. 


surface in the neighborhood of the conical point the sign of 0, may be obtained, 
and so also the sign of the Gaussian curvature. Since 6=0 at the conical point, 
an elliptic point corresponds to conical-supersonic flow, and a hyperbolic point 
to conical-subsonic flow. Inspection of the meridian curves, given in Fig. 10, 
then shows that for a continuous transition the flow goes through M,=1, while 
across a limit cone the flow remains conical-subsonic or conical-supersonic. 


This result may also be obtained by considering the variation of M, along 
the streamline. In order to use equation (67) the acceleration along the stream- 
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line may be found from (60) and (125) to be 


a? VM2—1 


8+ Gtihr M” 


(126) 


and, incidentally, the acceleration normal to the streamline may be obtained from 
equations (62) and (125): 
a M?—1 


See (y+1) 7 M 2 (127) 
while the radius of curvature of the physical streamline is equal to 
M2 
ey te (128) 
The variation of M, along the streamline is then equal to 
aM \ 6 oat 
( al 20 (129) 


which is exactly the negative of the value found at a parabolic point for M,=1, 
equation (91). Thus, in contrast to the case for parabolic points, where the local 
flatness of the hodograph surface leads to flow patterns with a close resemblance 
to two-dimensional plane flows, the infinite Gaussian curvature here leads to 
opposite tendencies. In fact, (129) shows that M, decreases in an expanding flow 
and increasés in a compressing flow *. 


Since M,=1 at the Mach cone, the same conclusions as given above with 
regard to the change of M, across this cone are then obtained. 


The foregoing discussion permits some conclusions with regard to the existence 
of a second conical-sonic line, as discussed first by FEeRrR1 [17] and later by 
BULAKH [22]. 


In [17] the supersonic flow around a triangular conical wing with sharp 
supersonic leading edges and a flat outboard region was discussed. It was shown 
by means of a series expansion in the vicinity of the conical-sonic line and by 
numerical calculations that adjacent to the conical-supersonic parallel flow in 
the outboard region another conical-supersonic region existed, terminated by 
a second conical-sonic line. Since it may be expected that in this case the flow 
will be expanding behind the downstream Mach cone considered, type Ib of the 
flows just discussed will prevail, giving a conical-subsonic flow in the adjacent 
region. As was shown by BULAKH [22], the series expansion used in [/7] was 
not valid in this case, so that conclusions can not be based on it. Also, BULAKH 
showed that in the rotational part of the flow field the conical streamlines are 
characteristics along which discontinuities in the derivatives of the velocity are 
possible. Such discontinuities actually occur along the conical streamline which 


* At elliptic and hyperbolic points the Gaussian curvature is finite and different 
from zero, and both tendencies occur simultaneously, which is expressed by the fact 
that the sign of the variation of 7, is not determined by the sign of the acceleration 
alone. 
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separates the rotational from the irrotational flow, and along which the jump 
of vorticity, introduced by the discontinuity in the curvature of the shock wave, 
is carried downstream. Most likely these discontinuities were not accounted for 
in the numerical calculations. It may also be noticed that in several of the flow 
patterns given in [J7] the conical-sonic lines on the convex body surfaces are 
inclined downstream instead of upstream as would follow from the previous 
discussion on conical-sonic lines. 

If the transition from or into the parallel flow is continuous, the foregoing 
considerations, resulting in flows sketched in Fig. 10, show that a conical- 
supersonic flow adjacent to a parallel flow is possible in an expanding flow in 
front of a downstream Mach cone (type Ia) and in a compressing flow behind 
an upstream Mach cone (type IIb). Actually the latter type of flow may be 
found in the well known Busemann diffuser flow [11]. These flows will in general 
be bounded by another conical-sonic line. 

If the transition occurs across a limit cone, Fig. 10 reveals that a conical- 
supersonic flow adjacent to a parallel flow occurs in compressing flow if the 
Mach (limit) cone is downstream with respect to the parallel flow (type IIIb) 
and in expanding flow if the Mach (limit) cone is upstream with respect to the 
parallel flow (type IVa). 


The foregoing discussion does not pretend to give a complete analysis of the 
singularities in the transformation. For example, the properties of the trans- 
formation if a->0, both when U +0 and when U--0, deserve further attention. 
Also, the existence of branch type singularities, where possibly the hodograph 
surface forms an edge line along which two sheets on the same side of this line 
meet, may be investigated. 

It is hoped, however, that the consideration of the local properties of the 
hodograph transformation for irrotational conical flow from the point of view 
of differential geometry may have proved to give valuable information on the 
structure of conical flow. Nevertheless, the solution of a particular flow problem 
can be completed only by numerical calculation. Such a calculation may be 
performed either in the physical space, using a numerical method as described 
in the references cited above, or in the hodograph space. It would be useful 
to construct both the flow field and the hodograph surface while performing 
these calculations. 

As an example of how the qualitative nature of the flow may be determined 
by means of the hodograph transformation, the flow around a flat delta wing 
with supersonic leading edges will now be discussed. 


9. Supersonic flow around a flat delta wing with supersonic leading edges 


An example of the flow around a delta wing with supersonic leading edges was 
first computed in the non-linear theory by MAsLEN [14] and since then has been 
given special attention by Fowe tt [74] and BuLAKH [21]. Recently, BRoox [29] 
also gave a discussion of the flow field. Consider a flat delta wing with semi-vertex 
angle 6, placed at an angle of attack «% and zero yaw in a uniform supersonic 
stream with a Mach number M. The leading edges will be chosen to be super- 
sonic; that is, 6>arc tan /M?—1 for « not too large. Take the origin of the 
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right-handed co-ordinate system at the center of the conical flow, which is the 
apex of the wing, and the x and y axes in the plane of the wing, the x axis along 
the center line of the wing, positive in the downstream direction, and the y axis 
perpendicular to it, positive on the starboard side of the wing. The z axis is 
normal to the wing surface, and the angle of attack « is chosen positive if the 
positive z axis is on the upper side of the wing. A sketch of the flow field and the 
hodograph surface is given in Fig. 11. All curves and points on the unit sphere 
are projected centrally from the center of the conical flow onto a plane normal 
to the center line of the wing. The flow corresponds approximately to «=5°, 
0=45° and M=3. 

Since the leading edges are supersonic, the flow on the upper and lower 
surface of the wing do not interfere and may thus be treated independently. 

Consider first the upper side of the wing, where the flow is expanding around 
the leading edges. As sketched in Fig. 1 the conical characteristics are straight 
in the parallel flow. The region of flow influenced by the wing is bounded by 
the downstream part B’C’ of the conical characteristic through the point of 
intersection of the leading edge with the unit sphere and by the remaining part 
of the conical-sonic line, being the envelope of the straight conical characteristics 
in the parallel flow (which is not a characteristic curve as stated in [21]). The 
parallel flow is given on the hodograph surface by the point A with the velocity 
vector g=0.802*, corresponding to M=3, pointing upwards at an angle of 5° 
with the u, v plane. The cone of the normals to the hodograph surface in 2 
is also sketched in Fig. 11. Adjacent to the straight characteristic B’C’ the flow 
will be a simple wave flow, and since at C’ the conical physical streamline exhibits 
a corner, the flow is centered in C’. In order for the streamlines in the parallel 
flow to turn around the leading edge they should have their concave side towards 
the wing surface, so that the straight characteristics are inclined downstream 
on the convex side of the streamlines; the flow is thus expanding. The flow 
expands until the velocity vector is parallel to the wing surface and another 
region of parallel flow in between the last characteristic of the simple wave 
flow C’F’ and the wing surface occurs. This parallel flow region maps onto the 
point P, in the hodograph space. The conical simple wave flow is given in the 
hodograph space by an edge surface, lying in the plane through BR, normal to 
the leading edge of the wing O’C’, and extends from R, to R. The curve coinciding 
with the edge surface is thus a plane curve, the binormal of which at every point 
is parallel to O’C’. The line of parabolic edge points corresponds to conditions 
in the center C’ of the expansion waves, where, as was discussed before, the 
acceleration is infinite and the radius of curvature of the streamline is zero. 
Since the flow is expanding, the edge surface contains only elliptic edge points. 
The radius of curvature of PB, may be obtained from equation (115) when 
the velocity is resolved into its components along the tangent, normal and 
binormal of the space curve RP. Since PP, is a characteristic, the velocity 
component gq, along it is sonic. From (2) it follows, when y=1.4, that 


g=a=T (1 — g*) =0.2 — 0.29%. (130) 


* Throughout this discussion the maximum speed will be taken as unity. 
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Furthermore, the velocity component along the binormal g, may be found by 
projection of the velocity in B.. Thus 


gi = (0.802 - cos 5°- cos 45°)? = 0.3192. (131) 


For the velocity component along the principal normal g, then follows 


qu = — Ve — GF — = — 1.29? — 0.5192, (132) 


and by (4115) the radius of curvature of P,P, is given by 


0 = 0.834 /1.2g? — 0.5192. (133) 


Starting at P, in the direction normal to the plane 0’ B’C’, the curve B PB, may 
now be constructed*. The point B, then appears to have a velocity vector 
g=0.810, corresponding to M=3-10 and deflected towards the center line of 
the wing through an angle of approximately 3°. From (117) it may be seen 
that along the first characteristic B’C’ the acceleration increases in the direction 
from Bb’ to C’ with the inverse of sin? to become infinite at C’. In B’, M,=1; 
thus U=a=q,, and W=|/q?—U?=Jq?+¢?, from which it follows that 
Qn (sin 3) += W=a\VM*—1. The acceleration along the streamline in B’ then 
is shown from (117) to be equal to 


DOR VM2—1 
ViteAy ae, MM 


Bs = 0.0565 -- (134) 
This is exactly twice the value of the acceleration in the flow adjacent to the 
conical-sonic line A’ B’, which may be seen from (126) to be 

a? VuMe—1 


8s“ opir M = 0.0283 —. (135) 


The difference in nature of the simple wave flow, determined completely 
by the leading edge having a two-dimensional character, and the flow adjacent 
to the Mach cone of the undisturbed stream, influenced by the apex of the wing 
and having a three-dimensional character, is thus expressed by a jump in the 
acceleration in B’. It may be recalled that in the linear theory there is a dis- 
continuity in the velocity across B’C’, while along A’ B’ the velocity gradient 
in the bordering flow becomes infinite [6]. 


The influence of the apex of the wing, noticeable along A’ B’, enters the 
simple wave flow at the point B’ and carries downstream along the downstream 
conical characteristic B’F’. This boundary was also taken in the solution of 
MASLEN [/4] and given special attention by BULAKH [2/] when commenting on 
a paper by FoweEtt [15]. It was pointed out by FoweE tr [/5] that it should 


* It should be noted that P P, is not the Prandtl-Meyer epicycloid shrunk by 
the factor //1— qj, as is stated in [73] in relation to centered simple wave flow around 
a swept-forward leading edge. 
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also be possible to construct the region of influence of the apex of the wing in 
the spatial flow field by considering the characteristic surface which emanates 
from the apex and intersects the unit sphere along a conical-sonic line. This 
line was then taken as boundary of the region of influence. As was noted in 
Section 2, however, disturbances originating at the apex of the wing admittedly 
travel initially along this characteristic surface but then continue along the 
characteristic surfaces from points on it. The envelope of these surfaces intersects 
the unit sphere along the above-mentioned conical characteristic B’F’*. In the 
hodograph space the shape of this characteristic is obtained from the known 
shape of PB P, and the normal at P, corresponding to B’ (direction O’ B'). With 
the aid of (103) the normals to the edge surface along the image of characteristic 
B’F’ may then be determined, starting from the known values of M, and (%,)char, 
at PR. The boundary between the parallel flow in the outboard region of the 
wing and the inboard region is the continuation of the characteristic B’F’ down- 
stream along the straight characteristics F’ E’. This boundary was also given by 
Masten [14] and Buraxu [21]. It was pointed out by Buraku [2/] that in 
addition to this line the conical-sonic line E’D’ to which F’ E’ is tangent con- 
cludes this boundary line, because when the angle of attack is decreased, the 
inboard region would otherwise fill the whole space influenced by the total wing. 
This argument, however, applies only to small angles of attack, for which a 
conical-sonic line E’D’ indeed may be expected, but the example given by 
MASLEN and several numerical calculations presented by Brook [29] show that 
at higher angles of attack the straight characteristic continues up to the wing 
surface **. 


It is natural to assume that the flow in the inboard region in the vincinity 
of the downstream Mach cone A’ 5’ is expanding, and the transition across 
the conical-sonic line will then be continuous since flow of the type Ib occurs. 
Thus, there is no ground to expect a shock wave to be formed near this line, 
as was done by BuLAKkH [21]. Along A’B’ the flow becomes conical-subsonic, 
and on the hodograph surface points close the conical point P, are hyperbolic. | 
Because of the symmetry of the flow with respect to the line A’G’, the stream- 
line A’G’ will not be curved out of the plane of symmetry; thus B=a,=0. On 
the image of the streamline A’G’ on the hodograph surface RG there will be 
hyperbolic points, where conditions as sketched in Fig. 3a prevail. The point 
G' is a conical stagnation point, where all the streamlines on the upper side of 
the wing come together, and it is mapped onto the orthogonal hyperbolic point 
G in the u,v plane. In the vincinity of the other downstream Mach cone D’ E’ 
the flow may be assumed to be compressing. This may be justified by noting 
that in the parallel flow along the wing surface the streamlines are directed 
towards the center of the wing and have to be deflected to become parallel to it. 
Thus the streamlines must be curved with the concave side to the leading edge, 


* Another conical flow solution where the region of influence of the origin has 
been incorrectly chosen was given in [30], where the supersonic flow near the junction 
of two wedges was studied. 

** In the case that a conical-sonic line D’E’ appears, in [29] the characteristic 


tangent to D’E’ in D’ is incorrectly taken as the boundary between the regions 
determined by the leading edge and the wing apex. 
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and the pressure gradient is therefore opposite to the direction of flow (and 
lying in the x,y plane, normal to the radius). Stated more freely, the two- 
dimensional expansion around the leading edge would be too severe for a three- 
dimensional flow, and as soon as the influence zone of the apex is reached, it 
appears that an overexpansion has taken place which has to be neutralized by 
a compression. From the considerations given in the section on parallel flow 
it follows, however, that the transition to compressing flow along a downstream 
Mach sone leads to flow type IJIb, such that the Mach cone is a limit cone. 
The curve D’E’ is therefore, for all angles of attack, a conical limit line. The 
flow near D’E’ in the inboard region is conical-supersonic. Since the streamline 
D'G' along the wing surface is not curved out of the x, y plane, it follows that 
B=0, and on the hodograph surface D’ H’ is mapped onto a curve of elliptic 
points FP, H lying in the wu, v plane where conditions as sketched in Fig. 6d prevail. 
In order for the streamline to arrive at G’ it should return once more, while since 
M,=0 in G’, the flow should change from conical-supersonic to conical subsonic. 
Since B=0 all along the streamline on the wing surface, the latter means also 
that on the hodograph surface at the same point the points along the hodograph 
streamline change from elliptic to hyperbolic. This transition occurs at the 
parabolic point H, where since B=0, M,=1 (Fig. 7d), while the image point 
H’ is a point of a conical limit line of the first type. The streamline H’G’ is 
mapped onto the line of hyperbolic points HG in the wu, v plane, where conditions 
exist as given in Fig. 3e. Along P,GBP, the hodograph surface is normal to the 
u,v plane. At H the parabolic line is normal to the uw, v plane, while at H’ the 
conical limit line is normal to the x, y plane. On the hodograph surface the 
parabolic line runs to some point K, and the whole of this surface containing 
elliptic points (excluding the elliptic edge points on P, 2) is mirrored with respect 
to the origin. In Fig. 14 the area to be reflected is enclosed by a dotted line. 
Evidently, the disturbance introduced at B’ propagates along the downstream 
characteristic B’E’ and the conical sonic line E’ D’ in such a way that down- 
stream of some point K’' a conical limit line of the second type is formed. The 
exact location of K’ can be found only by a numerical calculation of the whole 
flow field. In fact, the hodograph surface, excluding the edge surfaces BP, 
and P,P; only, contains regions where the differential equation is either elliptic, 
parabolic or hyperbolic, depending on an elliptic region.' For reasons which will 
become apparent later, it will be assumed that K’ lies on B’F’ somewhere in 
between B’ and F’. The straight characteristic F’E’ is then a limit line. The 
flow adjacent to it must have straight characteristics and thus will be a conical 
simple wave flow. Furthermore, the waves will not be centered, because other- 
wise the flow would be completely determined by the outboard region. In the 
hodograph space this conical simple wave flow is mapped onto an edge surface 
PL, which coincides with a space curve, since its torsion is different from zero. 
In P, the tangent to BL is normal to the plane O’E’F’, and since the flow is 
decelerating, the edge surface consists of hyperbolic edge points, which are 
obtained by mirroring the edge surface containing elliptic edge points. The 
region of conical simple wave flow will be bounded by the curved downstream 
characteristics from E’ and F’, which arrive at their point of intersection L’ 
as upstream characteristics. Since both characteristics are curved and only two 
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characteristic directions are possible, the direction of the straight characteristic 
must coincide with at least one of them at L’; thus M,=1, and both curved 
characteristics are tangent to the straight characteristic. The edge surface PL 
therefore ends on the conical-sonic line which runs from H to the conical point F.. 
Along F’L' and E’L' propagates a discontinuity in the acceleration, caused by 
the junction of a conical simple wave 
flow and an analytically different re- 
gion. It will be assumed that along the 
limit line H’ K’, H’ is the only conical- 
sonic point. The point L’ then lies 
downstream of H'K’, and the charac- 
teristic F’ L’ and E£’ L’ intersect the limit 
line in J’ and I’, respectively. Since the 
transition across F’J’ occurs continu- 
ously, the image of the region F’ J’ K’ on 
the hodograph surface P, J K consists of 
elliptic points, as do the edge surface 
in the region BJ J and the hodograph 
surface in the region P,J H, which is 
the image of region D’E'I'H’ on the 
unit sphere. The parabolic line Hk 
enters the edge surface at J and con- 
tinues as a line of parabolic edge points 
I J, leaving this surface at J to become 
a parabolic line J/ K again. The image 
of IJ on the unit sphere J’J’ is the 
envelope of straight characteristics. 
if Conditions in the inboard region 

d along B’F’ may be deduced from the 
hodograph surface by assuming that B 
varies continuously along the charac- 
teristic separating the edge surface PB, 
from the hodograph surface. Typical 
points along this curve are chosen in 
eage surface Fig. 12 to show the different situations 
Fig. 12a—e. Typical conditions on the characteristic B’F’ which prevail. The image of the stream- 
line through 5’ on the side of F’ leaves 

the edge surface P, P, in PB along the U axis, which is also the direction of Bit; 
at this point, and B=o,=0, where M,=1. From B’ to F’ along B’F’, M, increases; 
thus on the edge surface | B| increases when going from FP, to P;. On the hodograph 
surface |8| increases more rapidly, in accordance with the fact that at hyperbolic 
points the hodograph streamline should lie in the region which the hodograph 
characteristics enclose around a direction normal to the U axis. It follows then 
that at these points the flow is still expanding if |B] < 4a, as sketched in Fig. 12a, 
and when || is increased more, the acceleration along the streamline becomes 
Zero for || =32 (Fig.12b), while the flow is compressing if |B| exceeds this 
value (Fig.12c). In all these cases the transition across B’F’ is continuous, 
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since if the flow is reversed, the streamline on the edge surface is seen to be 
directed along the downstream characteristic*. Thus in the neighborhood of B’ 
in addition to the small two-dimensional expansion a spatial expansion occurs, 
and further downstream along B’F’, a continuous compression. If B’F’ is 
followed further downstream, || changes enough to bring the streamline into 
the other characteristic direction when the parabolic point K is reached. The 
deceleration has then become infinite, and the situation is sketched in Fig. 12d. 
It may easily be shown, from the transformation, that the conical limit line of 
the first type in K’ touches the characteristic B’F’, which from this point on 
becomes a conical limit line of the second type. Still further downstream the 
hodograph streamline enters the region which the characteristics enclose around 
the U axis, and the point is thus elliptic, which also may be deduced from the 
fact that the parabolic line is crossed. The transition is now across a limit line, 
since if the flow is reversed, the streamline on the edge surface is directed along 
the upstream characteristic, while the flow changes from expanding to compressing. 
This situation is sketched in Fig. 12e. The course of the physical conical stream- 
lines and the hodograph streamlines on the expansion side of the wing may 
now easily be understood from the considerations given above. 


Because of the appearance of limit cones at all angles of attack it may be 
concluded that no continuous solution for the flow on the expansion side exists **. 
In order that a physically possible flow pattern can occur, in the region of the 
limit cones a shock wave must be formed. The formation of this inboard shock 
may be expected to start with zero strength at the point K’ tangent to the 
two limit cones and to continue until it hits the wing surface normally ***. 


Indeed, as was shown by FoweEtt [15] in experiments, an inboard shock 
wave was found. The flow downstream of this shock wave does not necessarily 
have to be conical-subsonic, and it will be rotational if the shock is curved. 
In G’ a Ferri singularity in the entropy distribution will then occur. 


On the lower side of the wing a compression takes place, and an isentropic 
solution involving limit cones may also be constructed. Yet, obviously, unless 
detachment occurs, a plane shock wave C’M’ may be seen to be attached to 
the leading edge, downstream of which a region of flow parallel to the wing 
surface exists. This region of parallel flow is bounded downstream by a Mach 
cone, showing up as a circular conical-sonic line M’ P’ on the unit sphere. Since 


*x This continuous transition may be shown not to be in contradiction with the 
fact that the characteristic separates regions of elliptic edge points and hyperbolic 
points, because the edge surface lies on the same side of this separation curve as 
the hodograph surface. 

*x This conclusion was also reached in [75] and [29], where it was remarked that 
the simple wave regions from the leading edges at some angle of attack overlap, 
thus making a continuous solution impossible. As was pointed out before, however, 
this picture of the flow is based on an incorrect presentation of the spatial influence 
of the apex of the wing and does not justify conclusions with regard to the possibility 
of a continuous solution. 

x*xk Tn [21] and [29] it was assumed that the shock wave begins at D’. There is 
no reason, however, to expect that the shock will extend beyond K’ in the region 
where a continuous transition across B’F’ was found. 
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in the parallel flow the streamlines are directed away from the wing center 
line, an expansion occurs along O’ P’ in order to turn the streamlines towards 
the direction of the center line. The type of flow Ib therefore prevails, giving 
a continuous transition across C’M’ into a conical-subsonic flow. Thus, there 
is no reason to expect a shock wave to be formed in the neighborhood of C ‘M' 
as was done by BurakH [27]. In the example given by MASLEN [/4], along 
C’M’ the flow was found to be conical-supersonic. As was pointed out by 
BuLAKH [22], however, the jump in curvature of the shock wave at the junction 
of the plane shock C’M’ and the curved shock M’N’ introduces a jump in vorti- 
city, propagating downstream along the streamline M’Q’, which is most likely 
not accounted for in the numerical calculations in [14]. Since downstream of 
the curved shock the flow is rotational, a velocity potential cannot be defined, 
and the hodograph transformation breaks down. Also for this reason, the hodo- 
graph surface for the compression side of the wing is not given in Fig. 11. Again 
a Ferri singularity in the entropy distribution occurs at Q’. 


This completes the present description of the supersonic flow around the 
delta wing with supersonic leading edges. 


A note may be added concerning the supersonic flow around the tip of a 
flat rectangular wing at incidence, which was discussed in connection with the 
delta wing by BuLaxu [2/]. Since this flow pattern essentially may be obtained 
by turning one of the leading edges of a delta wing with 90° apex angle perpen- 
dicular to the undisturbed flow while the other then becomes a subsonic edge, 
the flow on the side of the leading edge may be expected to be similar to that 
for the delta wing as discussed before. In fact, it may be shown that on the 
expansion side two limit lines again appear, which are tangent at their point 
of contact on the characteristic terminating the Prandtl-Meyer expansion, and 
thus a continuous solution is not possible. Again an inboard shock will be formed, 
starting with zero strength at the point of tangency of the limit lines. On 
the compression side of the wing the parallel flow is again bounded by a down- 
stream Mach cone along which a continuous transition into expanding flow 
takes place and no shock occurs. The hodograph surface for this flow pattern 
as given by BUSEMANN [23] may be modified, according to the results found 
for the delta wing, to include regions of elliptic points, in which an edge surface 
is imbedded, and separated from the region of hyperbolic points by a parabolic 
line. The conditions on the subsonic edge are given by a parabolic line, coinciding 
with a characteristic, which also is a streamline. 
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